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ABSTRACT 

This study investigates the maneuvering of ships in- 
volved in the replenishment at sea operation under calm 
water conditions. 

Two sets of linear differential equations of motion 
of a ship in three degrees of freedom are implemented for 
an analog— digital simulation (hybrid operation) . Mainly 
a two phase hybrid simulation is carried out. 

In the first part the real-time dynamic response of 
a single ship is obtained. Small perturbations, result- 
ing from small changes in rudder angle and propeller RPM, 
are studied. 

In the second part interaction forces and moments 
are applied only to the leading ship as it is overtaken 
by the tracking ship. Thus the response of the leading 
ship is of primary interest on this phase. 

Finally interaction forces and moments are applied 
to both leading and tracking ship. Thus the complete 
response of the UNREP operation is obtained. 
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= Stability criterion for stability in 
straight line motion 


I ,I ,I 
X y z 


= Moments of inertia about the X, Y, Z axes, 
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= Rolling, pitching, and yawing moments, 
respectively 
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= Ship length between perpendiculars (LBP) 
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= Mass of ship 


p^q/r 


= Angular velocities of roll, pitch, and yaw, 
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= Initial equilibrium velocity component 
(ahead straight line motion at constant 
speed with rudder at amidships) 


Au 


= u •- U2 


• • • 

U, V, w 


= Acceleration components of the origin of 

the body axes relative to the fluid (longi- 
tudinal, transverse, and normal components, 
respectively 


X,Y,Z 


= Hydrodynamic force components in ship body 

(longitudinal, lateral and normal components) 


x,y,z 


= Coordinate axes fixed in ship. Origin of 
axes system need not be at the center of 
gravity of the ship (positive direction 
forward, starboard, and downward, respec- 
tively 
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Coordinates of the center of mass of the 
ship relative to body axes 
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II 


Coordinates relative to the coordinate 
system fixed in the earth 
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°G °G °G 


Coordinates of the center of mass of the 
ship relative to the coordinate system 
fixed in the earth 
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Angle of drift. Lateral side to side 
separation distance betv;een ships 
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Angular displacement of a control surface, 
normal to the rudder angle 
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P 


Mass density (mass density if sea water is 
1.9905 lb— sec 
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0 
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Velocity components of the ship axes 
relative to the space coordinate system 


a = 


Longitudinal separation distance between 
midships 


An 


Small change in RPM propeller speed 


AR 


Small change in rudder angle 
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I. 



INTRODUCTION 



At the present time the operational procedure of re- 
plenishment at sea underway is described by two manuals 
of the Department of the Navy [1] and [2] . This procedure 
includes collision hazard. Therefore an analysis of the 
replenishment operations may provide a mean for improve- 
ment of ship maneuvering and ship control. The maneuver 
of replenishment at sea involves six factors: course, 

speed, distance between ships, the approach, station keep- 
ing and departure. 

This thesis describes the analysis of underway replen- 
ishment (UNREP) using an analog computer simulation. This 
analog computer simulation model will describe two Mariner 
Class merchant ships. An effort is made to define the 
control parameters necessary to provide information for 
the conning officer and/or helmsman v/hich will increase 
the safety conditions of the UNREP operation. According 
to the tactical requirements the UNREP operation is inter- 
preted as follows: 

a. The replenishing (leading) ship is responsible 
for course keeping only. 

b. The receiving (tracking) ship is responsible for 
both course and station (distance between ships) 
keeping. 

For this reason one ship is constrained to move in a 
straight line course at constant speed. The rudder and 
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propeller of the second ship are controlled in order to 
keep it in a parallel course relative to the first ship. 

The analog computer output is presented and analyzed 
to determine the necessary information for the ship con- 
trol parameters required during these simulated UNREP 
operations . 
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INTERACTION EFFECTS 



II . 

When underway, the pressure distribution on the ship 
hull varies (as is shown in Figure 1) due to the variable 
velocity of the water flow, the so-called venturi effect. 







Jucr/oj\' 

/■' 'jy' -' 

' I \ 

Figure 1. Pressure distribution on the hull. 

When two ships are alongside underway, this venturi 
effect is increased and becomes more complicated due to 
the mixing of the pressure areas of two ships. 

Figure 2 shows ships which are in dangerous position 
because they are being acted on by radically different 
pressures. Reduction of speed increases the danger of 
collision. Also the pressure effects are more exaggera- 
ted and extra care is required in maneuvering (in depths 
less than 20 fathoms) . 

It is understood that to maintain station during 
UNREP operation a certain amount of rudder is required 



■~suc.r/oA’ 
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which depends on size, load, sea conditions, speed and 
ship separation. But this change of rudder deflection 








^ocr/oAi' 







Figure 2. Ships alongside underway. 



decreases handling capabilities of the receiving ship. 

The problem of interaction effects has been studied 
theoretically by Silverstein [3] and experimentally by 
Newton [4] with both approaches showing agreement in the 
general conclusions. 

When two ships are underway with a certain separation 
and on parallel courses, the pressure fields mix which 
results in an unbalanced force and moment on each ship. 
This interaction moment must be cancelled by the rudder 
action in order to maintain station. However, there is 
a position where the rudder force tends to add to the 
force of attraction, see Figure 3. Thus in these posi- 
tions a bigger rudder deflection is required so the yaw 
angle produced in this way creates an outboard force 
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Figure 3. Relative position where both interaction 
forces and moments tend to draw one ship 
toward the other. 



which balances both the rudder and interaction 
force . 

Newton's experiment, both with models and full scale 
trials in open sea, prove one important thing; It is the 
process of approach or departure from the abeam or "fuel- 
ing" position that includes the maneuvering risks. It 
should be noted here that an increase in velocity or de- 
crease of separation between ships results in the increase 
of the interaction effects. 
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III. 



DERIVATION OF THE LINEAR MATHEMATICAL MODEL 



A. GENERAL CASE-SIX DEGREES OF FREEDOM 

It is well understood that the coordinate— axes system 
fixed on the ship's body is an Eulerian system. It is as- 
sumed that at t=0 the Eulerian system and the space co- 
ordinates system coincide, as shown in Figure 4. The geo— 




Figure 4. Orientation of fixed axes and moving axes. 



graphic coordinate system then describes the motion of the 
ship through the fluid medium in six degrees of freedom. 
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v> 



The Newton's laws of motion in vector form are: 



' > 

F (External Force) = ^ (momentum) 

M (External Moment )= ^ (angular momentum) (III— 1) 

Then the six equations (six degrees of freedom) describing 
the ship's motion have been found [5] 

2 2 

X=m [u+qw-rv-X|^ (r +q ) +y^ (pq-f ) +z^ (pr+q) ] 

2 2 

Y=m [v+ru-pw-y^ (r +p ) (qr-p) +x^ (qp+v) ] 

2 2 

Z=m [w+pv-qu-z^ (p +q ) +x^^ (rp-q) +y ^ (rp+p) ] 

K=plj^+ (l 2 -Iy)qr+m [y^ (w-qu+pv) -z^ (v-pw+rv) ] 

M=qly+ (Ij^-l2)pr+m(Zg(u-rv+qw)-Xg(w-qu+pv) ] 

N=f 1 2 + ( I ^-1 ) pq+m [x g ( v-pw+r u ) -y q ( u-r v+qw ) ] (II 1-2 ) 



Using (III— 2) the following equations can be satisfied: 

^ /N ^ 

F = iX + jY + kZ 



(III-3) 

^ /S * /\ /N 

M = iK + jM + kN 

and where: m — mass of the ship 

X, Y, Z — components of force in the x,y,z directions 

K, M, N — components of applied moment about x,y,z 
axes 

u, V, w — components of velocity in the x,y,z 
directions 

Xq, Zg — distances of origin from center of gravity 
in x,y,z directions 
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p, q, r — components of angular velocity about the 
x,y,z axes 

I ,I ,I — moments of inertia about the x.y.z axes 
X y z ' 

Equation (III— 2) describes the reaction of the rigid 
body of the ship itself as a function of its geometric 
and physical characteristics. Note that they do not in- 
clude any external moments such as due to fins. 



B. HORIZONTAL PLANE MOTION 

It is well understood that the ship's motion in calm 
waters is described only by the following four equations: 

2 2 

X = m[u+qu-rv-x^(r +q ) +y^ (pq-f ) +z^ (pr+q) ] , surge 

2 2 

Y = m[v+ru-pw-y^(r +p ) +z^ (qr-p) +x^ (qp+r ) ] , sway 

N = fl +(I — I ) pq+m [x (v— pw+ru)— y (Cl— rv+qw) ] , yaw 
z y X Cj ^ 

K = pi +(I — I ) qr+m (y^ (w— qu+pv)— z (v— pw+rv) ] , roll (III— 4) 
X z y Cj Cj 

because under the calm water assumption it is true that 

Roll = Pitch = Heave = 0. That is the horizontal motion 

implies: p=p=q=q=w=w=0. Hence the equations 

(III— 4) can be written: 

2 

X = m[u-rv-x^r 

2 

Y = m[v+ur+x^r-y r ] 

Cji Cj 

K = m[-Zg(v+rv) ] 

N = fl +m [x (❖+ru)— y (u— rv) ] (III— 5) 

Z (j Lj 

Equations (III— 5) can be further simplified assuming 
that the center of gravity (CG) is placed at the origin 
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0. 



of the X, y, z coordinate system then x = y = = 

G G G 

Hence, neglecting the roll eguation we have: 

X = m[Ci~rv], surge 

Y = m[v+ur], sway 

N = fl , yaw (III— 6) 

z 

There are the reduced equations for steering and maneu- 
vering of a ship. It is noted that the left hand side of 
equations (III— 6) represents the forces and moment along 
the coordinate axes, and the right hand side shows the cor- 
responding dynamic response terms on the horizontal plane 
of motion. 

C. LINEARIZATION THROUGH TAYLOR'S SERIES EXPANSION [5], [6] 
The forces and moments on the left hand side of equa- 
tions (III— 1) through (III— 6) can be expressed as functions 
of properties of the body, properties of the fluid and 
motion. 

Since steering and maneuvering are of interest, forces 
and moments are also considered as functions of rudder 
(control surface) deflections and the change in r.p.m. 

(An) of the propeller shaft. Thus: 

Forces = f (properties of motion, rudder deflection, etc .) = 

MomentsJ = f (x , y , z , c|) , 0 , ij; ,u , v, w,p , q , r , u , v , w, p , q , f , 

<— > < : > 

orientation motion parameters 

parameters 

6 , 6 , 6 , An, etc . ) 

< > 

control surface parameters (III— 7) 
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For a surface ship moving on the horizontal plane no 
forces or moments are due to orientation changes. Thus 
the forces and moments will be functions of three degrees 
of freedom of motion parameters, rudder deflection para- 
meters and changes in r.p.m. Hence: 



X 

Y 

N 



1 



= f (u, V, r ,u, V, r , 6 , 6 , 6 , An, etc.) 



(III-8) 



In general we linearize a function f(x) by the use of 
Taylor's series expansion. 






Ax^ d”f(x) 



n; 



dx 



n 



where 



Ax = X — X, 



For small values of Ax the second order terms can be 
neglected and thus considering only the following expres- 
sion for f (x) : 
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This is the linearized Taylor's series expansion. The 
same principle can be applied for small perturbations in 
the equation (III— 8) , which is a function of many 
variables . 

Since the Taylor expansion is written for a particu- 
lar point we choose this point to be an equilibrium posi- 
tion. An equilibrium position is that of straight ahead 
motion, at constant speed with rudder amidships. The 
hydrodynamic forces and moment have been found [6] to be: 

X=XAu+Xv+Xr+ X-u + X-v + X'r + X.6 + X An 
u vruvron 

where Au = u — u^ and An = n — n^ and 

the subscript 1 is referred to the values of the variables 
at the initial equilibrium condition and where all the 
partial derivatives are evaluated. For Y and N similar 
expressions hold. Equating the linearized expression for 
X, Y, and N with equation (III— 6) results in the linear- 
ized equations of motion for steering and maneuvering: 



X Au+X v+X r+X-u+X. v+X-r+X.6+X An = mu 
u vruv ron 

Y Au+Y v+Y r+Y-u+Y- v+Y-r+Y^6+Y An = m(v+ru, ) 
u vr uvron 1 

N Au+N v+N r+N-u+N-v+N-r+N^6+N An = I r (III-IO) 

u vruvron z 



Note that the term mvr in the right hand side of the first 
of equations (III— 6) has been dropped since the ship was 
assumed in straight ahead motion. The derivatives 



X , 

V 



X- , 

V 



X , 
r 



X- , 
r 



X, 



Y 



U' 



Y- 

u 
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and 



^u' vanish for any symmetrical port and star- 

board shape of ship (symmetry about the xz— plane) . This 
has the effect of decoupling surge from sway and yaw. 
Thus equations (III— 10) become: 



(X*-m)u+X Au+X An = 0 
u u n 

(Y*— m)v+Y v+ (Y — mu, ) r+Y. f +Y „6+Y An = 0 
V V r 1 r 6 n 

(N--I )f+N v+N.v+N r+N,6+N An = 0 (III-ll) 

r z V V r 0 n 



Using as a basis the linear mathematical model of 
equations (III-ll) the criterion for dynamic stability in 
straight line motion [5] is evaluated as 



C 



Y N - N (Y -mu, ) 
V r V r 1 



> 0 



It is obvious that the treatment of the linear mathema- 
tical model, for the study of a dynamically stable ship, 
can be applied to predict maneuver ing^'and control only for 
small deviations from the original straight line motion 
due to small rudder deflections and small changes in r.p.m. 

D. NON-DIMENSIONAL EQUATIONS 

Table I gives the dimensionalized , nondimensionalized 
quantities and their respective conversion factors. The 
prime notation corresponds to the nondimensionalized quan- 
tities. In equations (III— 10)the force equations are 

2 2 . 3 2 
divided by (p/2)L u^^ and the moment equation by (p/2)L u^^. 
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The resulting nondimensional equations are: 



(X! -m' )u'+X'A’+X'A'n= o 
u u u n 



(Y! — m' ) V ' +Y ' V ' + (Y ' 
V V r 



-ra'u’ )r '+Y'r’+Y'6'+Y’A'n = 0 
1 r 6 n 



(Nl -I’ )r '+N’v'+Nlv’+N'r'+NU '+N' 6 ’n= 0 
r z' V V r 6 n 

It is noted here that 



u 



1 




1.0 



for small perturbations. 



(IV-1) 
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TABLE I 



NONDIMENS lONAL NOMENCLATURE 



Symbol 


Non-Dimensional 

Form 


X- 

u 


X* 

1 _ u 

u 1,2 

J PL 





X 


X 


t 

II 

X 


u 


u 1 ^ z 




2 ^1 



Definition 



Derivative of longitudi- 
nal force component with 
respect to longitudinal 
acceleration component u 



Derivative of longitudi- 
nal force component with 
respect to longitudinal 
velocity component u 



Y 



V 



Y' 

V 



V 



1 

2 f'' “l 



Derivative of lateral 
force component with 
respect to transverse 
velocity component 



Y* 



v 




Derivative of lateral 
force component with res- 
pect to transverse 
acceleration component 



Y 

r 



Y' 

r 



1 

2 



^3 

pL u. 



Derivative of lateral 
force component with res- 
pect to yaw angular 
velocity component 



Y. 



r 




Derivative of lateral 
force component with res- 
pect to yaw angular 
acceleration component 



Y 



6 



Y ' 

.6 



1 ^2 2 
2 “l 



Derivative of lateral 
force component with res- 
pect to rudder angle 
component 
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•Hlrsi I«h|cni IrHiCM Mi(N 



1 

2 



N 



V 



^3 

pL u- 



Derivative of yawing mo- 
ment component with res- 
pect to transverse 
velocity component 




Derivative of yawing mo- 
ment component with res- 
pect to transverse 
velocity component 



— Derivative of yawing com— 

pL ponent with respect to 

yaw angular velocity 
component 



^ Derivative of yav/ing mo— 

p L ment with respect to yaw 

angular acceleration 
component 



Ng 

— 2 — 2 Derivative of yaw moment 

pL Uj^ component with respect 

to rudder angle component 



X 

n 



1 t3 

2 ""l 



Derivative of longitudi- 
nal force component with 
respect to change in 
propeller rpm 



• Y 

2 — Derivative of lateral 

2 pL u, force change in with res- 

pect to propeller rpm 



N 

n 



1 

2 “l 



Derivative of yawing com- 
ponent with respect to 
change in propeller rpm 



rL 

U 




Yawing angular velocity 
component 



Yawing angular accelera- 
tion component 



u 



U' = 1 



V 



V ' = 



V 

U 



vL 



X 



X' = 



X 



1 ^2 2 
2 



Y' = 



1 ^2 2 
2 P^ 



m 



m' = 



m 



i P 



u 



u’ = 



u 

U 



u 



u* = 



uL 

u2 



I’ 

z 



m 



1 t5 

2 P L 



Velocity of origin of 
body axes relative to 
fluid 

Transverse velocity compon- 
ent of origin of ship axes 
relative to fluid 

Transverse acceleration 
component of ship axes 
relative to fluid 



Hydrodynamic longitudinal 
force (positive direction 
forward) 

Hydrodynamic lateral force 
(positive direction to 
starboard) 

Mass of the ship 



Velocity of origin of 
ship's axes along the 
X— axis 



Acceleration of origin of 
ship's axes along the 
X— axis 

Moment of inertia of the 
ship with respect to the 
z— axis 
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IV. AXES FIXED RELATIVE TO THE 
EARTH-LANGRANGIAN SYSTEM 



To convert equations (III— 6) from the Eulerian system 
to the Langrangian system, see Figure 4, the stationary 
coordinate system, the total forces in the x and y direc- 
tion in the ship coordinates must be expressed in terms 
of the stationary coordinates as follows: 

X = X cos ip + Y sin \p 

Y = Y cos ih + X sin ip ' (V— 1) 

o o 

Also the center of gravity possesses velocity and ac- 
celeration with respect to the stationary coordinate sys- 
tem expressed as follows: 

k ■ = u cos — V sin ip 
°G 

y = u sin ip + v cos tp (V— 2) 

°G . 

and 

x = u cos ^ — V sin ip — (u sin + v cos ip)^ 

°G 

•• • 

y = u sin ij; + v cos ip +(u cos — v sin ip^'P 
°G 

(V-3) 

Equations (V- 2) can be used to calculate the trajec- 
tory and velocity of the ship with respect to stationary 
coordinate system. Equations (V- 3) can be used to calcu- 
late the acceleration along this trajectory. 
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V. ANALOG COMPUTER PROGRAMMING FOR 
EACH SHIP INDEPENDENTLY 



A. PRELIMINARY STUDY 

In this section we will discuss the simulation of the 
ship dynamics in open calm sea. The equations of motion 
for each ship are programmed independently and no interac- 
tion forces and moments are considered. 

1 . Linearized Equations of Motion 

Equations (IV— 1) are repeated here for convenience 
as follows: 

(m'-x:)u' = X'u' + X'6n, Surge 
u u n ^ 

(m'-y:)v' = Y'v' +{y;-^'u' )’ l’’+Y.>'+Y'6 '+Y' 6R' , Sway 

V V \l> o ^ n n 6R 

(I'-N.',)ij}' = N'v' + Nli'Nlv'+N' 6R'+N' 6n' , Yaw (V-l) 
z ifi ^ V V 6R n 

where 

r = 
r = ip 

^6 = ^6R 
An = 6n 

u = u, 
o 1 

It is noted here again that all the variables corres- 
pond to small perturbation variables. It is first chosen 
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the equilibrium condition as follows: 

Nominal values: = 15 knots = 25.32 ft/sec = u^ 

V = 0 
o 

= 0 
^o 

2 . Magnitude Scaling of the Linearized Equations 

Following the usual procedure each equation of 
equations (V— 1) is scaled by letting each variable be re- 
presented by its maximum value times the scaled variable. 
That is : 

• • • 

u = u u 
m 

and in non-dimensional form: 



m 



where : 



u^ is the maximum expected value 



u is the scaled variable whose value ranges from 
—1 to +1 computer units (i.e.+ 100 volts for 
Cl— 5000 analog computer) . 

Makint this substitution for all the variables in 
equations (V— 1) and solving for the scaled variables u, v, 



and we obtain: 

X' u' 

u m 



X' 6n’ 
n m 



— • r u. lu T I . r 1 r I 

• = [”7 — i — < rr- \ — Ju* * + [- 7 — i — ■ , V > -rJ 6n' 

n Wm'-yMn' ‘‘ (m ' “X I ) U ' 



u 



(m'-XDu 

urn 



Y' v' 

V m 



u m 



(Yl“-m * u ) ij; * — 

\p o 



i(j 



m 



V m 



V m 



-] 6n ' 



Y' <5R' 

V m 



Y' 6n' 
n m 

(m'-Y:)v* 
V m 
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— N ' V ' N I ' 

j;. = [ !L_S_] V' + [ i--!!::— ] 3T' 



(I'-Nv !{;' 



( I ' -N ; ) ij; ' 



+ [■ 



Nl v' Nl._6R' 

■" 1 V’ + I ? -1 Sr' 



( I ' -N ' ) ' 

z m 



(I ' 

z 



N'6n' 

I r ^ 1 X“ I 

+ [ n— ] <Sn' , 

z i|; m 



(V-2) 



Letting each term in brackets be a coefficient (k) , 
equations (V— 2) can be written as: 



u' 








V' 


= -k^^v' 




' -k3_9<Sn 


TT 






' +k^ T 6n 
115 



(V-3) 



where the appropriate sign has been extracted so that each 
coefficient is now positive. Table IV gives the nondimen— 
sional values necessary for equation (V— 2) and (V— 3) . The 
following data are also given in Table II. 



TABLE II 



= 25.32 ft/sec (nominal ship's velocity) 

u^ = 640.0 (ft/sec) p = 2 lbs sec^/ft^; L = 527.87 ft 

= 9.428E+10; ^pL^u^ = 17.86E+07 
o 2 o 

jpL^= 4.095E+13; ipL^=l . 4707E+07 ; ipL^=2 . 7 862E+05 
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3 . Calculation of Potentiometer Coefficients 

In obtaining a set of scaled equations, it is best 
to work with dimensional quantities. Actually, scaling the 
nondimensional equations is equivalent to scaling the di- 
mensional set. Tables III through V give all the necessary 
data for the scaling procedure. Referring to dimensional 
quantities equations (V— 2) can be rewritten as: 

Xu X (5n 

• r ^ 1 • f ^ 1 I* 

" = ‘ (m-x.)u ‘ (n.-X.)u 

u m urn 



Y V ( Y , —mu ) 

. , V m ^ f if; o y 

U = [7t:=^r-T7r-] V + 



(m-Y. ) V 
V m 



(m-Y . ) 

V m 



^ ^ (m-Y , ) V ^ ^(m-Y.)'Cr ^ 

V m V m 



Y 6n 

V m 



— N V N, if; N. ❖ 

, _ r V m 1— . r m iT . r V m 

^(I ^ (I -N“)ip ^ ^ (I ^ 

' z i|j ^m ' 2 2 ^m 



Ntt, R N n 

I r u Xv m 1 jc I r m -] 

* ‘TT -n-;)? I- (I ' ' 

z Ip z ip ^in 



(V-4) 



Let D = (m— X.)u = mu — X.u 
Hum m urn 

= (1.166 X 10^ slugs x 0.844 ft/sec^ 

+ 1.0295 X 10^ Ibs/ft/sec^ x 0.844 ft/sec^) 

= 0.984 X 10^ + 0.0866 x 10^ = 1.07 x 10^ lbs 



D- „ = (m— Y.)v = mv — Y.v 
12 V m m V m 

= 1.66 x 10^ X 0.844 + 1.043 x 10^ x 0.844 

= 0.984 X 10^ + 0.880 x 10^ = 1.864 x 10^ lbs 
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D^-, = (I -N")i|j = I ip -N-;ip 

13 z t|j z^m 

= 1.43 X 10 ^*^ lbs ft/rad/sec^ x 0.87 x 10 2 

sec 

= 1.244 X 10 ^ + 1.533 x 10 ^ = 2.777 x 10 ^ Ibs-ft 



Then from equations (V— 3 ) and (V— 4 ) 

_ f— 8.465 X lO^lbs/f t/sec x 5.064 ft/sec, 
^ 11 1.07 X 10 lbs 



— 42.866 X 10 T r ,“3 n rv^rv^o-. 

= — 40.0623 X 10 = — 0.040623 

1.07 X 10 



X 6 n 



5 r ec/sec 

„ _r'^n'^“mT _ rl .72 X 10 Ibs/rev/sec x 30 rpm x 60 rpm t 

^ 11 1 . 07 X 10 lbs 



1.72 X 10 X 0.5 
1.07 X 10 ^ 



= 0.08037 



Y V 4 

ir _ r m, f— 8.76 x 10 Ibs/ft/sec x 5.064 ft/sec, 

f'l c ~ In ^ ■* 

^12 1.864 X lO'" lbs 

4 

—44 36 X 10 —2 

= = - 23.798 x 10 = - 0.23798 

1.864 X 10 



( Y , —mu ) Jp Y, ip—mu ^ 

t, _ r ip. o'^m, _ r > 1 ' o ^m, 

^16 D D 

12 ^12 

,-2 



.6 2 



_ ^ 1 . 042 x 10 lbs/rad/secx 3 . 49 x 10 rad/sec— 1 . 166 x 10 ft/sec x 
1.864 X 10 ^ lbs 



—2 

x 25 . 32 ft/sec x 3 . 49 x 10 rad/sec 



] 



_ r 3 . 6365 x 10 ^- 29 . 5 x 10 ^x 3 . 49 x 10 _ ^ 3 . 6365 x 10 ^- 10 . 303 x 10 ^^ 

I g J - 1 g J 



1.864 X 10 
.5 



1.864 X 10 



_ 6 . 66 x 10 ^ - 3.5762 x 10 ^ = - 0.35762 

1 . 864 x 10 
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K 



17 



_ _ r~2 . 0 96x10 ^radTsec^xO . 8 7xl0~’rad/sec'^ T 

~ l~K J — I c. J 



«7— 2 



3 — — 2 



D 



12 



1.864 X 10 lbs 



1.82 X 10 ^ -0.9764x10 ^ = -0.009764 

1.864X 10 



K 



18 



_ ,4.823xl0^1bs/rad x 0.35 rad, 
l-^ J - I g J 

^12 1.864 X 10 



K 



19 



1.688 X 10~ 
1.864 X 10^ 



= 0.090558 



Y Sn 



lbs 



rev/sec 



= [■ 



n iTi j _ 1.936x10 rev/sec x 30rpm 60 rpm j 

^12 1.864 X 10^ lbs 



-1.936x10^x0.5 -0.968x10*^ 



= -0.005193 



1.864 X 10 



1.864x10 



K 



111 



= [ 



_ r~l • 307xlo”^f t/f t/sec x 5.064 ft/sec, 
I'. J ~ t n J 



13 



2.777 X 10 Ibs-ft 



K 
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6.6186x1^ = -2.3834 or -0.23834 into a gain of 10 



2.777x10 



M 1 Q lbs— ft _2 

r r~'^ • 462xl0^rad/sec x 3.49x10 rad/sec, 

D 7 ^ 

13 2.777 X 10 Ibs-ft 



K 



113 



K 



114 



-15.5723x10 

2.777x10^ 

N. •O’ 



= —5.6076 or —0.56076 into a gain of 10 



lbs— ft 
>7 t2 



_ ^ V iH j _ ^ —1.552x10 ft/sec^ x 0.844 ft/sec j 



D 



13 



-1.30988x10 

2.777x10^ 



= ^ -1» 188x10 



2.777 X 10 Ibs-ft 
= -0.47169 

lbs— ft 



8 



rad 



X 0.35 rad 



D 



13 



2.777 X 10 Ibs-ft 



-4.158x10 

2.777x10 



= —1.497 or —0.1497 into a gain of 10 



30 



Ibs-ft 



N 6n 



rev/sec 



K 



115 



r n 110x10 rev/sec x 30RPM 60RPM , 

I D J ” ^ 7 ^ 

^13 2.777 X 10 Ibs-ft 



5.11 X 10^x0. 5 2.555 x 10^ 



= 0.0920 



2.777 X 10 



2.777 X 10 



The calculated potentiometer values are summarized in 
Table VI . 

Making use of the k's coefficients, equations 
(V— 4) can be written as follows: 

~ ~ (V-4a) 



u = kj^j^u + 








i + 


ki J + 




ki8<SR + kj^g6n 


* “ ’'aiP* 






k^l4«R+kii56n 



Analog programming of these equations will give 
the left hand side quantities but sign inverted, because 
summer amplifiers and integrators invert the sign of 
their input quantities. Hence for analog programming 
equations (V— 4a) are written as follov\;s next; 

(V-4b) 



u = u 


+ 


kj2 6Til 








* = -lki5 V 


4 


^16 * * 


ki7 i 


+ 6R + 


6n] 




+ 


kii2» + 


*=113* 


+ + 


kii5(5n] 



Equations (V— 4b) can be written as follows after 
taking care of the actual sign of each coefficient as 
already have been calculated previously; 

~ (V-4c) 



5 = ki4u- 


ki 2 6n 








« = ki5 V + 


*=16 * ■* 


ki7 t - 


^18 ^ ^ ^ 


*=19 


^ ' ’=111'' •* 


’=112’’' ■*■ 


'=113* * 




*=115 
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Equations (V— 4c) are finally programmed in the 



analog Cl— 5000 computer. 



TABLE III 



DERIVATIVE 



DIMENSIONAL VALUE 





1.042 X 10 




-2.096 X 10 


Y 

V 


-8.76 X 10^ 


Y. 

V 


-1.083 X 10 


N 

V 


-1.307 X 10 


N. 

V 


-1.552 X 10 


Y 

n 


-1.936 X 10 




-1.762 X 10 


^6R 


4.823 X 10 


''6R 


-1.188 X 10 


^u 


-8.465 X 10 


X. 

u 


-1.0295X 10 


X 

n 


1.720 X 10 


N 

n 


5.110 X 10 


N, 

'I' 


-4.462 X 10 
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TABLE IV 



DERIVATIVE NON-DIMEN . FACTOR NON-DIMEN . VALUE 



(X. - 
u 


m) 


1 t3 
jpL 


-850 X 10“^ 


<^v- 


m) 


1 t3 
jpL 


-1500 X lO”^ 




m u^) 


1 t 3 
•^pL u 
2*^ o 


-510 X lO”^ 


%- 


L) 


1 t5 

jpL 


-77.9 X 10“^ 






1 ^3 
■oPL u 
2*^ o 


.0028 






2PL 


-. 0027 


Y 

V 




1 t2 

■^pL u 
2^ o 


-.01243 


y. 

V 




1 t3 
2PL 


-0.0071 


N. 

V 




1 t4 

jpL 


-.0002 


N 

V 




1 t 3 
■opL u 


-.0035 


Y 

n 




1 t 3 
■^pL u 
2^ o 


-.0000052 






1 t3 
2PL 


-.00043 


^6R 




1 ^2 2 
• "^PL u 
2*^ o 


- .0027 


^6R 




1y^t 3 2 
o-pL u 
2*^ o 


-.00126 


X 

u 




1 t2 

Ts-pL u 
2^ o 


-.0012 


X. 

u 




jpL 


-.0007 






1 t 3 
^pL 


.00005 


N 

n 




ipL^u 
2*^ o 


-.0000026 


N, 

'I' 




1 t4 

2P^ 


-.00227 


m 




1 t3 
2PL 


.0079 
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TABLE V 



PARAMETERS AND MAXIMUM VALUES 



m, = 1.26 X 10 (slugs) 



— mass of ship with entrained 
water 



m 

I 

2 



m 



m^/1.08 = 1.166 (slugs) — mass of ship itself only 
1.43 X 10^ Inertia (slugs — FT2) 

0.26 Yaw Angle (RAD) = 14.9 (degrees) 



ip = 0.0349 Yaw Vel . (RAD/SEC) = 2 DEG/SEC 



'm 



'Prr, = 



' m 



X = 



om 



0.00087 Yaw Accel. (RAD/SEC^)= 0.05 DEG/SEC^ 

5.064 = y (FT/SEC) 

■'em 



om 



^om 

6n = 
m 

y = 

■'om 
L = 



= 0.844 = y 



om 



(FT/SEC ) 

550 (FT) 

20 Rudder -Defelction (DEG) = 0.35 (RAD) 

30 Propeller Speed (RPM) = (RAD/SEC) = ^(REV/SEC) 

160 (FT) 

527.85 Ship Length (FT) 
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TABLE VI 



COEFFICIENT VALUES 



*^11 


= 


0.0404623 










*^12 


= 


0.08037 










*^15 


= 


0.23798 










^16 


= 


0.35762 










ki7 


= 


0.009764 










*^18 


= 


0.090558 










•^19 


= 


0.005193 










^111 


= 


0.23834 (into 


a 


gain 


of 


10) 


*^112 


= 


0.56076 (into 


a 


gain 


of 


10) 


^113 


= 


0.47169 










^114 


= 


0.1497 ' (into 


a 


gain 


of 


10) 


’^IIS 


= 


0.0920 











For ship No 2 the following relationships hold: 

^21 " ^11 

^22 ~ ^12 etc. 
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TABLE VII 



MARINER CHARACTERISTICS 



Length, ft. 


528.5 


Beam, ft. 


76.0 


Draft, ft. 


29.75 


Displacement, tons 


20.900 


Block Coefficient Cj^ 


.6125 


Prismatic Coeff. C 

P 


.6246 


Midships Section Coeff. C^ 


.9807 



4 . Analog Patching Configuration 
a. Static Test 

(1) Ship No 1 . Figure 5 shows the analog 
patching used for the static test. The 9300 digital ma- 
chine was used for this job. Computer program lA con- 
tains the source deck of the program as well as the 
print out. Table VIII gives the values of parameters as 
used. 

(2) Ship No 2 . Figure 6 shows the analog 
patching used for the static test. Computer program IB 
gives the source deck and Table IX gives the values of 
parameters as used. 

The static test was satisfactory since 
the printed values were the correct expected ones. As 
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After the coordinates transformation 




TABLE VIII 



Potentiometer Adjusted Parameter Assigned I.C. 

Address Value of integrators 

for static test 



POOO 


0.1000 


An 


AOOl 


3. 


volts 


pool 


0.0300 


IC-AOOl 


A003 


6. 


volts 


P002 


0.0803 


''12 


. A005 


9. 


volts 


P003 


0.0600 


IC-A003 


A013 


12. 


volts 


P004 


0.0404 


"ll 


A017 


15. 


volts 


POOS 


0.0052 


Ki9 


A051 


3. 


volts 


P007 


0.0900 


IC-A005 


A057 


45.4 


volts 


POlO 


0.2379 


"l5 








poll 


0.2383 


"ill 








P012 


0.4717 


"ll3 








P013 


0.0250 


[ip />P ] 

^m' ^m 








P014 


0.1200 


IC-A013 








P015 


0.5607 


"ll2 








P016 


0.1340 


^m ^m 








P017 


0.1500 


IC-A017 








P020 


0.0920 


"ll5 








P021 


0.0097 


"l7 








P022 


0.3576 


"l6 








P053 


0.0905 


"l8 








P054 


0.0500 










P055 


0.1497 


"ll4 








P050 


O’. 0337 










P051 


0.3000 


IC-A051 








P056 


0.0092 


[* An,] 
m m-* 








P057 


0.4540 


IC-A057 








P437 


0.1000 


AR 








P006 


0.1660 


[u /u ] 
m' m 
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TABLE IX 



Potentiometer 

Address 


Adjusted 

Values 


Parameter 


Assigned I.C. to 
integrator for 
static test 


P024 


0.1000 


An 


A035 


3. volts 


P025 


0.0300 


IC-A035 


. A027 


6. volts 


P026 


0.0803 


^22 


A033 


9. volts 


P027 


0.0600 


IC-A027 


A041 


12. volts 


P030 


0.0404 


r- 

(N 


A043 


15. volts 


P045 


0.0920 


^215 


A065 


30. volts 


P031 


0.0052 


''29 


A061 


54.5 volts 


P023 


0.0905 


■<28 






P034 


0.1660 








P033 


0.0900 


IC-A033 






P046 


0.4717 


"213 






P032 


0.2379 


"25 






P035 


0.2383 


"211 






P047 


0.1497 


"214 






P036 


0.0097 


K27 






P037 


0.0250 


^m ^m 






P040 


0.1200 


IC-A041 






P041 


0.1340 


['1' ] 
^m ^m 






P042 


0..1500 


IC-A043 






P043 


0.5607 


"212 






P044 


0.3576 


"26 






P407 


0.0337 








P425 


0.3000 


IC-A065 






P426 


0.5000 








P052 


0.5000 








P417 


0.0092 








P427 


0.5450 


IC-A061 






P416 


0.1000 


AR 
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far as the logic board is concerned Figure 6c shows the 
wiring necessary for the amplifier's operation where RT 
means real time and is connected to a capacitor of IpF. 
This Figure was made for amplifier AOOl and A003, but 
this is true for all the amplifiers used. Also, Figure 
6d shows the logic operation for the amplifiers AOOl and 
A003, which is again true for all the amplifiers used. 




B. CHARACTERISTIC LINEAR RESPONSE OF MARINER 

First of all, amplifier AOll was used as a timer for 
t h 

the 8 channel of the chart recorder. The wiring was as 
shown above . 

So by this configuration, a pulse was recorded every 
1 sec, i.e., every one time unit. 

Figures 7 and 8 show the analog patching diagrams used 
to obtain the responses of ship 1 and 2, respectively. 
These response curves are shown in a later section. At 
this time interaction forces and moments are not included. 
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Figure 6c. Logic patching diagram for the amplifiers used. 



R1 

o- 






/ 




— O A001-R 
— O f\ooi-R 




P\OOi-H 
AOO'l- H 



RT O 



TS O 



O AP 01 - i.O 

-Q f\ooi - \.0 
"O AOO'l - o.oA 

_Q f\0O*>- O-O'l 



0.1 Kc 




RT=;TS -0 



IW Rt.P\U TirAt. 

rw Timt lCAi-£ 



tion regarding the feedback capacitor, 
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f 

I 

I 

I 




Figure 7. Analog diagram for study of 
of ship No. 1. 










Figure 8. Analog diagram for study of dynamic behavior 
of ship No. 2. 
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m 









The dynamic behavior of the ships separately is of inter- 
est. Also note that all the initial conditions were set 
equal to zero. 

1 . Dynamic Test 




_ 1 _ 1 

RC (lMf2) (lyF) 



Amplifier AOOl— A026 



®o ^RC ®in RC 
s s 

RC~ ^in “ RC^ ®o 
s s 

_ 1 _ 1/RC 

e . 1+RC S+l/RC 
in s 



1 sec 



Settling time = = 



4 sec 



Amplifier A003-A027 




®o RC ^in RC ^o 

s s 

% = 1 ^ 1/RC 

e . RC +a S+a/RC 
in s 

_ 1 _ 1 _ 1 

a/RC 0.0404 0.0404 

(lMf2) (lyF) 



Settling time = T^ = 100 sec 



= 24.8 
sec 



Amplifier A017— A043 




1 



1 _ 1 
RC (iMfi) (lyF) 



1 sec 



T = 4 sec 
s 
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Amplifier A013— A041 




For the force free case it is: 



^2 ° 



Hence : 



and 



e = 103 e 
o o 






. 1 

®o ^^®o *^^o 

so X + :^^(103ct)x = 0 



Thus the characteristic equation becomes 
1 



S + 



Rj^C 



(103a) = 0 



Hence : 



T = 



R^C 



1 _ioo 

10a3 (10) (0.5607) (0.0250) 14 



= 7.15 
sec 
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Thus : 



T = 4 X 7.15 = 28.6 = 30 sec 
s 



By another way the same results can be obtained, i.6 



-e = 
o 



a 



R,C 



and 



e = 
o 



e = 
o 



^ r •• 



1 

- lOBe 



substituting this back gives: 



. _ a 

R,C 



log / dt 



and 



e = 
o 



lOag 

R^C 



/ e dt 
o 



or 



_ IQgg E (s) 
(s) Rj^C s 



Hence 



= 1 + iQgg 



sRj^C+lOag 

sRj^C 



sR^C 

TMg 



+ 1 



SR^C 

lOag 



and 



T = 



Rj^C 

lO^ "" (10) (0.5607) (0.0250) "" OTTT 



1x1 



Tg = 4 X 7.15 = 30 sec 



7.15 
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Amplifier A005— A033 

Again for the force free case only the closed loop 
is examined. 

- — / wdt 



= -[f + k2l0x] 



= -[k^lOx] 

or X = — k^w = — k^ [k 2 l 0 x] 

X + k^k^lOx = 0 

The characteristic equation is: S+10k^k2 = 0 

X =z - = 1 = i =25 

■^12 10(1.66x10 -^) (2.379x10 3.94x10 sec 




or for f=0: w 



T =4x2. 5= 10 sec 
s 

Figure 9 shows next the response when approximately 
+ 5° of rudder angle is applied. Note' that + 1 computer 
unit corresponds to + 20°. So a pot setting of P 437 for 
ship 1 and P416 for ship 2 equal to 1.0 corresponds to 20°. 
The sign of AR will be determined of course from the sign of 
the reference voltage. In fact for every curve taken by 

the chart the following relationship holds: 

+ 100 _ A [divisions] x s [volts/divisions] 

^m^^^ units] ~ x[in units] 
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Mariner 




Figure 9b. Dynamic test of Mariner. AR = + 5° 
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Figure 10a. Dynamic test of Mariner. 




Figure 10b. Dynamic test of Mariner. AR = + 5° 
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Figure 11a. Dynamic test of Mariner. 




Figure 11b. Dynamic test of Mariner. AR 



+ 10 ° 
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where: x = maximum expected value in actual units used 

m 

in magnitude scaling the equation 
A = divisions measured 
s = sensitivity in volts/division 
X = actual value of parameter x in actual units 
100 volts = 1 computer unit 

In Figure 10 the time scaling is changed and the 
linear response of Mariner is obtained in terms of u , tj; , i]^ , lO 




Figure 12. Steering Control used for the dynamic test 
of Mariner. 



when AR = 15° is applied. In Figure 11 the response for 
AR = + 10° is shown. The time constants estimated pre- 
viously are the same as the obtained ones. 

2 . Steering Control 

a. Analysis of Steering Control 

Up till now a simplified model for the repre- 
sentation of the steering control was used, which produces 
step type commands + or — according to what reference 
voltage is chosen. The steering control configuration is 
shown in Figure 12. 

One question which arises immediately after 
looking at Figure 12 ; is this a physically realizable 
steering control? The answer is no for several reasons. 
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Primarily the actual steering has a certain time lag in 
its response. Also for a given rudder change angle by the 
helmsman the actual rudder accelerates at the beginning 




and then decelerates at the end of its final position. So 
to be closer to the physical situation of a steering con- 
trol a definite time lag for the rudder action is needed 
as well as certain limits on the rate of change of the 
rudder angle, but without changing the initial gain set. 
The time lag can be produced by the above shown circuitry. 
The feed back gain is: k = 0.175 x 0.286 = 0.05 

The time constant and settling time are next calculated 

a/10 k/10 

®o RC ®in RC ^o 
s s 

^o _ d/lORC 
®in s+k/lORC 

• k = 0.5 

T = 1/tk/lORC] = 2 sec 
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ri! 




following can be deduced. The desired time lag is obtained 
as well as the big initial slope of AR. The desired res- 
ponse is the dotted line curve AR^ seen in the previous 
sketch. To achieve that a limiter is used which limits 
the output e^ as it can be seen from Figure 13a. 




The function of the limiter is similar to the 
saturation effect as next shown. Initially the error sig- 
nal [6e/6e^] hits the + 7° line v/hich corresponds to +100 
or —100 reference voltage, since the feedback is zero. 
Right after that due to the feedback the error signal is 
reduced almost linearly to zero. The time required for 
this is the desired time lag. Figure 13b illustrates the 
various steps made for the improvement of steering control 
in terms of the output rudder command. Figure 14 shows 
the complete analog configuration of the steering control 
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2. Command with time lag 




3. Command with both time lag and rate of change 
limitation 



Figure 13b. Improvement steps of rudder commands. 
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Figure 14. Steering control for ship No 



b. Time Lag and Settling Time of Improved 
Steering Control 



From Figure 13a it can be seen that: 



e . = — /a e, dt 
out 1 



e, = e, 

1 le 



im 



+ — 2. e + — ^ e . 
10 out 10 in 



where 



(e, I = lOO, -100 < <+100 

^im 



j . = — /a [e^ + e + e. ] dt 

out le. 10 out 10 in 

im 



For the force free case and within the specified limits 
it is : 



= — /a e .. dt 



'out 



10 out 



e . = —a tt— e . 
out 10 out 



or 



e . + a e . = 0 

out 10 out 



The characteristic equation is: 



from 



and 



T = 1| [ap/10] 



Now for: 6R = 20° 
m 



6e = 
m 



6R 



P = Y = 



a = 



m 



106e 



m 



k6e 



m 



6R 



7° 

k = 0.5 
3 = 1 

= 0.286 
= 0.1750 



for real time simulation and 



m 
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it is obtained: T = 2 sec 

and T =8 sec as before, 
s 

Figure 15a and 15b shov; the responses for input commands 
AR = + 5° and An = + 9 RPM but with different time scaling 
respectively. Figure 16a and 16b show the responses for 
input commands AR = + 10° and An = + 9 RPM but again with 
different time scaling of the chart recorder. It is ob- 
served that the expected time constants and settling times 
can be measured. 

C. TRANSFORMATION OF COORDINATES 
1 . Introductory Discussion 

The responses observed in the previous sections 
v;ere taken with respect to the Eulerian system., i.e. moving 
coordinate axes on the ship. This makes difficult the 
measurement of the separation distance between the two 
ships during the UNFIEP operation. It is desired therefore 
to use the space coordinate system, which for ease of com- 
puter print out is made to move ahead at 15 knots. It is 
also noted that the ship initially moves straight ahead at 
15 knots. Figure 17 illustrates the above discussion. 

From the inspection of Figure 17, equations (V— 2) 
can be recalled, i.e. 

x^ = u cos — u sin rp 

y^ = u sin + u cos \p 
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Figure 15a— 1. Dynamic test of Mariner. AR = + 5° , An = + 9 RPM. 











Hi 



m 



w 



iii 










Figure 15a— 2. Dynamic test of Mariner. AR = + 5°, 
An = 9 RPM. 
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Dynamic test of Mariner. AR=+5°, An=+9 RPM. 
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ih 
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Figure 15b— 2. 



Dynamic test of Mariner. AR 
An = + 9 RPM. 



+ 5 % 



65 




66 



Figure 16a— 1. Dynamic test of Mariner. AR = + 10°, An = + 9 RPM. 



im 














Figure 16a— 2. Dynamic test of Mariner. AR = + 10°, 
An = + 9 RPM. ~ 
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Figure 16b— 1. Dynamic test of Mariner. AR = + 10°, An - + 9 RPM. 




■'35 







Figure 16b-2. Dynamic test of Mariner. AR = + io®. 
An = + 9 RPM. “ 
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At this point computer unit XDS— 9300 is used in hy- 
brid operation, needed for the coordinates transformation, 
with the analog unit Cl— 5000. Computer program IVA and IVB 
show the source decks used and Figure 18a shows the control 
statements used for both digital and analog computers. 

Figure 18b shows the communications trunk lines used between 
digital and analog computers. Figures 19a and 19b show the 
analog patching diagrams used for the coordinates transfor- 
mation for both ships respectively. Table X contains the 
potentiometer values used. It should be mentioned here that 
for the control of the hybrid operation for ship No 1 the 



logic 


switches 


"0" 


and 


"1" were used. 


For 


ship No 2 the 


logic 


switches 


"2" 


and 


"3" were used. 


The 


hybrid operation 






U-o- 


^ 


1 1 1 1 - 








/ 

r 



V 

Figure 17. Space coordinates system 

was performed with both the switches being in the middle 
position . 
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Figure 18a. Flow chart for control of analog and digital 
unit . 
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Figure 18b. Communication trunk lines between digital and 
analog computer used for the coordinates 
transformation . 
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(After coordinates transformation section) 

Figure 19a. Analog diagram used for hybrid simulation 
for ship No. 1. 
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(After coordinates transf ormation section) 

Figure 19b. Analog diagram used for hybrid simulation 
for ship No 2 . 
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TABLE X (continued) 
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2 . Computed Linear Response of Mariner after 

Coordinates Transformation 

Figure 20 shows the linear response of the Mariner 
to approximately + 10° changes referred to midships in the 
rudder angle. 

Initially the Mariner is moving straight ahead at 
15 knots with the rudder at midships (6=0), and then 
AR = + 10° is applied for 25 seconds. Three phases are 
observed for a turn. 

a. Initial Approach Phase in which 
v=0, f=ijj=0, v=0, r=ij;=0, \p=0 

b. First Phase in which 
f>0, f=ip<0, v=0, r=4;=0 

This phase starts at the instant the rudder 
is turned, and ends before the rudder reaches its full 
deflection angle. The rudder force rudder mo- 

ment (Ng6) are dominant and produce accelerations which 
are opposed only by the inertia of the Mariner. 

It should be noted here that the rudder force 
Yg6 is towards starboard, since the rudder is at the 
stern of the Mariner, which corresponds to a port (nega- 
tive) turn although the transverse velocity v is positive 
(towards starboard) . 

c. Second Phase in which 

v>0, v>0, r=ifi<0, r=}p<0 

It must be mentioned here that a force Y^^v 
is created towards the center of the turn (towards port) 
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Figure 20—1. Characteristic linear response of Mariner at AR 




Figure 20—2. Characteristic linear response of Mariner 
at AR = + 10° 
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which is balanced eventually with the outwards centrifu- 
gal force of the Mariner. The acceleration v is zero 
when this balance is accomplished. 

The third phase of the turn is not developed 
because the rudder is given a negative —10° angle and 
the Mariner starts to turn starboard. It should also be 
noted here that at the instant the rudder is at —10° the 
center of the gravity of Mariner is at maximum port (nega- 
tive) position. The process of turning +10° and —10° 
angle to the rudder was performed several times. 

Figure 21 shows the linear response of the 
Mariner due to + 5° angle of the rudder. The difference 
with the response obtained with + 10° rudder is basically 
the smaller amplitudes. 

Figure 22 shows the linear response of the 
Mariner due to a positive change of R.P.M. of the pro- 
peller speed of approximately 9 R.P.M. with respect to 
the equilibrium R.P.M. (corresponding to 15 knots). 

The reference axes moving at a 

velocity 15 knots initially straight ahead as is the 
Mariner and consequently the (x,y) axes. The response 
is recorded by curves of calculated changes in the para- 
meters i|),ilj,i(i,y^,y^,x^,x^. Initially the origin of the 
(x,y) axes in the Mariner with respect to the(x^,y^) axes 

is 250 ft on the x axes and 48 ft on the y axes. Note 

o o 

here that the computer unit, i.e. 100 volts, equals 160 ft 
for y^ axes and 550 ft for x^ axes respectively. 
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Figure 21—1. Characteristic linear response of Mariner at AR 





Figure 21—2. Characteristic linear response of Mariner 
AR = + 5° 
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Figure 22-1. Characteristic linear response of Mariner at An = + 9 RPM 



-- 








































^ 1 




















T 

1 

1 




1 J 


111; 








I 






1 ■ 
























\ 




1 


1 




;>■ 








1 1 \ 




\ ' ' 


yi 






































";i ! 


’ 




1 ! 




1 1 


1',' 




1 


































iSii ' 






1 ' I ; 
1 ! 


'll: ' 


i yi 


1 






■ 


L 




























! 1 1 1 


1 ■ 1 


il , 






i'" i 


.iL. 


T‘ 








1 


























i 




ill' ' 


* ’ ; • 






" ' i 

■ 1 ■ 1 ' 


!l;.: 


ij" 


I'll 'i 




■ i 


I 


























:;i 


III 


'll- 

ll'l 


ll i 




111 




' Ii’' 








'■ 






























1 


!i' 


1' 


■' i 


Ip 


!'i' ' 




I y 








y 


















. 








■ ) 


"it"H 


1 


1 ; 




iiii 


il'M 


r- ‘ 1- ' 1 ' 






































' 1 
'li! 


jiii 


1 




ii,i 


ill; 


ji ‘ 


1 1_ . 








, 1 1 ' 




























’I;' 


'll' 1 

-l-Lf : 


1, 

1. 






y : 


'' 






i 








' y' 


- — P 






























V' 


yi ! 












1 

i 


. , 


























I'll 

i|il 


!; ! 


! ’ll 




y ' 


li;' ' 
















■ 


















1 




-- 


ill 


'|i! 


I'li 1 


lll' 








■ ' 












-i - 






1 . 




o: 

} 

i 

> 


s' 

1 

i 




















1 

1 

T 

T 

! 

1 

1 

It 


:!i 


ilil 


in' il 


h ill: 




1 1; 


ii:; ; 


;ii 111 

y ■ 


' 1 1 




' 'll’ 


1 




' ' ' 1 1 1 


i‘ 'i 




~ 


;;i 














r 


:r 


I' ■ 
















I; 






1 












• , 












r 












































■ 1 




y 








I’i 


i 


















1 


i * m 

\(/} 








< 








I < 


1 — 
(J : 


i 

1 






' 1 ' 


' ' i 

1 




A 




'ji4 










j 


( ' ^ 














-U 


rJ/> £>\ 


— n 






1 ' ' ; ' ' 


. i 
























ij 


\ 


1 




! : ! 1 


'Im 




. ' 




















j 


1 














I 




1 

1 








i 


























— 












i 






'1'': 


■ 


















Uil 


7^ 
















-4 


















3 








■ 

T 








1 1 




























5-S 






o 








'rC. 


012^ 














r.5 






i 


o 




i 


’7.: 

r 




















1 




1 •! 


h 




1T1 


1 








T. 




i 



Figure 22—2 



Characteristic linear response of Mariner at 
An = + 9 RPM 



The movement in the direction is due to the term 

in equation (V— 1) . The origin of the (x,y) axes in the 
Mariner moves to the origin of the axes in 110 sec 

approximately. The ship also yaws —1.5° approximately. 

It should be noted here that only the force 
and moment due to the RPM propeller change are simulated. 
The force and moment due to the propeller equilibrium 
speed of 15 knots is neglected. 

Figure 23 shows the characteristic linear res- 
ponse of Mariner obtained in terms of the calculated para- 
meters An,x^,x^,y^,y^,iJ;,ip,\|^ and caused + 5 RPM change of 
propeller speed. Again the Mariner is assumed moving at 
15 knots at equilibrium position. The initial position of 
the origin of the (x,y) axes was the same. 

Figure 24 shows the linear response of the 
Mariner to a rudder change of approximately +0.562° (for 
P437=+0 . 0281) and a simultaneous +9RPM propeller speed 
change in terms of the calculated parameter perturbations 
x^, x^, y^, y^, ij> , ip , An, AR. Again initially the Mariner 
is moving straight ahead at 15 knots with rudder at mid- 
ships. Then both the rudder and propeller speed change at 
+0.562° and +9 RPM respectively. This results in yawing 
moments which cancel each other (i]j= 0) . The only obtained 
perturbation is a change of location of the origin of the 
(x,y) axes on the Mariner but only along the x^ axes. 

Thus it is concluded that a rudder angle of about +0.562° 
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Figure 23—1. Characteristic linear response of Mariner at An - + 5 RPM 




Figure 23-2. 



Characteristic linear response of Mariner at 
An = + 5 RPM 
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Figure 24—2. 



Characteristic linear response of Mariner at 
AR = + 0.562°, An = + 9 RPM 
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is enough to create a yawing moment necessary to cancel 
the yawing moment created by a change of +9 RPM of pro- 
peller speed. Figure 25 shows the response of the 
Mariner obtained in terms of the calculated perturbations 
of the parameters An, x^, y^, y^, i}j, \p , AR and 

created by a simultaneous change in rudder angle and 
propeller speed at + 5® and + 9 RPM respectively. These 
curves can be analyzed analogously with the preceding 
curves. It is mentioned here that the same response was 
obtained for both ship No 1 and ship No 2 in all the pre- 
vious cases. 
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Figure 25-1. Characteristic linear response of Mariner at AR 
An = + 9 RPM 




Figure 25—2. 



Characteristic 
AR = + 5° An 



linear response of Mariner at 
+ 9 RPM 
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VI . 



INTERACTION FORCES AND MOMENTS ARE INCLUDED 



A. INTRODUCTORY DISCUSSION 

During the UNREP operation the tracking (receiving) 
ship is usually keeping course to avoid collision and 
maintaining station relative to the leading (replenish- 
ing) ship. The leading ship must keep a steady course 
and keep oscillations about this course to a minimum. 

For this reason the remainder of this study is subdivided 
into the following sections, namely: 

1. Phase I: Interaction forces and moments are 

included but are applied only to the 
leading ship as it is overtaken by 
the tracking ship. 

2. Phase II: This phase considers the complete 

hybrid simulation of the UNREP 
operation including all the inter- 
action forces and moments. 

It is noted that in phase I the response of the lead- 
ing ship is of interest. Further work in extending phase 
I will be the response of the tracking ship as it over- 
takes the leading ship so interaction forces and moments 
are applied to the tracking ship. This will be better 
understood after reading of next section, since it is 
related with how the interaction data curves are given. 
Briefly it can be mentioned that these data correspond 
to that case in which the interaction force and moment 
act on ship A as ship B overtakes it. So as a first case 
ship A is chosen to be the leading ship. 
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B. INTERACTION DATA CURVES-INTERPOLATION PROGRAMMING 



Figures 26 and 27 show the nondimensional N moment and 
Y force respectively versus the longitudinal separation 
(a) between midships and the lateral separation (3) • The 
generation of these curves by the analog computer exceeds 
its capability. Therefore these curves are made piecewise 
linear and stored on a digital computer namely the XDS— 
9300 machine. It should be mentioned here that in this 
stored array the points which correspond to 600 ft longi- 
tudinal distance (a) are included for both N— moment and 
Y— force. Hov^ever these points are not shown on Figures 
26 and 27 respectively. Computer program III shov;s the 
program necessary for a two-dimensional interpolation in 
two separate arrays for N— moment and Y— force respectively. 
It is mentioned that the longitudinal separation distance 
(a) can give negative or positive values but the lateral 
distance between ships (3) is restricted to positive 
values only. Figure 28 shows the geometry as far as 
longitudinal and lateral separation distances are con- 
cerned. It should also be mentioned that in phase I 
interaction moments and forces are applied only to 
ship A (see Figures 26, 27, 28) which corresponds to 
the leading ship because the data curves for the 
interaction effects are such. To apply interaction 
moments and forces to ship B, which corresponds to the 
tracking ship, the longitudinal separation distance 
should be sign inverted in order for it to give the 
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Figure 26. Dimensionless N Moment vs. Longitudinal Separa- 
tion . 
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Figure 27. Dimensionless Y Force vs. Longitudinal 
Separation. 
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Figure 28. Geometry of longitudinal and lateral separa- 
tion distances. 



appropriate value for M and Y respectively. M and Y 
are vectors with opposite direction in the case of ship 
B as compared with those applied to ship A. Ships A 
and B are defined in Figures 26, 27, 28. 



C. MEASUREMENT OF LONGITUDINAL AND LATERAL 
SEPARATION DISTANCE 

1 . Longitudinal Distance 

To measure the longitudinal distance, a, between 

midships the follov.'ing equation must be used: 



or 



a = X — X 

°1 °2 



—a = X — X 
°2 



(VI-1) 



substituting in this equation: 



a = 


a 


a 

m 




X 


= 


X 


X 


°2 




°2 


^2m 


X 


— 


X 


X 


°1 




°1 


Ot 

Im 



where x ,x , a are the maximum expected values of 
o_ ' o, m 

2m Im _ _ 

these variables respectively, and x ,x ,a are the 

°1 °2 
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scaled variables v/hose values range from —1 to 
equations (VI— 1) become: 

—a a = X X — X X 



m 



o„ o 



2m 



o, o 



Im 



+1. Hence 



(VI-2) 



or 


—a = 


X X 

On Ot 

r 2m. — r Im. — 

[ ] X — [ ] X 

a o„ a o- 

m2 ml 


(VI-3) 


Choosing : 


a = 
m 


550 ft 






X = 

Oo 
2m 


550 ft 






X = 

Im 


550 ft 





equation (VI— 2) becomes: 

—a = X — X (VI— 4) 

02 

Figure 29 shov;s the analog patching diagram necessary for 
equation (VI— 4) 




Figure 29. Longitudinal distance analog diagram. 



2 . Lateral Distance 

To measure the lateral side— to— side distance be- 
tween ships the following equation must be used: 
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or 



B = y — y 
°2 °1 



■3 = y - y 

°1 °2 



(VI-5) 



substituting in this equation: 

3=33 
y 



o. 



m 

= y 



O, o, 

1 Im 



y = y y 

°2 °2 °2m 

where: y , y ,3 are the maximum expected values of 

O <1 ITlI 

Im 2m _ 

these variables respectively, and y , y ,3 are the 

• O 2 

scaled variables v;hose values range from —1 to +1. Hence 



equation (VI— 4) can be written as: 



— 3 3= y y ~y y 

m -^o, -"o, -^o„ -^o„ 

Im 1 2n 2 



(VI-6) 



or 



Choosing : 



o r 1^1 - 

-3 = [-^ — ] y 

% °1 



r 2 m, - 

m A 



3 = 100 


ft 


m 




= 160 


ft 


2m 




= 160 


ft 


°lm 





equation (VI— 7) can be written as: 
-3 = [1.60]y - [1.60]y 



o. 



(VI-7) 



(VI-8) 
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Figure 30 shows the analog patching diagram needed for 
equation (VI— 8) . 

It should be noted at this point that because 

the coefficients of y and y in equation (VI— 8) are 

°1 

greater than 1, the pots are set as follows: 

P040 = 0.1600 into a gain of 10 
P052 = 0.1600 into a gain of 10 
During the simulation 3 can take neg-ative values 
but the data curves are given for positive values of 3 
only. Thus in some v;ay an absolute value should be pro- 
vided either by digital subroutine ABS or by analog— patch- 
ing configuration. The latter is chosen since the use of 
digital machine adds time delay in the simulation which 
may give rise to an error roughly of 5%. 




Figure 31a shows the analog patching diagram neces- 
sary to obtain the absolute value of 3 • 

It should be noted that in Figure 31a in the second 
input of A024 summer the gain is 10 but since: P405 = 0.2000 
the gain becomes finally equal to 2. Also the limiter L 10 
limits the output of A066 to negative values. More speci- 
fically — |3| corresponds to +1 computer unit limit and 
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+ |3| corresponds to —1 computer unit limit respectively. 
However according to the following illustrated fact: 




there is a better way of getting a gain of 2 since two 
resistors of value 1.0 M in parallel gives a value of 0.5 
as value of input resistor which corresponds, to a gain of 
2 as shown next in Figure 31b. 




inpu-t 



Figure 31b. Absolute value patching diagram of 3. 



Finally the way of Figure 31b is chosen being more compact. 

D. APPLICATION OF INTERACTION MOMENT N AND FORCE Y 
TO SHIP A (LEADING SHIP) 

1 . Introductory Discussion 

Ship No 2 is chosen arbitrarily to be the ship A, 
that is the leading ship, on which moments and forces due 
to interaction effects would be applied. Figure 32 illus- 
trates the way N moment and Y force are interpreted by 
hybrid operation of the machines. 

The digital program which takes as data the values of 
3f a, in the computer mode of hybrid operation, finds the 
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Figure 32a. Interpretation of N— moment and Y— force for 
hybrid computer operation. 



corresponded values of N and Y and then gives back to the 
analog the values of N and Y, wliich are scaled quantities, 
takes care of: 

a. Rescaling of 3/ ct since 3 = 1.0 and a = 1.0 
corresponds to 100 ft and 550 ft respectively. 

b. Making linear interpolation in the two- 
dimensional array. So it finds the corres- 
ponded nondimensional quantities of N' and Y'. 

c. Since, from the data curves, Figures 26, 27, 



N' = 44 . 5 X 


-5 


10 


max 


—5 


Y ' = 90 . 0 X 


10 


max 





hence, finding the dimensional values of 

and Y ' namely 
max 



N 

max 



N' X 4pb^u^ = 44.5 X 10 ^x 9.428x 10^^ 
max 2 o 

4.195 X 10^ Ibs-ft 
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and 



. 12 2 

Y = Y' X ^pL u 

max max 2^ o 

= 90.0 X 10”^x 17.87 X 10^=1.608 x lO^lbs 
= 72.9 tons 

Finding from N' and Y' the dimensional values 

of N' and Y' respectively as follov/s: 

N = N'x = N' X 9.428 x 10^^ 

2 o 

Y = Y'x ipL^u^ = Y' X 17.87 x 10^ 

2^ o 

Scaling both N and Y quantities taking in 
consideration the fact that the subroutine 
DAC multiplies each number given to the analog 
terminal by a factor of 100 i.e. it converts 
in computer units, so: 



NSCALE = N/N 

max 

YSCALE = Y/Y 

max 

Computer program VA contains the source deck as 
well as the potentiometer values and amplifier addresses 
used for phase I. 

Computer program VB contains the source deck used 
for phase I, but ship A (leading ship) was in that case 
the second model of the Mariner. The difference was only 
in the trunk lines which were used in the previouse case. 

In Figure 32 the trunk lines t 426, T427 belong to the ship 
model No 1. 

2 . Dynamical Representation of Interaction Moment 

and Force— Analog Programming 

Recalling Newton's laws of motion, equation (III— 1) 
and since the study is done for the ship to be in 3 degrees 
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of freedom, equations (III— 6) , (V— 1) , (V— 4) , the applica- 

tion of the external interaction N— moment and Y— force are 



going to be an additional term in the yaw and sway equa- 
tion respectively. Considering the dimensional form of 
the linearized equations of motion, equations (V— 4) will 
become for ship No 2 



Xu X 6n 

“T" •- 1,1 rn. -1 •- n m -i *7^ ^ 

u = — ] u + [— p:; ] 5u , Surgo 



D 



21 



D 



21 



Y V Y,-mu - Y.,ip _ 

,, r V m, — r w o m. , , r ifj m. 

U = [— — ] V [ — ^ + [_x — ] if, 



D 



21 



D 



22 

Y, 



Y;,„5R Y;, 6n 

, r OK m, . r on IUt , 

+ [— ] 6R + [ — 7 ^ ] 6n -f - 



22 

Y 



D 



22 



D 



22 



D 



22 



I 



Sv;ay 



additional term 



N V N.il _ N.v 

• 7 - r V m, — , r ni. , , r V m. -r 

Jf = [-- ] V + [-^ ] ^ [-^r ] V 

23 23 *^23 



N.^SR N 6n 

+ [-4^ — -] 6R + — -] 6n + 



N 



D 



23 



D 



23 



1^23 



Yaw 



additional term 



(VI-9) 



where: = (m— X.)u = 1.07 x 10^ lbs 

21 u m 

= (m-Y.)v = 1.864 x lO^lbs 

2 2 V m 

= (I -N.Jip = 2.777 x lo"^lbs-ft 
Y = interaction force 
N = interaction moment 

It is clear from equation (VI— 9) that the surge 
equation remains unchanged to interaction effects and 
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only sway and yaw equations have each of them one addi- 
tional term respectively. Thus for the analog program- 
ming of equation (VI— 9) it is only necessary to add the 
circuitry shov/n in Figure 33a to the patching diagram 
of Figure 19b. Note in this Figure: 



P042 



N 

max 



D23XIO 



4.195x10^ 

2.777x10® 



0.1511 



Here there is a factor of division by 10 since the value 
(N yD^^) > 1 and so it is necessary to feed amplifier 
A036 in a gain of 10. Note that since in amplifier A036 
we need one more, not available, input gain of 10, a 
separate internal resistor of 0.1 is used through the 
summing junction (SJ) of A036 giving the required gain of 
10. Also: 

P017 = = 1- *^08x10 ^ 0.1581 

^^22 1.017x10 



cJ -la) 





DigCfcodL 

- C^) 






0W.U 

T H-.23 - («•) 



Figure 33a. Analog Diagram programming of interaction 
N— moment and Y— force for ship A (No 2) . 
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For the ship model No 1 the analog patching cir- 
cuitry of Figure 33b was added to this of Figure 19a 
whenever it was desired to be used as leading ship and 
ship model No 1 to be used as tracking ship respectively. 




Figure 33b. Analog diagram programming of interaction 
N— moment and Y— force for ship model No 1. 

It must be mentioned here that for Phase II where 

interaction effects have to be applied on both ships, both 

the configurations of Figures 33a and 33b must be used. 

3 . Obtained Responses for Phase I 

a. Stationary runs (i.e. both ships have same 
propeller speed) 

Figure 34a, b shows the obtained linear res- 
ponse of tlie leading ship (ship A) in terms of the calcu- 
lated parameters Y, N, y^, a, 3, v, v, i/i . 

Originally the leading ship was placed at the 
origin of the space coordinate system 

tracking ship was placed at = +524 ft and y^ = +70 ft. 
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Figure 34a. 



( 1 : 







Figure 34b. 



Linear response of leading ship for Phase I. 
(AR=0, 6n=0, a=+524 ft, 3= 70 ft) 
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Then both the ships were left to run with the same propel- 
ler speed which was 15 knots. The same speed of 15 knots 
was assumed for the origin of the space coordinate system 

The leading ship was forced to yaw to positive 
yaw angle due to the additional term of the interaction 
moment N. The interaction force Y produced a lateral ac- 
celeration and lateral velocity which forced- the leading 
ship to move laterally decreasing continuously the lateral 
separation distance. In the sway equation (VI— 9) not only 
the Y— force term Wcis acting but the derivatives Y^,Y^,Yj^|.^ 
were producing additives effects since due to the yaw 
angle a change of yaw angle and rate of yaw angle change 
were obtained. It is also noted that the interaction force 
Y and moment N are slightly increased during the run since 
the lateral separation was decreasing. Similar results 
were obtained for the second run in which the leading ship 
was again placed originally at the origin of (^q/Yq) axes 
and the tracking at = —524 ft and y^ = +70 ft. Again 
the origin of the(x^,y^) was moving at 15 knots as well as 
both the leading and tracking ship. The response was 
calculated in terms of the perturbed parameters, i|^,y,N,a, 
3,v,v as shovm in Figure 35. The lateral distance during 
this second run was decreasing faster since the values 
for interaction force— Y and moment— N were greater. Ob- 
viously from the data curves both force Y and moment N 
are not symmetrical about the origin. 
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Figure 35. Linear response of leading ship for Phase 
(AR=0, 6n=0, a= -524ft, 6= 70 ft) 



I 



Figure 36 shows the response of the leading 
ship when the longitudinal distance and the lateral dis- 
tance were set equal to zero and 70 feet respectively at 
the beginning of the run. Now because the longitudinal 
distance is zero, the Y— force and N— moment have not only 
greater values than previously but tlie leading ship yaws 
to negative yaw angles. This has as effect a negative 
lateral velocity and an acceleration initially which 
eventually reduces the interaction force— Y and moment— N 
drastically, but the lateral distance is increasing nov/ 
since the velocity along the y^— axis has opposite sign. 

For all the above three responses the rudder 
for both the ships was set at zero, i.e. = 0. It should 
be noted at this point that the hybrid simulation program 
(computer program V) was made to work between lateral 
separation distances 100 ft and 50 ft since data are given 
for this range. The increase of the lateral distance is 
due : 

a. To the negative yaw angle. From equation {V— 2) 
it seems that eventually the term (u sin ip) 
becomes predominant over the term (u cos iji) 
hence even if a positive value is obtained for 
the V velocity a negative value of y is 
actually applied to the leading ship. 

b. The negative yaw angle is obtained because for 
a=0, 3=70 ft a negative moment N is applied. 

It has been mentioned previously that the two 

models are identical. So for Phase I No 2 model has chosen 

to be ship A, i.e. leading ship. To complete our work for 
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Figure 36. Linear response of the leading ship for Phase 
(AR=0, 6n=0, a=0, 6=+70 ft) 




I 




Phase I and check the wiring for model No 1 the case was 
reversed and ship A was the No 1 model and ship B the No 2 
model. Remember ship B is always the tracking ship. For 
this last case a change in wiring on the analog board was 
made for the circuit representing the absolute value of 
the lateral separation distance (3)« The change was due 
to the lack of operation of the limiter L 10 as was shown 
in Figure 31a and 31b. 

The following configuration is used from now 
on for the absolute value of the lateral distance (3) es 
is next shown in Figure 37. 




Figure 37. Absolute value patching diagram of the lateral 
distance . 

Just briefly is mentioned that it is desired 

to keep the voltage of point A between the limits 0 and — <». 

This is because diodes D , D are connected in such a way. 

12 

Now when point C has a positive voltage is OFF and 
point B has a negative voltage hence point A has a positive 
voltage since D 2 is conducting. V7hen point C is negative, 
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voltage is ON and B has about 0 or +0.5 volts voltage 
which is not enough to break down hence point A is at 
zero voltage. 

Figures 38 and 39 show the obtained response 

for the leading ship, when the longitudinal distance is 

+400 feet and —400 feet and the lateral distance +70 feet 

from the tracking ship. The response was obtained in 

terms of the perturbated parameters ijj , y, N, , a, 

V, V . The ships were moving at the same speed of 15 knots 

with the leading ship placed at the origin of the 

axis. In Figure 38 the tracking ship position was at = 

+400 ft and y = +70. The interaction force — Y was almost 
o 

zero during the time of this run, while the interaction 
moment ~N had a positive value producing a positive yaw 
angle. Because y^ was positive the track of the leading 
ship reduced the lateral distance (B)- 

Similar discussion can be made for the response 
of Figure 39 where the tracking ship was placed at x^ = 

—400 feet and y^ - 70 feet. The lateral distance was again 
reduced continuously. 

Several "stationary" runs, i.e. at the same 
speed of 15 knots for both the ships, were made for the 
longitudinal distance (a) being 160 ft, 300 ft, 200 ft and 
the lateral distance (B) 70 ft, respectively. Figures 40, 
41, 42, 43, 44, 45 show these responses in terms of the 
perturbated parameters of the leading ship, , Y , N , y^ , a , B ^ 
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Figure 40. Linear response of the leading ship for Phase 
(AR=0, 6n=0, a= -160ft, 3= +70ft) 
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Figure 41—1. Linear response of the leading ship for Phase 
(AR=0, 5n=0, a=+160ft, 3=-+70 ft) 
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Figure 41—2. Linear response of the leading ship for 

Phase I. (AR=0, 6n=0, a=+160ft, 3=+70ft) 
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Figure 42. Linear response of the leading ship for Phase 
(AR=0, 6n=0, a=-300ft, 3=+70ft) 
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Figure 43. Linear response of the leading shio for Phase 
(AR=0, 6n=0, a=+300ft, 6=70 ft) 
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Figure 44—2. 



Linear response of the leading ship for 
Phase I. (AR=0, 6n=0, a=+200ft, g=+70ft) 
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v,v. The analysis of these response curves is analogous 
to the previous ones. The lateral distance in some cases 
is reduced eind others increased by the interaction effects 
depending on the relative longitudinal position of the 
ships since the lateral distances were kept constant at 
70 ft at the beginning of each run. 

From the observation of these "stationary" runs 
it seems that, there are positions where both the interac- 
tion force and moment acting on the leading ship tend to 
draw the leading ship towards the tracking ship. Such posi- 
tions seem to be at the longitudinal distances of 524 ft, 

400 ft and 300 ft. 

An interesting effect was observed when the 
longitudinal distance is +200 ft and the lateral distance 
+70 ft. Figure 44. Although the interaction force— Y and 
moment— N had both positive values they did not affect the 
lateral distance (g) throughout the time of the run (2 
minutes roughly) because the interaction effects counter- 
acted each other, due to the fact that the positive inter- 
action moment was producing a yaw positive angle such as 
to counteract the positive interaction force. The the y^ 
velocity of the leading ship was zero in this case and so 
the lateral distance could be kept constant. At this point 
for the ease of the manual control the input voltage into 
the steering control block and into the propeller speed 
blcok is changed as Figure 46 shows for ship model No 1 and 
Figure 47 shows for ship model No 2. 
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Figure 46. Rudder and propeller input commands diagrams 
for ship model No 1. 




Figure 47. Rudder and propeller input commands diagrams 

for ship model No 2. 

Those analog configurations shown on Figures 
46 and 47 for both ship models respectively replace the up 
till now used configuration. This is done primarily because 
one can have negative or positive rudder and propeller 
speed changes without necessarily changing the switch posi- 
tion. This obviously gives operative flexibility for the 
operator making the control of the ship more realistic 
besides the necessary required sensitivity of the rudder 
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wheel and propeller. Then for zero rudder or propeller 
speed change either the corresponding switch has to be open 
or the corresponding potentiometer has to be set at 0.5000 
value. For example for the rudder angle will be 5 turns 
for 0 to 20° and 5 turns 0 to —20° v;hich corresponds to 4° 
per turn. Likewise it will be 5 turns for 0 to + 30 RPM 
hence it will correspond to 6 RPM per turn. For the ease 
of operation the potentiometer values for P4_37, P4 06, P416, 
P406 should be set in "RESET" mode having in mind that 100.0 
volts corresponds to 20° rudder and 30 RPM of propeller 
speed respectively. This avoids the confusion of converting 
in "POTSET" mode the potentiometer values to be set in 
actual volts since now the 0.5000 position of the potentio- 
meters corresponds to zero volts. 

b. Runs with different propeller speed between 

ship A (leading ship) and ship B (tracking ship) 

Tv^?o kinds of runs were made, one with the track- 
ing ship have a speed of +10 RPM greater than the speed of 
15 knots of the leading ship. In the second run the differ- 
ence in propeller speed was 5 RPM. Note that no controls 
were' applied . The following assumption was made: both ships 
can be placed at any desired initial position at the begin- 
ning of the run. This was done because no complete run with 
the tracking ship passing the leading ship was possible to 
be made primarily for two reasons : 

(1) Saturation effects of the analog computer 
amplifiers . 

(2) Interaction effects can be applied only 
■ for the lateral distance being between 

50 ft and 100 ft. 
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This assumption holds for Phase II also. A matter of fu- 
ture investigation would be a rescaling of the analog com- 
puter unit to avoid the saturation effects. Also more 
data should be available for the interaction effects to 
be applied for lateral distances outside, the range of 50 
to 100 ft. Nevertheless the tracking ship v;as placed at 
different initial longitudinal distances from —524 ft to 
+524 ft. All the responses of the leading ship for Phase 
I were calculated in terras of the perturbated parameters 
of the leading ship , Y , N , , a , 3 , v , v . From the observa- 

tions of the obtained responses, Figure 48 to Figure 64, 
it can be seen that during the approach and the departure 
of the tracking ship the leading ship tends to reduce the 
lateral distance, 8. So a possibility of collision appears. 
More specifically the tracking ship starts at —524 ft longi- 
tudinal distance with 10 RPM propeller speed greater than 
that of 15 knots and tries to overtake the leading ship. 
Figure 43. Due to interaction moment N the leading ship 
yaws to a positive angle and so tends to reduce the lateral 
distance, 3/ although the interaction force Y is negative. 

In Figure 49, the next run, the tracking ship starts to 
overtake the leading ship at a longitudinal distance of 
—400 ft with 10 RPM propeller speed greater than that of 
the leading ship. It can be seen that originally the yaw 
angle is positive due to the positive interaction moment 
and the lateral distance is decreasing again. This in fact 
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Figure 49 . Linear response of the leading ship for 
Phase I (AR=0, 6n=+10 RPM, initially 
a=-400ft, B=100ft) 
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is happening because the bow of the leading ship yaws to- 
wards the oncoming tracking ship. When the longitudinal 
distance is reduced to —200 ft the leading ship starts to 
yaw at a negative angle; this is because both the Y— force 
and N— moment changed sign, hence the lateral distance be- 
ing unchanged for about 15 sec at 50 ft approximately, 
starts to increase toward 100 ft when the longitudinal 
distance is approximately 160 ft. 

In the next shown run. Figure 50, similar analy- 
sis can be made except from the fact that the tracking ship 
now starts at initial position at x = —300 ft and y = 

+50 ft. The lateral distance being 50 ft for about 15 sec 
starts to increase when the longitudinal distance is approxi- 
mately —200 ft. 

Figure 51 shows the response for an initial 
position of the leading ship at (0,0) and the tracking ship 
at (0,50 ft). The tracking ship possesses a 15 RMP greater 
propeller speed than that of 15 knots. The leading ship 
immediately starts to yaw negatively and thus increases the 
lateral distance drastically. 

In Figure 52, where initially a = +160 ft and 
B = 70 ft it is seen that the lateral distance increases 
until the longitudinal distance is +200 ft. Then due to 
the interaction effects the leading ship reduces the 
lateral distance, since a positive y^ velocity has been 
created . 
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Figure 50. Linear response of the leading ship for Phase 
(AR=0, 6n=+10 RPM, initially a=-300ft, 3=50ft) 
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Figure 52. Linear response of the leading ship for Phase 
(AR=0, 6n=+10 RPM, initially a=+160ft, 6=70ft) 



Figures 53, 54, 55 show the response of the 
leading ship when a and 3 were initially set at +300 ft, 

+400 ft, +524 ft, and 70 ft respectively, at 5 RPM differ- 
ence in propeller speed. It seems that the best position 
for two ships of the same size is exactly abeam v;hen along- 
side. The increase of the lateral distance when the longi- 
tudinal distance is betv;een approximately —200 ft and +200 
ft implies that the leading ship yaws. Hence control has 
to be applied for it to keep a constant course. It is also 
obvious that the approach and departure of the tracking 
ship implies the use of control on the leading ship in order 
for it to avoid collision. 

Figures 56 till 64 show the linear response of 
the leading ship for different initial positions of the 
tracking sliip with a difference of 5 RPM in propeller speed. 
The tracking ship overtakes the leading ship, which has a 
15 knots speed. A similar analysis can be made for these 
responses as was done for the previously mentioned I'uns. 

4 . Obtained Responses for Phase II 

Computer Program VI shows the source deck used for 
Phase II. In Phase II interaction forces and moments are 
applied on both ships. The circuitry of both Figures 33a 
and 33b is used. It should be mentioned here that the 
digital machine computes for the leading ship the inter- 
action force and moment corresponding to those longitudinal 
distance, a, and lateral distance, 3/ which the analog com- 
puter passes through the trunk lines to the digital computer. 
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During the same cycle of computation the digital computer 
inverts the sign of longitudinal distance, a, and for the 
same lateral distance, 3/ gives another pair of Y— force 
and N— moment to the analog computer, v/hich are the inter- 
action effects applied to the tracking ship. 

a. Stationary Runs (i.e. both ships have same 
propeller speed) 

Figures 65 to 75 shov/ the linear response of 
both the leading ship and the tracking ship alongside at 
the same speed of 15 knots at various longitudinal and 
lateral distances respectively. The responses were calcu- 
lated in terms of the perturbated parameters ' 

a, 3/ y / 'I'd where the indices A and B mean the leading 
Or B 

ship (ship A) and tracking ship (ship B) respectively. 
Depending upon the relative position of the ships there 
are cases of attraction between them and repulsion between 
them. An interesting case was observed at the exactly 
abeam position, while the two ships were alongside at 70 
ft lateral separation distance. The lateral distance was 
kept constant at 70 ft throughout the time of recording 
although both ships were yawing at negative angles. 

b. Runs with Different Propeller Speed between 
Ship A (Leading Ship) and Ship B (Tracking Ship) 

The assumption made for Phase I that both ships 
can possess any desired initial positions with respect to 
the axes is carried throughout phase II, since 

again the tracking ship was never able to completely over- 
take the leading ship during the progress of a single run. 
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Figure 72. Linear response to interaction effects for Phase II. 
(AR=0, 6n=0, a=-300ft, 3=50ft) 
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The obtained responses were calculated in terms of the 
perturbated parameters Y^, N^, Yq / ot / 3 , 

where the indices A and B mean ship A and ship B respec- 
tively. With no controls applied on both ships the 
tracking ship was chosen to have in the first case 10 RPM 
propeller speed and in the second case 5 RPM propeller 
speed, greater than the speed of 15 knots, which was the 
leading ship's speed. 

In Figure 76 the initial position of the track- 
ing ship is at (—524 ft, 70 ft) and that of the leading 
ship at (0,0). Obviously the longitudinal distance de- 
creases since the tracking ship has a greater speed. The 

lateral distance, decreases also since the y velocity of 

°A 

the leading ship is positive and greater than tlie y^ veloc- 
ity of the tracking ship, which is also positive. This 
means that the leading ship is pulled towards the tracking 
ship. It is seen also that the leading ship yaws to posi- 
tive angles as v^7ell as the tracking ship but since the 
leading ship's yav-; angles are greater than those of the 
tracking ship the stern of the leading ship comes towards 
the bow of the tracking ship. 

Next in Figure 77 initially the longitudinal 
distance, a, and the lateral distance, 3/ were —400 ft and 
70 ft respectively. The lateral distance, 3/ was decreas- 
ing until the longitudinal distance became roughly —250 ft. 
At this point both the interaction force and moment on the 
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leading ship changed sign. This can be observed in Figure 
78 where initially a and 6 were —300 ft and 50 ft respec- 
tively. The above mentioned change of sign results in a 
negative leading ship's yaw angle, while the tracking ship 
continues to have positive yaw angle. The lateral distance 
starts to increase from the instant at which the longitudi- 
nal distance is roughly —250 ft. 

In Figure 79 a and B were initially —200 ft 
and 50 ft respectively. It can be seen that the leading 
ship yaws to negative yaw angles while the tracking ship 
maintains zero yaw angle. Although this is happening the 
y velocity is negative and the y almost zero. Hence 
the leading ship is pushed away from the tracking ship and 
the lateral distance is increased. 

The same analysis can be done for the response 
of Figure 80a where a and B v.'ere initially —160 ft and 70 
ft respectively as v;ell as for the response of- Figure 80b 
where initially was set a = —160 ft and B = 50 ft. 

Figure 81 shows the response obtained for the 
exactly abeam position of the two ships being alongside. 

It is seen although both ships are yawing excessively to 
negative angles the lateral distance is slightly increas- 
ing. Hence this is the best position for station keeping. 

Figure 82 sho^^7s the obtained response for a 
and B being initially +160 ft and 100 ft respectively. 

It is seen that the lateral distance B starts to decrease 
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Figure 78. Linear response to interaction effects for Phase II. 

(AR=0, 6n=+10 RPM, initially a=-300ft, 3=50ft) 
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when the longitudinal distance becomes roughly +200 ft, 
while the leading ship yaws slightly negatively and the 
tracking ship yaws excessively negatively. These facts 
bring the bow of the leading ship towards the tracking 
ship's stern. At +250 ft longitudina], separation between 
midships the leading ship reverses direction of yawing 
while the tracking ship continues to yaw negatively. 

This can be seen more easily in the response 
of Figure 83 where initially a and 3 set at +200 ft and 
100 ft respectively. Due to these opposite yaw angle 
directions the sterns of both ships are pulled toward each 
other wliile the lateral distance is decreasing. In Figure 
84 a and 3 initially were set to be equal +300 ft and 70 
ft respectively. The lateral distance is decreasing due 
to the fact that the leading ship's y^ velocity is posi- 
tive and greater than that of the tracking ship. 7ilso 
the leading ship yaws positively while the tracking ship 
yaws slightly negatively. This brings leading ship's 
stern towards the' tracking ship. Note that at roughly 
+350 ft longitudinal separation the tracking ship reverses 
the direction of yawing. This tends to bring the bow of 
the tracking ship towards the stern of the leading ship. 

Figures 85 and 86 show the responses for 
initially setting a and 3 +400 ft and +524 ft, 70 ft res- 
pectively. The lateral distance, 3? is increasing although 
both ships yav7 positively. 
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Figure 83. Linear response to interaction effects for Phase II. 

(AR=0, 6n=+10 RPM, initially a=+200ft, B=100ft) 
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Figures 87 to 97 shov/ the responses obtained 
for different initial positions of the ships. The analy- 
sis of these curves is analogous to the previous one, 
the difference being that the tracking ship has now a +5 
RPM greater propeller speed than that of 15 knots. Note 
that for the smaller difference in propeller speed there 
is a smaller rate of change of the lateral separation 
distance although both ships yaw to greater • values of yav/ 
angles since they have adequate time for that. 

c. Manual Control is Applied in Phase II 

From the previous obtained responses it is 
obvious that control must be applied in order to keep a 
desired course as well as a lateral separation distance, 

3. Up till now for lateral distance greater than 100 feet 
or less than 50 feet the digital programming \\?as made such 
that for these cases a zero force— Y and moment— N v^;as given 
to the analog computer. This fact made impossible for the 
tracking ship to completely overtake the leading ship dur- 
ing only one run.' At this point the following approxima- 
tion was made to override this difficulty, namely the 
digital programming. Computer Program VII, was made such 
that the digital compvater gives (for values of 3 outside 
the table limits) the values of Y— force and N— moment cor- 
responding to the extreme values of 3/ i.e. 50 ft and 100 
ft. Obviously in the actual case the interaction effects 
would be more pronounced for 3 less than 50 ft and less 
pronounced for 3 greater than 100 ft. But with this 
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approximation it was possilsle to get the response of Fig- 
ure 98 with control manually applied in terms of the per— 

turbated parameters ijj , AR, , N , y , a, 8, where 

A A A D B 

A and B indices mean ship A and ship B respectively. The 
tracking ship originally placed at = —200 ft, y^= 50 ft 
has a 10 RPM greater propeller speed than the leading ship. 

These responses do not give the actual needed 
rudder to maintain station and course unless the operator 
is practiced enough to know the behavior of both ships and 
hence is familiar with the amount of rudder actually needed. 
Nevertheless, these responses support the already known 
fact (from experience at sea) that underway replenishment 
is a controllable operation although experienced helmsmen 
are needed on both ships of the UNREP. 

E. CONCLUSIONS - REMARKS - SUMMARY 

Before any conclusion is derived it must be mentioned 
here that this study was done for the underway replenish- 
ment at sea of two equal sized merchant ships, namely 
MARINER class ships. The study also assumes calm v/ater 
conditions which consequently assumes motion of ships in 
3— degrees of freedom only. The hybrid simulation was done 
primarily for the no control plant although some remarks 
can be deduced for the nature of rudder (control surface) 
needed. Most of the observed results of this simulation 
study are already known from experience of real replenish- 
ment at sea operations. But this agreement assures that 
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this simulation was done correctly and gives the neces- 
sary background for the investigation and solution of 
the control problem. 

During the approach and departure of the tracking 
ship there is a risk of collision since at these posi- 
tions the lateral separation distance is decreased due 
to the interaction effects created from the pressure dis- 
tribution fields of the hulls of both ships. 

To be more specific, with the leading ship placed at 
(0,0) and having a 15 knots speed, i.e. the same speed 
with that of the space coordinate system, also v/ith the 
tracking ship placed at (—524 ft, 3) and having a An 
RPM greater propeller speed than that of the 15 knots 
ship the following can be deduced. 

a. During the approach the stern of the leading 
ship and the bow of the tracking ship tend 
to attract, while the lateral distance is 
decreased until the longitudinal separation 
distance is roughly —250 ft. Also both ships 
yaw in the same direction, but at different 
rates and values depending on their speed 
difference . 

b. For a decrease in longitudinal separation 
from —250 ft to zero ft the lateral separa- 
tion increases although the ships yaw to 
different directions and values. 

c. At zero longitudinal distance, i.e. exactly 
abeam position and 50 ft lateral distance 
both ships yaw negatively by the same amount 
while the lateral separation distance slightly 
increases throughout the time of recording. 

d. With both ships having the same speed of 15 
knots and being exactly abeam there is little 
change of the lateral distance throughout 
the time of recording although both yav; 
negatively . 
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e. During the departure from the exactly abeam 
position to until the longitudinal distance 
of roughly +350 ft. 

(1) A decrease of the lateral distance takes 
place . 

(2) Initially the bow of the leading ship 
tends to be attracted to the stern of 
the tracking ship, while both ships yaw 
negatively at different angles. 

(3) At roughly +250 ft longitudinal separa- 
tion distance the leading ship changes 
direction of yawing v/hile the tracking 
ship maintains negative yawing. This 
tends to bring the stern of the tracking 
ship towards the leading ship. 

f. For longitudinal separation distances greater 
than roughly 350 ft, the lateral distance 
increases while botli ships yavv/ positively 

at different angles and rates. 

One basic assumption was carried out throughout the 
UNREP simulation. That is both ships could be initially 
placed at any desired position. 

Fi'om the above discussion it is seen that the best 
position for station keeping and collision avoiding for 
tv;o ships of the same size is while keeping the same speed 
exactly abeam, although a certain amount of rudder has to 
be applied on both ships for course keeping. This is in 
agreement v;ith reference [1] . It is obvious that dujring 
the departure and approach the rudder to be applied must 
counteract both the interaction moment as well as create 
a yav; angle sufficient to counteract bot)i the attraction 
force and the rudder force. By these means collision 
should be avoided provided that there is enough initial 
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lateral separation distance so the time elapsed between 
the instant action is applied by the operator (helmsman) 
and rudder action takes place is adequate. It can be 
seen also that during the approach and departure both 
ships yaw from positive to negative angles and vice versa. 
This implies that in a short period of time the operator 
has to apply negative to positive rudder angle provided 
that the operator has been trained enough so that he knov;s 
the exact time at which rudder has to be applied . Ob- 
viously this is not easy. 

This study probably would be more effective using the 
interactive graphics available on modern computer systems. 

A matter of future investigation should be the imple- 
mentation of automatic controls in the hybrid sim.ulation 
as well as the inclusion of sea states. 
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VII . 



INVESTIGATION OF THE CONTROL PROBLEM 



The leading ship usually is responsible for course 
keeping and the tracking ship is responsible for both 
course and station keeping. These requircrnents imply two 
different controls for the course and the station or dis- 
tance keeping. 



Figure 99. Course keeping loop. 



A. COURSE KEEPING LOOP 

Figure 99 shows the control loop for course keeping. 



first order lag function. In the section B2 the time lag 
was estimated to be of’ the order of 8 sec. Hence the 
transfer function should be as follows: 




K r‘- 




G 
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The course keeping loop may not contain the yaw angle 

reference, since during the UNREP operation course 

changing is not desired but course keeping. Hence for 

the beginning at least of the simulation with automatic 

control included the reference yaw angle should be 

set equal to zero. The ships dynamics block contains the 

6 R 

transfer function G = — which can be found by Laplace 

transforming equations (V— 1) and solving for the above 

mentioned ratio. It should be noted here that both the 

transfer function G, and G„ have been built for the analog 

il) 2 ^ 

simulation of Phase I and II. Obviously in question must 
be the order of magnitude of the feedback gains and K. 
The values of and K can be designed by normal feedback 
control techniques. 

B. DISTANCE OR STATION KEEPING LOOP 

Figure 100 shows the distance keeping loop where 3^.^^ 
is the ordered lateral separation distance and 3 the actual 
lateral distance between the two ships. Note that this 




,e 



Figure 100. Distance keeping loop. 
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control loop contains the ships dynamics G2 including 
hybrid simulation for obtaining the lateral distance. 

In question must be whether the feedback function 
should contain the rate of change as v/ell as the change 
of the lateral distance. The steering control transfer 
function is again: 



Note that for the distance keeping loop the ships 
dynamics transfer function G2 involves hybrid simulation 
since transformation of coordinates is implied for the 
quantities being measured with respect to the space co- 
ordinates system. 
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COMPUIFR PROGRAM I A 



'k 

k 

* STATIC TEST FOR SHIP NO I 
-.JOB GALL225 
-iFCRTRAN LS,GO 

C SECTION CHECKING POT AND AMP VALUES FCR SHIP Ul 

DIPENSIDN POTAUOOI ,POTVAL(30) ,AMPAD(30) ,APPVAL(30) 
INTEGER POTADjAMPAC 
DATA N/26/ 

READ! 5, 100) (POT ADI 15 ,POTVAL( 1 ) , I 1 , N ) 

IOC FORMAT! A4, FI 0.4) 

DO 10 I^1,N 

10 CALL SETPQTIPOTADI I ) , POTVAL ( I ) ) 

WRITE! 6 ,101) !POTAD! I ) .POTVAL! I ) , I = 1,N) 

101 FORMAT! A4, FIO. 4) 

READ!5tl02) ! AMPAD ! I ) , 1 = 1 , 28 ) 

102 F0RMATIA4) 

CALL RESETIIOOO) 

DC 12 I=U20 

12 CALL SC AN! AMPAD! I ) , AMPVAL ! I ) ) 

WRITE! 6, 103) ! AMPAD! I ) , AMPVAL! I ) ,1=1,28) 

103 FORMAT! A4.F10. 6) 

DO 20 1=1 tN 

20 CALL SCAN! POTAD! I ) , POTVAL ! I ) ) 

WRITE! 6, 200) ! POTAD! I) , POTVAL! I) ,I=1,N5 
200 FORMAT! A4, FIO. 6) 

END 

-.LOAD XR,MAP 
-.DATA 



POOO 


0.1000 


POO 1 


0.0300 


P002 


0.0803 


P003 


0.0600 


POO 4 


0.0404 


P005 


0 .0052 


P0C6 


0.1660 


POO? 


0 . 0900 


POlO 


0.2379 


poll 


0.2303 


P012 


0.4717 


P013 


0.0250 


P 0 1 4 


0.1200 


P015 


0.5607 


P016 


0 . 1340 


PCI? 


0.1500 


P020 


0 ,C920 


P02 1 


0 .0097 


P022 


0.3576 


P053 


0.0905 


P054 


0.0500 


P059 


C. 1497 


P0900 


.0337 


P051 


0.3000 


P056 


0.0092 


POST 


0.4540 


AOOl 




A003 




A005 




A013 




A017 




COO 1 




D003 




D005 




0013 




0017 




A006 




AO 16 




A002 




AO 10 




AO 14 
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A022 
A203 
A 112 
A201 
A050 
A056 
A052 
A060 
A212 
A05 1 
A057 
D051 
D057 



COMPUTER PROGRAM IB 



* 

* STATIC TEST FOR SHIP NO 2 
•nJOB GAL 1225 
•nFORTRAK LS,GG 

C SECTION CHECKING POT AND AMP VALUES FCR SHIP 112 

DIMENSION POT AD ( 30 ) ,POTVAL ( 30 ) t AMPA0I30 ) , AMPVAL ( 30) 
INTEGER POTAD.AMPAO 
DATA N/22/ 

READ! 5 , 10 0 ) ( POTADC I ) , POTVALl I) , I = 1,N ) 

100 FORMAT! A4 ,F10. 4) 

DC 10 I = ltN 

10 CALL SETPOTI POTAD( I ) , POTVAL ( I ) ) 

WRITE! 6, 101) (POTADI I ) , POTVAL! I ) » 1 = 1 ,N) 

10 I FORMAT ( A4» F10.4) 

READI5.102) ! AMPAD! I ) , 1 = 1 ,27 ) • 

102 FORMAT I A4) 

CALL RESETIIOOO) 

DO 12 1=1,27 

12 CALL SC AN I AMPAD! I ), AMPVAL ( I ) ) 

WRITE { 6, 103 ) ! AMPAD! I ) , AMPV AL ! I) , I = 1 , 27 ) 

103 FORMATI A4,F10.6) 

DO 20 1=1, N 

20 CALL SCAN! POT AD! I ), POTVAL ! I ) ) 

WRITE (6 ,200) ! POTAD! I ) , POTVALl I) , 1 = 1 ,N) 

200 FCPMAT { A4, FL0.6) 

END 

-.LOAD XR,MAP 
•nDATA 



P0 2 4 


0.1000 


P025 


0.0300 


P026 


0.0803 


P027 


C.0600 


P030 


0 .0404 


P045 


0.C920 


P031 


0.0052 


P023 


0.0905 


P034 


0 . 1660 


P033 


0.0900 


P046 


0 .4717 


P032 


0.2379 


P035 


0.2383 


P047 


0 . 1497 


P036 


0.0097 


P03 7 


0.0250 


P040 


0 . 1200 


P041 


0 . 1340 


P042 


C. 1500 


P043 


0 . 5607 


P044 


0.3576 


P052 


0.5000 


A035 




D035 




A027 




D027 




A033 




D033 




A04 1 




D041 




A043 




C043 




A034 




A036 




A026 




A044 




A030 




A 04 2 




A 104 




A 06 5 




D065 
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A061 
D06 I 
A072 
A076 
A066 
A20A 
A21 1 
A214 
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COMPUTER PROGRAM III 



* LINEAR INTERPOLATION IN 2-DIMENSIONS ARRAY CF N-MOMENT 
AND Y-FORCE 
-.JOB GAL 1225 
-.FORTRAN LS,GO 

REAL NMAXt NDI M, NSCAL 
REAL NBARrNtNBARL 
DIMENSION Y( 30 ,30 ) ,N(30,30 J 
C STCRAG OF N ONO I ME NS 1 ONAL ARRAY FOR N AND Y 
DO 10 J=l,25 

1C READ! 5, 1001 ( Y ( J , I ) , I = 1 ,6 » 

100 F0RMAT(6E1 1.4> 

DO 12 J=l,25 

12 READ(5, 100) (NU, I ) ,1 = 1,6) 

WRITE(6,200) 

20C FORMAT!' ',‘THIS IS Y(J,I)‘,//) 

DO UJ=1,25 

WRITE (6,201) (Y( J, I ) , 1=1,6) 

201 FORMAT!' ' , 1 OX , 6 ( E 1 1 . 4 , 2X I , / / ) 

11 CONTINUE 

WRITE (6,300) 

30C FORMAT!' ','THIS IS N(J,I)',//) 

DO 13 J=l,25 

WRITE (6, 301) !N( J, I ) ,1 = 1,6) 

301 FORMAT!' • , lOX , 6 ( E 1 1 . 4 , 2X ) , / / ) 

13 CONTINUE 

C INTERPOLATION IN 2-DIM. ARRAY. PASS BETBAR AND ALPBAR 
C REMEMBER BETBAR = 1.0 CORRESPONDS BFTA = 100.0 FEET=?** 

C ALPBAR=1.0 CORRESPONDS ALPHA=550.0 FFET 
C ALSO BETBAR AND ALPBAR ARE PAST IN COMPUTER UNITS 
READ( 5, AOOIALPBAR, BFTBAR 
400 FCRMAT(2F11.6) 

W R n E ( 6 , 4 0 1 ) 

AOl FORMAT!' ',10X,'THIS IS AL PBAR-BETBAR ' , // ) 

WRITE (6 ,402) ALPBAR , BETBAR 

402 FORMAT!' ' , l.OX , 2 F 1 1 . 6 ) 

C RESCALING OF ALPHAB AR= ALP BAR AND BETABAR = BE TBAR 
ALPHA=550. 0"'-ALP8AR 
KRITE(6,403) 

403 FCPMAT(' ',10X,'THIS IS ALPHA',//) 

WPITE(6, 501) ALPHA 

BETBA=ABS( BETBAR) 

BETA=100.0*BETBA 
WRITE! 6 ,404) 

404 FORMAT!' ',10X,'THIS IS BETA',//) 

WRITE(6,501)BETA 
IF(BETA.LT.50.0) GO TO 803 
IFIBETA .GT .100.0) GO TO 801 
IF!ALPHA.LT. -550.0) GO TO 801 
IF!ALPHA.GT. 550.0) GO TO 801 

GG TO 800 
801 NEAR=0.0 
YBAR=0. C 
GC TO 802 

C SEARCH FOR IMDET I F 1 C AT I ON OF ARRAY ELEMENT 

C FIND CORRFSPGNDED Y AND N -LINEAR INTERPOLATION IN 2-DIMEN 
80G IE1=BETA/10.0 
IB=IBl-4 
BETAL=r Bi>no 
WTE=! BETA-BFTAl ) /lO.O 
IF! ALPHA.lt. 0.0) GO TO L 
IF! ALPHA. GT. 0.0) GO TO 2 
C ALPHA =0 
IA=13 
WTA=0 .0 
GC TO 3 

C ALPHA LESS THAN ZERO 
! ’ Al=- ALPHA/50.0 

ALPhAl = -IAr-'50 
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KTA=-( ALPHA- ALPHA I )/50.0 
I A=13-I AL 
GC TO 3 

C ALPHA GREATER THAN ZERO 

2 IAL=ALPHA/50.0 
I A=13+IAl 
ALPHAL=IAl*50 

WTA={ ALPHA-ALPHAl ) /50.0 

3 IE2=IB4 I 
IE(IfWEQ.6) IB2 = 6 
I A2 = ] A 1 

IF(IA.EQ.25) IA2=25 
WR1TE(6, 101) IBl , 1 B, !B2, I A1 

101 FCRMAT(' ' vlOX,' IB1 = ' , I3,5X,' IB=S I3,5X, 'IE2=S I3,5X,' 
WRITE{6,102) IA,iA2 

102 FCRMAK' ' , lOX , ‘ I A= ' , I 3 , 5 X , ' I A2 = M 3 , / / ) 

WRITE (6 ,202) BETAl ,WTB,ALPHA1 

202 FORMAT! > ' , ' B ET A 1 ' , E 14 . 8 , 2X , ‘ KT B ■ , E 14. 8 , 2X , • ALPHA 1' ,E1 

WR1TF(6,203) WTA 

203 FORMAT!' ' , ' WTA= ' , E 14 . 8 , // ) 

I F( ALPHA.lt. 0.0) IA2=IA-1 
IFUA.EG.l) IA2^1 
IF(WTA.EQ.0.0)IA2=1A 
IF(WTB.EQ.O.O) IB2=I3 
NGAR1=N(IA,IB2) 

WRITE {6, 602) 

602 FORMAT!' ','THIS IS NBARl',//) 

WR1TE!6,502) NBARl, IA,1B2 
NPAR=N! I A, IB) 

KRITF (6 ,603) 

603 FORMAT!' ','THIS IS NBAR',//) 

WRITE ( 6. 50 2) NBAR, I A, IB 
BARN1= ! NBARl-NBAR) *WTB+NBAR 
WRITE! 6,604) 

604 FORMAT!' ','THIS IS BARNl',//) 

WRITE(6,501) BARNl 

NBARI=N( IA2, IB2) 

WRITE 16,605) 

605 FORMAT!' ','THIS IS NBARl',//) 

WRITE!6,502) NB AR 1 , I A2 , I B2 
NBAR=N! IA2,IS) 

WRITE! 6,606) 

606 FORMAT!' ','THIS IS NBAR',//) 

WRITE! 6, 502) NBAR, I A2, IB 
BARN2:=1 NBARl-NBAR) >^WTBi-NBAR 
WRITE (6,607) 

607 FORMAT!' ','THIS IS BARN2',//) 

WRITE! 6,501) BARN2 

NOAR= ! BAPN2-BARN1 )«WTA-^BARNl 
WRITE! 6,500) 

500 FORMAT!' ','THIS IS FINAL NBAR',//) 

WRITE!6,501) NBAR 

501 FORMAT!' ' , lOX , E 14 . 8 , // ) 

YBARl-Y ! I A, IB2) 

WRITE !6 ,608) 

608 FORMAT! ‘ ','THIS IS YBARl',//) 

WRITE!6,502) YBARl, IA,IB2 
YBAR=Y! I A, IB) 

WRITF!6,605) 

609 FORMAT!' ‘,'THIS IS YBAR',//) 

WRITE! 6,502) YBAR, lA, IB 
BARY1= ! YBAR 1-YbAR ) ’^WTB + YBAR 
WRITE (6 ,610) 

61C FORMAT!' ','THIS IS BARYl',//) 

WRITE !6,60l)BARYl 
YBAR1=Y!IA2, IB2) 

WRITE! 6, 61 1) 

611 FORMAT!' ','THIS IS YBARl’,//) 

WPITE!5,502) Y3A R 1 , I A2 , I B2 
YBAR=-Y! IA2, IB ) 

W P I T r ! 6 , 6 1 2 ) 

612 FORMAT! ‘ ‘,‘THIS IS YBAR',//) 
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WRITE(6,502) YBAR,IA2,1B 
BARY2= { YBARL-YBAR ) ’^WTB+YBAR 
lvRITE(6 ,613) 

612 FQRMAT(* ',«THIS IS BARY2*,//) 

KRITE(6t60L) BARY2 

YBAR={ 3ARY2-BARYI ) >'-WTA+BARYl 

WRITE (6 ,600) 

600 FCRMAK* ','THIS IS FINAL YBAR',//) 

KRITE(6,601) YBAR 

601 FGRMAK' ' , 1 0 X , E 1 4 . B , / / ) 

5C2 PORMAK » • , I OX , E 16 . 3 , 5X , ' R 0W = ' , I 3 , IX , ‘ C OL = ' , 1 3 ) 

C REREMBER SIMULATION IS DONE IN COMPUTER UNITS SC NO NEED T 
DIVIDE BY 100.0 SINCE SUBR. CAC CONVERTS BY ITSELF IN 
GO TO 2000 

802 NEAR=66 .5E-05 
WPITE(6,500) 

WRITE(6,501) NBAR 
YEAR=90 .0E~05 
WRITE (6 ,600) 

WRITE(6,601) YBAR 
GO TO 2000 

802 KRITE(6,300) 

WRITE (6, 501) NBAR 
WRITE (6 ,600) 

WRITE (6, 601) YBAR 
2000 END 
LCAC XRfMAP 
DATA 
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-0. 1600E-03-r0 .1263E-03-0, 1160E-03-0. 1 058E-03-0 . <5 75 3E-0A-0 .<30 
-0.2600E-03-0 .2112E-03-0 .1897E-03-0, 1681E-03-0. 15C8E-03-C. 13 
-0. 3 2COF-03-0.2 762E-03-Oc243 3E-03-0.2 106E-03-0 . 1861E-G3-0 . 16 
-0. 38C0E-03-0 .3225P-03-0.2793F-03-0. 2362E-03-0 .20 16E-03-0 . 17 
-0.4000E-03-0 .3 337E-03-0 .2840E-03-0 .2542E-03-0 . 1945E-03-0 . 15 
-0. 2 50 0E-03-0 .2 870E-03-0.23 97F-03-0. 19 25E-03--0 . 1547E-03-0 . 12 
-0. 26 00E-02-0 .2079E-03-0 . 16 89E-03-0 . 1 2 99E-03- 0 . 9 8 63E-04-- C . 70 
“0. 1300E-C3-C.9384E-04-0 .667 lE-04-0 .3959E-04-0 .1789E-04 0.20 
0.4000E-04 0.5507E-04 0.6627E-04 0.7767E-04 0.8671F-04 C.95 
0.2700E-03 0.2508E-03 0 .2364E-03 0.2221E-03 0.2 1C5E-03 0.20 
O.5C0CF-03 0.4507E-03 0.4137E-03 0.3767E-03 0.3471E-03 0.32 
0.7C00E-03 0.6225E-03 0.5643E-03 0.5062E-03 0.4596E-03 0.41 
0.8500E-03 0.7500F-03 0 .6750E-03 0.6000E-03 0.54C0E-'03 0.48 
C.9C00E-03 0.7844F-03 0.6977E’-03 0.6110E-03 0.5416E-03 0.47 
0.8200E-03 0.7159E-03 0.6378E-03 0.5597E-03 0.4972E-03 0.44 
0.6250E-03 C.5579F-03 0.5075E-03 0.4572E-03 0.4169E-03 0.38 
0.4500F-03 0.4C89E-03 0.3781E-03 0.3473E-03 0.2226E-03 0.30 
0.3000E-03 0.2781E-03 0.2616E-03 0.2452E-03 0.2221E-03 0.22 
0.1700E-03 0.1673E-03 0.1652E-03 0.1632E-03 0.1615E-03 0.16 
0.6000E-04 0.6822E-04 0.7438E-04 0.8055E-04 0.8548E-04 C.90 
-0.2000E-04-0 .3562E-'05 0 .8767E-05 0.2 110E-04 0.3C96E-04 0.40 
-G. 5C00E-04-0 .3493E-04-0.2363E-04-0. 1233 E-04-0 .32 88 E-05 C .50 
-0.6C00F-04-'0 .49 04E-'04-0 .4082E- 04-0 . 3 2 60E-04- 0 . 2 6C3E-04- C . 20 
-0. 5C00 E-04-0 .46 99E-04-0 .447 3 E-04-0 .4 247 E-04-0 .4 06 5E -04 - 0 . 39 
-C. lCOOE-04-0 .2096E-04-0.29L8E-04-0.3740E-04-0 .4297E-04-C.50 



* 

* 

* 

* 

0.6000E-04 
0. llOCE-03 
0. 1390F-03 
0 . 133C:-03 
0. 14C0E-03 
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0 .6000E-04 
0, 1036F-03 
0. 1278E-03 
0 . 13 49E-03 
0.I197E-03 



0 .6000E-04 
0 .9T86F-04 
0.U76F-03 
0 .1186P-03 
0 . I012E-03 



0.60CCE-04 
0 .9280E-04 
0. 1087E-03 
0. lC^'^p-03 
0.8497E-C4 



0.6CC0E-04 0.60 
0.8890F-04 0.85 
0.LC19E-03 0.95 
r 

6. 7 2 49 E -04 0.60 
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0.7000E- 
0.5000E- 
0. I ?5(3E- 
0.3000E- 
^ ICOE- 
A420E- 
4350E- 
3 74 0E' 
28C0E- 
l 300E- 
O.-^ECCE- 
0.'^ ?C0E- 
0. 1) OOE- 
0. L^OOE- 
0. 1320E- 
0. 12C0E- 
0. 1020E- 
0.7500E- 
0. A5C0F- 
0.9000E- 



OA 0 
04-0 
C3-0 
03-0 
03-0 
03-0 
03-0 
03-0 
03-0 
03-0 
03-0 
04 0 
03 0 
03 0 
03 0 
03 0 

03 0 

04 0 

04 0 

05 0 



.5728E 
,46 1 8F 
• 1 536E 
.2598E 
.3520E 
.3S80E 
.3785E 
.33 OOE 
.2520E 
. 1373E 
.43 73F 
.3640^ 
.9728E 
. I 247E 
. n 93E 
. lObOE 
. 875CE 
.6355E 
.3864E 
. 10 53E 



-04 0 
-04-0 
-03-0 
-03-0 
-03-0 
-03-0 
-03-0 
-03-0 
-03-0 
-03-0 
-04-0 
-04 0 
-04 0 
-03 0 
-03 0 
-03 0 
-04 0 
-04 0 
-04 0 
-04 0 



.4 572E- 
.4272E- 
. 13 42F- 
.2238E- 
.300 8E- 
.3294E- 
.3272E- 
.2900F- 
.2266F- 
. 12 57E- 
,4257E- 
.3132E- 
.8572E- 
. 1.1 C9E- 
.L077E- 
.9329E- 
.7432E- 
.5315E- 
.3286E- 
. 1 19 lE- 



04 0 
04-0 
03-0 
03-0 
03-0 
03-0 
03-0 
03-0 
03-0 

03- 0 

04- 0 
•04 0 
04 0 
■03 0 

03 0 

04 0 
04 0 
■04 0 
■04 0 
04 0 



.3561E- 
. 3968E- 
. I 172E- 
. 19 13E- 
.2552 E- 
.2 824E- 
.2823E- 
.2550E- 
.2043E- 
. 1 156E- 
.4 156E- 
.2687E- 
.7561E- 
.9873E- 
.9761E- 
.8217E- 
. 6279E- 
.44C5E- 
.2780E- 
. 13L3E- 



04 0 
04-0 
03-0 
03-0 
03-0 
03-0 
03-0 
03-0 
03-0 

03- 0 

04- 0 
•04 0 
04 0 
04 0 
04 0 
•0 4 0 
04 0 
04 0 
■04 0 
04 0 



. 2780E-04 

.3734E-04- 

. 1041 E-03- 

. 1667E-03- 

.22C1E-03- 

.2462E-03- 

.2476E-03- 

.2280E-03- 

. 1872E-03- 

. 1C78E-03- 

.4078E-04- 

.2343E-04 

.6780E-04 

.8936E-04 

.0980E-04 

.7358E-04 

.5390E-04 

.37C2E-04 

. 2 390 E -04 

. 14C6E-04 



0.20 
0.35 
0 .91 
0.14 
C. 18 
0.21 
0.21 
0.20 
0.17 
0.10 
0 .40 
0.20 
0.60 
0 .80 
0.82 
C.65 
0.45 
0.30 
0.20 
0.15 
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COMPUTER PROGRAM IVA 



* COORDINATES TRANSFORMATION FOR SHIP NO 1 
-^JCB GAL1225 
-^ASSIGN 6=NQ 
-.FORTRAN LS,GO 

C SECTION CHECKING POT AND AMP VALUES FOR SHIP #1 

D I MENS I ON POT AD (30 ) , POTVAL (301, AM P AD ( 30 ) , AMFVAL (30) 
INTEGER PGTAO.AMPAD 
DATA IFLAG/0/, IFLAGL/0/ 

DATA Nl/24/ 

READ! 5, 10 0 ) ( POT AD (I ) , POTVAL { I) v 1 = 1 ,NU 

100 FORMAT ( AA , FIO .A) 

DO 10 I=1,M 

10 CALL SETPOT(POTAD( I) , POTVAL! I ) ) 

WRITE (6, 101) (POTACd ),P0TVAL(1 ) ,I = 1,N1) 

ICl FORMAT! AA, FIO. 4) 

READ(5,102) ( AMPAD( I) , I=U29) 

102 FORMAT (A4) 

CALL RESET! 1000) 

DC 12 1=1,29 

12 CALL SCAN( AMPADI I 5 , AMPVAL ( I ) ) 

WRITE! 6, 103) (AM PAD (I ) ,AMPVAL(I ) ,1=1,29) 

103 FCPMAT(A4, FI0.6) 

DC 20 I=l,M 

2C CALL SCAN(POTAD( I ) tPOTVALd ) ) 

WR ITE (6 ,200) ( POT AD! I ) , POTVAL! I ) , 1 = 1, N1 ) 

2CC FORMAT! A4,F10. 6) 

C WHFN READY 

30 OUTPUT (101 ) 'READY ' 

C TRANSFORMATION OF COORDINATES FOR SHIP HI 

11 CALL COMPUTE 

111 CALL ADK!0,U1BAR, 1 ,V1BAR,2,PSIBAR ) 

WRITE! 6,201) Ul B AR , V 1 B AR , P S I BAR 

20 1 FORMAT ( * « , 'U IBAR ' , E14.8, 2X, ' V IBAR* ,E14.8,2X,'PSIBAR', 

C RESCALING OF PSIBAR 

PSIMAX=0. 26 
PSI=PSI MAX^PSIBAR 
A = SIN! PSI ) 

B=COS ( PSI ) 

WR1TE!6,106) A,B 

IC6 FORMAT!' ‘ ‘ A = ' , E 1 4 . 8 , 2X , ' B= ‘ , E 1 4 . 8 ) 

WRITE (6,202) PS I 

202 FORMAT! ' ' , 'PSI = ' , E14.8) 

XC 1D0T=U1 PAR^'-B--V1BAR*A 
YC1D0T=U1 BAR-^A+V1PAR*B 
WPITE!6,107) XQ1D0T,Y01D0T 

107 FORMAT! • ' , ' XQ1D0T = • , P 14 . 8 , 2X , ' YO ID0T= ' , E 14 . 8 ) 

CALL DACd ,X01D0T , 2,Y01D0T ) 

C MANUAL DIGITAL SWITCH'D' AND *1' 

15 IF!TFST d ) .GT .0) GO TO 16 
CALL HOLD 

17 IF!TEST d) .LT.O) GO TO 17 
CALL COMPUTE 

16 IF!TEST (3) .GT.O) GO TO 111 
CALL RESETdOOO) 

PAUSE 
GO TO 11 
END 

-.LOAD XR.MAP 
-.DATA 

POCC C.3000 
P025 0.2860 

P002 0.0803 

P023 0.2860 

P0C4 0,0404 
P005 0.0052 

P0C6 0.1660 
POlO 0.2T79 
Poll 0.2383 



207 




■■■As 



P012 


0 . A717 


P013 


0.0230 


P015 


C.5607 


P016 


C. 1340 


P017 


0 .000 


P020 


0.0920 


P02 1 


0.C097 


P022 


0.3576 


P052 


•C.C905 


P027 


0. 1750 


P055 


0 . 1497 


P05G 


0.0337 


P05 1 


0.3000 


P056 


0.0092 


P057 


C .4540 


AGO 1 




DOOL 




A002 




A003 




D003 




A 201 




A006 




A005 




D0C5 




AGIO 




AOIA 




A016 




A0L3 




D013 




A017 




D017 




A022 




A203 




A 112 




A060 




A053 




D053 




A206 




A200 




A051 




D05 1 




A057 




D057 




A202 





COMPUTER PROGRAM IVB 



* COORDINATES TRANSFORMATION FOR SHIP NO 2 
-'JOB GAL1225 
-.ASSIGN 6 = N0 
-.FORTRAN LS,GG 

C SECTION CHECKING POT AND AMP VALUES FOR SHIP H? 

DIMENSION POT ADI 30 ) ,POTVAL ( 30 ) , AMP AD { 30 J , AMPVALOO) 
INTEGER POTADfAMPAC 
DATA I FLAG/0/, IFLAGl/O/ 

DATA N2/22/ 

READ! 5 , 100) (POT ADI I ) ,POTVALI I ) , I = 1,N2) 

100 FORMATI AA, FIO.A) 

DC LO I=LtN2 

1C CALL SETPCTI POTADII ) , PQTVAL I I) ) 

WRITE (6, 10 1) IPOTADII ) ,POTVAL( I ) ,I = L,N2) 

101 FCPKATl AA,Fi0.4) 

READ(5t102) IAMPADII ) ,1=1,27) 

102 FCPMATIAA) 

CALL RESETIIOOO) 

DC 12 1=1,27 

12 CALL SCAN I AMPADI I ) , AMPVAL I I ) ) 

WRITE (6, 103) IAMPADII ) , AM PV AL { I ) , I = 1 , 27 ) 

102 FORMAT! A4,F10. 6) 

DC 20 1=1, N2 

2C CALL SCAN! POTAD I I ) , POTVAL ( I 1 ) 

WRITE I 6 , 20 0) I POT ADI I ) , POTVALI I ) , I = 1,N2) 

20C FORMAT! A4,F10. 6) 

C WHEN READY 

20 OUTPUT { 10 I ) 'READY ' 

C TRANSFORMATION OF COORDINATES FOR SHIP #2 
11 CALL COMPUTE 

112 CALL A0K(4,U2BAR, 5,V2BAR,6,PS2BARJ 
C RESCALING OF PS2BAR 
PS2MAX=0. 26 
PS2=PS2MAX*PS2BAR 
A2=SIN( PS2 ) 

B2=C0SI PS2) 

XC2D0T=U2EAR*B2-V2BAR* A2 
YC2D0T=U2BAR-A2+V2EAP*E2 
CALL DACI 5,Y02D0T, 6,X02D0T) 

C MANUAL DIGITAL SWITCH '2' AND *3' 



18 


IFITESTI2) .GT.O) 
CALL HOLD 


GO 


TO 


19 


21 


I FITEST 12 ) .LT.O) 
CALL COMPUTE 


GO 


TO 


21 


19 


IFITEST 14) .GT.O) 

CALL RESETIIOOO) 

PAUSE 

GC TO 11 

END 


GO 


TO 


112 


-'LOAD 

-.DATA 


XR,MAP 








P02A 


C.2000 








P09A 


0.2860 








P026 


0.0803 








POCl 


C. 1750 








P02C 


C .C404 








P04 5 


0.C920 








P03 1 


0.C052 








P023 


0.0905 








P03A 


0 . 1660 








P0C3 


0.2860 








P046 


0 .A717 








P032 


0.2379 








P035 


0.2383 








P047 


0 . 1A97 








-P036 


0 .C097 








P03? 


C , C250 








P04 1 


0.1240 
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P042 

P0A3 

P007 
P014 
A035 
D035 
A026 
A027 
D027 
AOAA 
A034 
A033 
D023 
A030 
A036 
AOA 1 
COA I 
AOA3 
DOA3 
AOA2 
AlOA 
A052 
A007 
DOC7 
A205 
A20A 
A065 
D065 
A21A 
A061 
D06 I 



C.OOOO 
C.5607 
0 .2576 
0.0337 
0.0092 
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COMPUTER PROGRAM VA 



PHASE I . INTERACTION FORCE AND MOMENT ARE APPLIED TO THE 
LEADING SHIPISHIP MODEL NO 2) 

-.JOB GAL1221) 

-^FCRTRAN LS,GO 

REAL NMAXjNDIM^NSCAL 
REAL NBAR,N,NRARl 
DIMENSION Y(30,30 ) tN(30,30) 

DIMENSION POT AD ( 70) , POTVAL ( 70) , AMP ADI 70) , AMFVAL (70 ) 
INTEGER POTAD,AMPAC 
DATA I FLAG/0/, IFLAGl/0/ 

DATA NO/48/ 

C STORAGE OF NCNDI MENS I GNAL ARRAY FOR NAND Y 
DO 10 J = 1 , 25 

10 READ(5,100) I YI J, I ) , 1 = 1,6) 

100 FCPMATI6E11.4) 

DC 12 J=1 , 25 

12 READI5, 100) (NU, I), 1 = 1, 6) 

X KPITE(6,200) 

200 FORMAT!’ ',«THIS IS Y(J,I)’,//) 

DC 11J=1,25 

X NRITE { 6,201 ) I Y (J , I ) , I =1 , 6 ) 

201 FORMAT! • ‘ , lOX , 6 ( E 1 1 . 4 , 2X ) , // ) 

11 CCNTINUE 

X NRITEI 6,300) 

300 FORMAT! • «,'THIS IS N!J, I )’,//) 

DC 13 J=l,25 

X NRITE(6,301) I N IJ , I ) , I = 1 , 6 ) 

301 FORMAT!* ‘ , lOX , 6 ( £ 1 1 . 4 , 2X ) , / / ) 

13 CCNTINUE 

C SECTION CHECKING POT AND AMP VALUES FOR SHIP Ul AND #2 
RE AD! 5, 1005) (POT AD! I ) , POTVAL! I I , 1 = 1, NO) 

1005 FORMAT! A4 ,F10. 4) 

DC 10 5 1 = 1, NO 

105 CALL SETPOT! POTAD! I) , POTVAL! I ) ) 

X NRITE! 6 , 1015) ! POT AD ( I ) , PO T V AL ! I ) , 1= 1 , NO ) 

1015 FORMAT ( A4,F10. 4) 

RE AD (5, 1025) ! AM P A D ! I ) , I = 1 , 66 ) 

1025 FORMAT! A4) 

CALL RESETIIOOO) 

DC 125 1=1,66 

125 CALL SCAN ( AMPAD! I ) , AMPVAL (I ) ) 

X WPrrE(6,103) (AMPAD! I ) jAMPVAL! I ) ,1 = 1,66) 

1C2 FORMAT! A4,F10. 6) 

DO 20 1=1, NO 

20 CALL SCANIPOTAD! I ) ,PQTVAL ! I ) ) 

X NR I TE (6,2005) (POTAD! I ) , POTVAL! I ) ,1=1, NO) 

2005 FORMAT! A4,F10. 6) 

C WHEN READY 

30 OUTPUT! 101 ) ‘READY’ 

C TRANSFORMATION OF COORDINATES FOR SHIP ^1 AND U2 
115 CALL COMPUTE 

111 CALL ADK! 0,U1BAR, 1 ,V1BAR,2,PSIBAR,4,U2BAR,5,V2BAR,6,PS 
6R ,3, ALP BAR ) 

C RESCALING OF PSIBAR AND PS2BAR 
PSIMAX=0.26 
PS2MAX=0.26 
PS1=PS1MAX=«PSIBAR 
A-SIN!PSI ) 

B = C '' S ! P S I ) 

X010QT = U1BAR-^B-V1BAR'''A 
YC1D0T=UIBAR*A+V1BAR*B 
PS2 = PS2MAX=^PS2BAR 
A2=SIN!PS2) 

B2 = C0S ! PS2 ) 

XG2mT = U?BARAB2-V2PAR>''A2 
YC2D0T=U2BAR*A2+V2BAR*B2 

C INTERECLAT lO:; !fl 2--DI''’. A P P AY .. A S S PET'^A'^ AN'r^ ALP^AR 
C REMEMEER BE-iBAR=l.0 CORRESPONDS bETA=lOO.O FtET>:''-^* 
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C ALFBAR=1.0 CORRESPONDS ALPHA=550.0 FEET 
C ALSO BETBAR AND ALPBAR ARE PAST IN CCRPUTER UMTS 
X WRITE(6,A0L) 

AOl FORMAT!* ',).OX,‘ThIS IS AL P B AR- B ET BAR ‘ , / / ) 

X VjRITE{6,A02) ALPBAR , BETBAR 

402 FCPMATI' « , lOX , 2 F I I » 6 ) 

C RESCALING OF ALPHABAR= ALPBAR AND BET A E AR=B ET B A R 
ALPHA=5E0. C^ALPBAR 
X WRITE(6v403) 

403 FORMAT!' ',10X,'THIS IS ALPHA',//) 

X WRITE(6,50!.)ALPHA 

BETA=100.0*BETBAR 
X RR1TE!6,404) 

404 FORMAT! • *vlOX,'THIS IS BETA',//) 

X TRITE! 6 , 501 ) BETA 

I F! BETA .LT .50.0) OUTPUTUOl) BE T A ; NSC AL = 0 . 0 ; YSC AL = 0 . 0 ; 
IF ! BETA.GT.IOO. 0 ) OUTPUT! 101 ) B'^T A ; NSCAL = 0 . 0 ; VSC AL‘-= 0 . 0 *,G 
I F! ALPHA.GT.550.0)0UTPUT( 1 0 1 ) AL PH A ; NSC A L = 0 . C ; YS CAL = 0 .0 
IF! ALPHA.lt. -550 . 0 .'OUTPUT! 10 U A L PH A ; NSC AL = 0 .0 ; Y SC A L = 0 . 
C SEARCH FOR I NOET I F I CAT I ON OF ARRAY ELEMENT 

C FIND CORRESPONDED Y AND N -LINEAR INTERPOLATION IN 2-DIMEN 
I E l = BETA/ 10. 0 
IB=lBl-4 
BETAl=l Bl*10 
WTE=! BETA-BET A I ) / lO.O 
IF!ALPHA.LT.O.O) GO TO I 
IF! ALPHA.GT.0.0) GO TO 2 
C ALPHA =0 
IA=I3 
WTA=0.0 
GO TO 3 

C ALPHA LESS THAN ZERO 

1 1 Al=-ALPHA/50.0 
ALPHA1=- I Al-^=50 

kTA=-! ALPHA-ALPHA! ) /50. 0 
IA=13-I Al 
GO TO 3 

C ALPHA GREATER THAN ZERO 

2 IAl=ALPHA/50.0 
I A = 13 + I A1 
ALPHA1=IAI=^50 

WTA=! ALPHA-ALPHA 1 ) /50 .0 
2 IE2 = IB4-1 

IFIIB.EC.6) IB2=6 
I A2=IA+1 

IF!1A.EG.25) IA2=25 

X TRITE 16,101) IBl, 18, IB2,IAl 

101 FORMAT ! ' ' , 1 OX , ' I B 1= ' , I 3 , 5X , ' I B= ' , 1 3 , 5 X , ' I B 2= ' , I 3 , 5X , ' 

X T'RITF! A , 102) IA,IA2 

102 FORMAT!' ' , LOX , ‘ I A = ' , I 3 , 5X , • I A 2= • , I 3 , / / ) 

X TR1TE!6,202) BETA 1 , WTB , AL PHAl 

20 2 FORMAT !' ','BETA1',E14.8,2X,'WTB',E14.8,2X, 'ALPHA I ' , El 

X WRITE 16,203) WTA 

2C3 FORMAT!' * , ' WTA = ' , E 14 . 8 , / / ) 

IF! ALPHA.lt. 0.0) IA2=IA-1 
IFUA.EC.l) IA2=1 
IF! WTA. EQ.C.O) I A2=IA 
IF !WTB.EQ.O.O) IB2=1B 
NEAR1=N ! I A, 1 B2) 

X TRITE!6,602) 

602 FCPHATI' ','THIS IS NBARl',//) 

X V^RITE!6,502) NBARl, IA,IB2 

NBAR=N! lA, IB) 

X WRITE!6,603) 

603 FORMAT!' ','THIS IS NBAR',//) 

X KRITE!6,502) NBAR, !A, IB 

6APN1= ! NBAR1-NBAR)"«WT E+NBAR 
X TRITE! 6,604) 

604 FORMAT!' ','THIS IS BARN!',//) 

X TRITE! 6,501) BARN I 

NEAP 1 = N ! ' A ?. I B2 ) 

X WRITE! 6, 605) 
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LC5 FCRKAK' ',‘THIS IS NDARl',//) 

X WRITE(6,502) NB AR 1 , 1 A ? , I B2 

l\'fiAR = N( 1A2. , IB) 

X URITE(6,606) 

606 FORMAT!' '.■'THIS IS NBAR',//) 

X URITE(6,502) NBAR,IA2,IB 

BARK2= ( KBARl-NBAR) Bi-NBAR 
X KRITE(6,607) 

607 FORMAT! ‘ ‘,'THIS IS BAPN2',//) 

X RRITL! 6,501) BARN2 

K'PAR= ! BARIM2-BARM1 ) ^WTA + BARNl 
X V«RrrE(6 ,500) 

50C FORMAT!' ','THIS IS FINAL NBAR',//) 

X RRITE!6,501) NBAR 

501 FORMAT!' ' , lOX , E 1 4 . 8 , / / ) 

YBAR1 = Y ! lA, I B2) 

X URITE(6,608) 

608 FORMAT!' ','THIS IS YBARl',//) 

X WRITE!6,502) YBARl, IA,IB2 

YBAR=Y( I A, IB) 

X WRITE!6,60S) 

609 FORMAT!' ','THIS IS YBAR',//) 

X R, RITE! 6,502) YBAR, lA, IB 

B ARY1= ! YBAR 1-YBAR ) *WT B+YBAR 
X WRITE!6,610) 

61C FORMAT! • ‘,‘THIS IS BARYl',//) 

X WRITE !6 ,601) BARYl 

YBAR1=--Y ! I A2, I B2 ) 

X WRITE! 6,611) 

611 FORMAT!' ','THIS IS YBARl',//) 

X WRITE!6,502) YBAR 1 , I A2 , I B2 

YBAR=Y! IA2, IB ) 

X WRITE! 6, 6 12) 

612 FORMAT!' ',‘THIS IS YBAR',//) 

X WRITE!6,502) YBAR,IA2,IR 

BARY2=! YBAR 1-YBAR ) >^WTe-!-YBAR 
X WRITE 16,613) 

613 FORMAT!' ‘,'THlS IS BARY2',//) 

X WRITE!6,601> BARY2 

YBAR=! BARY2-0ARYI )-WTA+BARYl 
X WRITE! 6, 600) 

6CC FORMAT!' ','THIS IS FINAL YBAR',//) 

X WRITE! 6, 601) YBAR 

601 FORMAT!' ' , lOX , E 14 8 , / / ) 

5C2 FORMAT!' ' , lOX , E 1 4 , 8 , 5 X , ' RO W= ' , I 3 , IX , ' COL = ' , 1 3 ) 

C REMEMBER SIMULATION IS DONE IN COMPUTER UNITS SC NO NEED T 
DIVIDE BY lOOcO SINCE SUBR. DAC CONVERTS BY ITSELF IN 
C MAGNITLDE SCALING OF N 
NMAX=4. 195E+07 
NDIM=NBAR*!9.428E+10) 

NSCAL=-NDIM/NMAX 
C MAGNITUDE SCALING OF Y 
YMAX=1 . 6C8E+05 
YDIM=YBAR=f( l7.87E-f 07) 

YSCAL=YDIM/YMAX 

X WRITE!6,900) NSCAL,YSCAL 

900 FORMAT!' ' , ' NSCAL= ' , E 14 .6, 5X, ' YSCAL= ' , E 14 . 6 , // ) 

X WRITE! 6, 901) NDI.M,YDIM 

901 FORMAT!' ‘ ' N DI M = ' , E 1 4 . 6 , 5 X , ' YD I M= • , E 1 4 . 6 , / / ) 

X WRITE!6,5001 ) XOIDOT, YOlOOT 

5CC1 FORMAT!' ' , • XOIDOT = ' , E 14. 6 , ' YOl DCT= • , E 1 4 . 6 , // ) 

X WRITE! 6,5000) XQ2DCT, Y02D0T 

■ 5000 FORMAT!' ' , ' X02D0T= ' , E 1 4 . 6 , • Y02 DOT= ' , E 1 4 . 6 , // ) 

C TRUNK LINES T426!7) FOR A076,T427!8) FCR A07 2 , TA22 ! 3 ), FOR 
C***TA23!4) FOR A207 

804 CALL DAC!1,XG1DOT,2,Y01DOT,5,Y02COT,6,X02CCT,3,NSCAL,4 
C MANUAL DIGITAL SWITCH'D' AND '1' 

15 IF!TEST!I) .GT.O) GC TO 16 
CALL HOLD 

17 IF!TEST ! U .LT.O) GC TO 17 
C/LL CO/'r^UT'^ 

16 1 F!TEST!:iJ .GT.O) GC TO 111 
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i 



CALL RESET( 1000) 

PAUSE 

GO TO 115 

END 

-iLCAC XPtMAP 
-iDATA 
❖ 

* 



THIS IS Y ARRAY 



i.': 

A 



-0. 1 AOOE- 
^0.2AC0E- 
-0. 220CE- 

-o.seooE- 

-0. ACCOE- 
-0. 3500E- 
“0. 2600E- 
-0. 1300E- 
0. AOCOF- 
0.2700E- 
0. 5C0CE- 
0.7COOE' 
0.85C0F- 
0. GCOOF- 
O.REOOE- 
0. A25CF- 
0. A5C0E- 
0.3C00F- 
0. 1 7C0E- 
0.6C00E- 
-0. 20COE- 
-0.5C00E- 
-0. 6C00F- 
-C. 5000L- 
-0. iOOOF- 



03-0 
03-0 
03-0 
0 3-0 
03-0 
03-0 
03-0 

03- C 
OA 0 
03 0 
03 0 
03 0 
■03 0 
03 0 
03 0 
C3 C 
03 0 
■03 0 

03 0 

04 0 

04- 0 
04-0 
04-0 
OA-0 
04-0 



. 1263E- 
.21 12E- 
.2762E- 
.3225E- 
.3337E- 
.2870E- 
.2079E- 
.9384E- 
.5507E- 
.2508E- 
.4507F- 
.6225E- 
.7500E- 
.7844E- 
.7159F- 
.5579E- 

.4089F- 

.2781E- 
. 1673E- 
.6822E- 
.3562E- 
.3493F- 
.4904E- 
.46 99F- 
.2096E' 



03- 

03- 

03- 

03- 

03- 

03- 

03- 

04- 
04 
03 
03 
03 
■03 
03 
03 
03 

03 
■03 
■03 

04 
•05 
■04- 
■04- 
04- 
04- 



0.1160E- 
0 .1897E- 
0 .24 33E- 
0 .2793E- 
0 .2840E- 
0 .2397E- 
0 . 1689E- 
0 .6671E- 
0 .6637E- 
0 .2364E- 
0 .4137E- 
0 .5643E- 
0 .67 50E- 
0.6977E- 
0 .6378F- 
0 .5075E- 
0 .3781E- 
0 .2616F- 
0.1652E- 
0 .7438E- 
0 .876 7F- 
0 .2363E- 
0 .4082E- 
0 .A473E- 
0.2918E- 



03-0 

03-0 

03-0 

03-0 

03-0 

03-0 

03- 0 

04- 0 
04 0 
03 0 
03 0 
03 0 
03 0 
03 0 
03 0 
03 0 
03 0 
03 0 

03 0 

04 0 

05 0 
OA -0 
04-0 
04-0 
04-0 



. I058E-03- 
. 1681E-03- 
.2 1C4E-03- 
. 2 3C2E-03- 
.2342F-03- 
. 1925E-03- 
. 1299E-03- 
.3959E-04- 
.7767E-04 
.2221E-03 
.3767E-03 
. 5062E-03 
.6000E-03 
.6 llOE-03 
.5597E-03 
.4572F-03 
.3473F-03 
. 2A52E-03 
. 1632E-03 
. 8055E-04 
.2 i lOF-04 
. 1233E-04- 
.32COE-04- 
.4247E-04- 
.374CE-04- 



9753E 
1 5C8E 
18A1E 
2C 16E 
19A5E 
1547E 
9863E 
17898 
0 .8671E 
0.21C55 
0 . 3 A 7 1 E 
0 .A 596F 
0 .5AC0E 
0 .5A16E 
0.A973E 
0 .A 169E 
0 . 3226E 
0 .2321E 
0 . 1615E 
0 -85A8E 
0 .3096E 
0 .3288E 
C.26C3E 
0 .A0Y5E 
0 . A397E 



-04-0.90 
-03-0. 13 
-03-0 . 16 
-03-0 . 17 
-03-0 . 15 
-03-0 . 12 
-04-C.70 
■04 0.20 
-04 0.95 
•03 0.20 
-03 0.32 
■03 0.41 
-03 0.A8 
-03 0.47 
-03 C.A4 
-03 0.38 
-0 3 C 
•03 C, 
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-03 0.16 
-04 0.90 
-04 O.AO 
-05 0.50 
-04-0.20 
-04-0.39 
-04-0.50 



A 

A 



THIS IS N ARRAY 



0, 

0 . 

0 , 

0 



6COOE 
1 ICOE 
1390F 
. 1530E 
0. 1400^ 
0.7C00E 
-0 . 5000E 
-0. 1750E 
-0. 3C00E 
-0. AlOOE 
-0. AA20E 
-0 .A350E 
-0. 3740 E 
-0.2800E 
-0. 15Q0E 
-0. A500E 
0. 4200 F 
0. HOOF 
0. lACOE 
0. 1320E 
C. 12CCF 
0. 1020E 
0.7500E 
C. A5C0F 
C.900C^= 
POOO C 
P04C C 



-04 0 
-03 0 
-03 0 
-03 0 
-03 0 
04 0 
-04-0 
-C3-0 
-03-0 
-03-0 
-03-0 
-03-0 
-C3-0 
-03-0 
-03-0 
-03-0 
-04 0 
-03 0 
-03 0 
-03 0 
-03 0 
-03 0 
-04 0 
-04 0 
-05 0 
1 5 1. 1 
1600 



.6000E- 
. 1036E- 
. 1278F- 
.1349E- 
. 1197E- 
.5728E- 
.4618E- 
. 1536E- 
. 2598E- 
.3528E- 
. 3800E- 
.3785E- 
.33C0E- 
.2520E- 
.1373E- 
.4373E- 
.3640E- 
.972SE- 
. 1247E- 
.1 193E- 
. 1060E- 
.8750P- 
.6355E- 
.3864E- 
. 1053E- 



04 

03 

03 

03 

03 

04 
■04^ 

03 
03- 
OS- 
OS- 
OS- 

05 
OS- 
OS- 

04 
04 
04 
03 

03 
■03 

04 
■04 
■04 
04 



0.6000E- 
0 .9786E- 
C.1176E- 
0 . 1 1P6E- 
0 . 1012F- 
0 .4572E- 
■0 .4272E- 
-0 . 1342E- 
-0 .2238E- 
-0 .30C8E- 
-0 .3294E- 
•0 .3272E- 
-0.29C0E- 
-0 .2266E- 
-0 . 1257E- 
-0.4257E- 
0 .3 132E- 
0 .8572E- 
0 . 11 C9E- 
0 . 1077E- 
0 .9329E- 
0 .7432E- 
0 .531 SE- 
O.SPBE^- 
0 . 119 lE- 



04 0 
04 0 
03 0 
03 0 

03 0 

04 0 
04-0 
03-0 
03-0 
03-0 
03-0 
03-0 
03-0 
03-0 

03- 0 

04- 0 
04 0 
04 0 
03 0 

03 0 

04 0 
04 0 
04 0 
04 0 
04 0 



.60C0E-04 
.9280E-04 
. 1087E-03 
. 10A2E-03 
. 8497E-04 
.3561E-0A 
.3 96 8E-04- 
. 1 172E-03- 
. 19 13E-03- 
.2552E-03- 
.2 824F-03- 
.2823E-03- 
.2550E-03- 
. 2C43E-03- 
. ) 156E-03- 
. A 156E-04- 
.2687E-04 
.7561E-04 
.9873E-04 
.9761E-04 
.8217E-04 
.6279E-04 
.4405E-04 
.278GF-04 
. 13 13F-04 



0.60C0E-04 0.60 
0.8G9CE-04 0.85 
0.1019E-03 0.95 
0.93C8E-OA C.82 
0.7249E-04 C.60 
0.2780E-0A C.20 
0 .3734E-04-0.35 
0 . 1 CAl E-03-C .91 
0. 1667E-03-0. 14 
0 .22C1F-03-C . 18 
0.2A62E-03-0.21 
0 . 2A76E-03-0 . 21 
0 .2280E-03-0 .20 
0. 1872E-03-0 . 17 
0 . 1C78E-03-0 . 10 
0 .4C78E-04-0 .AO 
0.23A3E-04 0.20 
0.6780E-04 C.60 
0.8936E-04 0.80 
0.8980E-04 C.82 
0.7358E-04 C .65 
0.5390=-04 C.A5 
0.37C2E-04 0.30 
0.2390E-04 C.20 
0.1AC6E-04 C.15 
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P052 


0 . 1600 


PC?.5 


0.2860 


P0C2 


0.0803 


P033 


0.2860 


PCOA 


C . 0A04 


POC5 


G .0052 


P006 


0 . 1660 


POlO 


0 . 2379 


PCI 1 


0.2383 


P012 


0 


P012 


0. C250 


P015 


C.5607 


PO 16 


0 . 1340 


P02G 


0 .C920 


P021 


0.0097 


P022 


0.3576 


POf)3 


C.C905 


P027 


0 .1750 


P055 


0 . 1497 


P05Q 


C.0337 


P05 1 


0.4380 


P056 


0.C092 


P057 


C .9540 


P02 6 


0.1581 


PG5^i 


C.2860 


PC26 


0.C8O3 


pool 


0 . 1750 


P03C 


C.C404 


P04 5 


0 .0920 


P031 


0.0052 


P023 


C.C905 


P03^< 


C . 1660 


P003 


0.2860 


P0A6 


0.4717 


P032 


0.2379 


P035 


0.2383 


P047 


0 . 1497 


P036 


0.0097 


P037 


0 .0250 


PO^ 1 


0 . 1340 


P043 


0.5607 


PO^^ 


0.3576 


P0C7 


0.0337 


POIA 


0.0092 


POA2 


0. 1511 


P017 


0.1581 


AOOl 




OCCl 




A002 




A003 




0003 




A201 




AC06 




A005 




D005 




AGIO 




AOIA 




A016 




A013 




0013 




A017 




00 17 




A022 




A203 




A112 




A060 




A053 




0053 




A 20 6 




A20C 




AOOl 




DO 51 





A057 

DC 5 7 

A202 

A03S 

D035 

A026 

A027 

D027 

A044 

AOS/i 

A022 

C032 

A030 

A036 

AOA 1. 

DOA 1 

AOA3 

C0A3 

A0A2 

A1C4 

A052 

A007 

DOC7 

A205 

A204 

AC65 

D06 5 

A214 

A06). 

D06 1 

A2U 

A 06 € 

A056 

A024 

A2U 

A213 

A207 

A064 

A21C 

A212 



If -if- 



COMPUTER PROGRAM VB 



PHASE I . INTERACTION FORCE AND MOMENT ARE APPLIED TO THE 
* LEADING SHIPISHIP MODEL NO 1) 

-JOB GAL1225 
-FORTRAN LS,GC 

REAL NMAX, NDI MtNSCAL 
REAL NBAR,N,MBAR1 
DIMENSICN Y(30,30) ,N(30,30) 

DIMENSION POT AD (70) ,P0TVAL(70) , AMPAD(70) ,AMFVAL(70) 
INTEGER POTAD,AMPAC 
DATA IFLAG/C/, IFLAGl/0/ 

DATA N0/A8/ 

C STORAGE OF NONDI MENS I CNAL ARRAY FOR NAND Y 
DO 10 J=^l,25 

IC REA0(5, 100) (Y( J, I ) ,1 = 1,6) 

IOC FORMAT ( 6E I 1.4) 

DC 12 J=l,25 

12 REAC(5, 100) (N{ J, I ), 1 = 1,6) 

C MANUAL SWITCH *2' UP FOR BYPASS SETPCTAND SCAN 

1 F(TEST { 2 ) .GT .0) CALL RESET ( 1000 ) ;CALL CCMFUTE;GO TO 1 
X \sPITE{6 ,200) 

2CC FORMAT!' ',‘THIS IS Y(J,I)',//) 



DC llJ=l,25 

X WRITE (6,201) ( Y ( J , I ) , 1=1 , 6) 

201 FORMAT! • ' , 1 OX , 6 ( E 1 1 . 4 , 2X ) , / / ) 

II CONTINUE 
X WRITE(6,300) 

300 FORMAT! • ','THIS IS N(J,I)',//) 

DO 13 J=l,25 

X WRJTE!6,301) ( N ( J , I ) , I =1 , 6 ) 

301 FORMAT!' ' , lOX , 6 ! E 1 1 . 4 , 2X ) , // ) 

13 CCNTINUF 

C SECTION CHECKING POT AND AMP VALUES FCR SHIP HI AND il2 
READ!5, 1005) (POTAD! I ) ,POTVAL! I ) , I = 1,N0) 

1005 FORMAT ( A4 , FIO. 4 ) 

DC 10 5 1 = 1, NO 

105 CALL SETPGT!POTAO! I ) , POTVAL ( I ) ) 

X WRITE(6 ,1015) ! POTAD! I ), POTVAL! I ), 1 = 1, NO ) 

1015 FORMAT! A4, FIO. 4) 

REA0!5,1025) ! AMPAD ! I ) , I = 1 , 66 ) 

1025 FORMAT !A4) 

CALL RESET! 1000) 

DC 125 1=1,66 

125 CALL SC AN! AMPAD! I ) , AMPVAL ( I ) ) 

X WRITE!6,103) ! AM P AD ( I ) , AMPVAL ! I ) , I = 1 , 66 ) 

103 FORMAT! A4, FIO. 6) 

DO 20 1=1, NO 

20 CALL SC AN ! POTAD! I ) , POTVAL ! I ) ) 

X WRITE !6 ,2005 ) ! POT AD ! I ) , POTV AL ! I ) , I = 1 , NO ) 

2005 FORMAT! A4, FIO. 6) 

C WHEN READY 

30 OUTPUT! 101 ) ' READY' 

C transformation of coordinates FOR SHIP /fl AND H2 

115 CALL COMPUTE 

111 CALL ADK(0,U1BAR, 1 ,V1BAR,2 ,PSIBAR,4,U2BAR,5,V2BAR,6, PS 
6R, 3, ALP BAR) 

C RESCALING OF PSIBAR ANC PS2BAR 



PSIMAX=0.26 
PS2MAX=0.26 
PSI=PSI MAX*PS I3AR 
A = SIN( PSD 
B = COS(PSI ) 



XC1D0T=U1 0AR*B-V1BAR*A 
Y01D0T=UIBAR*A-!-V1BAR*B 
PS2=PS2MAX*PS2BAR 
A2 = SIN( PS2 ) 

B2=C0S! PS2 ) 

X C 2 D 0 T = U 2 :: A * B 2 - V 2 B A R ^ A 2 
YC2D0T = L2BAr^-iA2 + V2BAR*B2 
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C INTFRPOLAT ICM IN 2-OlM. ARRAY. PASS DETBAR AND ALPBAR 
C RERERBEH BETDAR=L.O CCRRESPONDS BETA=100.0 FEET*** 

C ALPBAR=1.0 CCRRESPONDS AL PH !5 0 . 0 FEET 
C ALSC BETBAR AND ALPBAR ARE PAST IN COMPUTER UNITS 
X V.RITE(t) ,401 ) 

401 FCRMATI' ‘,10X,'THIS IS AL P B AR- 6 ET 6 AR ' , / / ) 

X WRITE (5, 402) AL PBAR , BETBAR 

402 FORMAT! ' ' , 10X,2F1 I .6) 

C RESCALING OF ALPHAB AR=ALPRAR AND BET ABAR=BETBAR 
ALPHA=550. C*ALPBAR 
X WRITE (6, 403) 

402 FORMAT!' ',10X,‘THIS IS ALPHA',//) 

X WRITE{6 ,501 ) ALPHA 

BETA=100.0*BETBAR 
X WRITE(6,404) 

404 FORMAT! ‘ ',1.0X,‘THIS IS BETA',//) 

X WRITE(6,501)BETA 

IE(8ETA.LT.50.0) OUTPUT(lOl) B E T A ; NSC AL = 0 . C ; YSC AL = 0 . 0 ; 
IF(BETA .GT. I 00.0) OUTPUT (101 ) BET A ; N S CAL= 0 . 0 ; YSC AL = 0 . 0 ; G 
IF( ALPHA.GT.550.0)OUTPUT! 10 1 ) AL PHA ? NSCAL= C , C ; YSCAL^O .0 
I F! ALPHA. LT.-550 eOOUTPUT ! 1 0 1 ) AL P FA ; NSC AL = 0 .0 ;YSCAL = 0. 
C SEARCH FOR INDET IF I CAT I ON OF ARRAY ELEMENT 

C FIND CORRESPONDED Y AND N -LINEAR INTERPOLATION IN 2-DIMEN 
I ei'-^BETA/lO.O 
I E--=IBi-4 
BETA1=I Bl*10 
WTB=( BETA-BETAl ) /lO.O 
IF! ALPHA.lt. 0.0) GO TO 1 
IFIALPHA.GT.O.O) GO TO 2 
C ALPHA =0 
IA=13 

wta=o.o 

GO TO 3 

C ALPHA LESS THAN ZERO 

1 I Al=-ALPHA/50 .0 
ALPHA1=-IA1*50 

WTA=-( ALPHA- ALPHA 1 )/50.0 
I A = 13-I A1 
GO TO 3 

C ALPHA GREATER THAN ZERO 

2 I Al=ALPhA/50.0 
IA=13-M A1 
ALPhAl= IA1*50 

WTA=( ALPHA-ALPHAl ) /50.0 

3 IE2=IB+1 
IFUB.FC.6) IB2 = 6 
I A2 = IA-i- 1 

IE!IA.EC.25) IA2=25 
X WRITE!6,101) IBl , I B, IB2, I A1 

101 FORMAT!' ‘ , lOX , ' I B 1= ' , I 3 , 5X , ' I B = ' , I 3 , 5X , ' I D2 = ' , 1 3 , 5X , ' 

X WRITE(6,102) IA,IA2 

102 FORMAT!* ' , 1 OX , ' I A = ' , 1 3 , 5 X , ' I A2 = ' , I 3 , / / ) 

X WRITE (6, 202) BETA I , WTB , ALPHA 1 

202 FORMAT!' • , • B ET A 1 ' , E 14 . 8 , 2X , ‘ WT B ' , E 1 4 . 8 , 2 X , ' ALPHA I ' , E 1 

X WRITE!6,203) WTA 

203 FORMAT! « ‘ » WT A= • , E 14 . 8 , / / ) 

IF! ALPHA. LT. 0.0) IA2=IA-1 
IF(IA.EQ.l) IA2=1 
IF(WTA.EQ.O.O)IA2=IA 
JF(WTB. EQ.0.0) IB2=IB 
N'EAR1 = N( lA, IB2) 

X WPITE(6,602) 

602 FORMAT!* ','THIS IS NBARI',//) 

X WRITE(6,502) N3AR1,IA,IB2 

N'BAR=N! I A, IB ) 

X WRITE!6,603) 

602 FORMAT!' ','THIS IS N8AR',//) 

X WRITE !6, 502 5 NBAR, lA, IB 

RARM= ( NBARl-NBAR ) *WTG + NBAR 
X WRITf- ( 6 ,604 ) 

6C4 FORMAT (« ','THIS IS BARNl',//) 

X WRITE!6,501) BARNl 
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NBARUN ( IA2, IB2J 
X WRITE(6,605) 

6C5 FCPKAK' ','THIS IS NBARl',//) 

X WRITE(6,I>02t NBAR 1, IA2t IB2 

NEAR = N( IA2 , IB ) 

X WRITE(6,606) 

606 FCPRAT{« S'THIS IS NBAR',//) 

X WPITE(6,602) NBAR,IA2,IB 

BAR(\2= (NBAR 1-N BAR) *WTB + NBAR 
X V>/RITL-(6,607) 

607 FORMAT!' ','THIS IS BARN2',//) 

X KRITE(6,501) BARN2 

NBAR=( BARN2-BARNI }«WTA+BARN1 
X RR1TE(6,500) 

500 FORMAT!' ','THIS IS FINAL NBAR',//) 

X RITE (6, 501) NBAR 

501 FORMAT! • ' , lOX , E 1 A . 8 , / / ) 

YBARI=Y (1 A, IB2 ) 

X KRnE(6,608) 

608 FORMAT!* *,'THIS IS YBARL',//) 

X WRITE(6,502) YBAR1,1A,IB2 

YBAR=Y( lA, IB) 

X WRITE (0,605) 

609 FORMAT!* ','THIS IS YBAR',//) 

X WRITE( 6,502) YBAR, lA, IB 

BARY1.= ( YBARl-YBAR) ^WTB-s-YBAR 
X WRITE (6, 610 

61C FORMAT!' ',‘THIS IS BARYl',//) 

X WRITE (6,601 ) BARYl 

YBAR1 = Y ( IA2, I B2 ) 

X WRITP(6,61.l) 

611 FORMAT!' *,'THIS IS YBAHl',//) 

X WRITE(6,502) YB AR I , 1 A2 , 1 B2 

YBAR=Y( IA2, IB ) 

X WRITE!6,612) 

612 FORMAT!' ','THIS IS YBAR',//) 

X WRITE (6, 502) YBAR, I A2, IB 

BAPY2 = ( YBARl-YBAR) ^-WTE+YBAR 
X WRITE (6, 613) 

613 FORMAT!' ','THIS IS BARY2',//) 

X WPITE(6,601) BARY2 

Y6AR= ( BARY2-BARY I ) *WTA+BARY1 
X WRITE !6, 600) 

6CC FORMAT!' ','THIS IS FINAL YBAR',//) 

X WRITE(6,60I) YBAR 

601 FORMAT!' ‘ , 1 0 X , E 1 4 . 8 , / / ) 

502 FORMAT!' ' , lOX , E 1 4 . 8 , 5 X , ' RG W = ' , 1 3 , IX , ' C OL = ' , I 3 ) 

0 REMEMBER SIMULATION IS DONE IN OGMPUTER UNITS SO NO NEED T 
CIVIDE BY 100,0 SINCE SUBR. DAC CONVERTS BY ITSELF IN 
C MAGNITUDE SCALING OF N 
NMAX=4, I95E+07 
ND I M=NB AR* ! 9 . 42 8 E+ 1 0 ) 

NSCAL=NDIM/NMAX 
C MAGNITUDE SCALING OF Y 
YMAX=1 . 608E+05 
YOIM.=YBAR*! 17.87E+C7) 

YSCAL=YDI M/YMAX 

804 CALL DAC! 1,X01D0T,2,YQ1D0T,5,Y02D0T,6,X02DCT,7,NSCAL,8 
C MANUAL DIGITAL SWITCH'D' AND 

15 IF!TEST(1) ,GT,0) GC TO 16 
CALL HOLD 

17 IF !TEST ! I ) .LT .0 ) GO TO 17 
CALL COMPUTE 

16 IF!TEST !3) .GT.O) GO TO 111 
CALL RESETllOOO) 

PALSE 
GO TO 115 
END 

-LOAC XR,MAP 
-,DAT A 

* 
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THIS IS Y ARRAY 



* 

:« 

* 

-0 . MOOE- 
-0. 2^.C0E- 
-0. 3200F- 
-0. 38C0E- 
-0. ACCOE- 
-0. 2500E- 
-C. 2 6C0E- 
-C. 12C0E- 
AOCOE- 
27CCE- 
500CE- 
7000E- 
85COE- 
9000E- 
82C0E- 
0. 6250E- 
O.AECOE- 
0. 3CCCE- 
0. 1700E- 
0.6C00E- 
-0. 20C0E- 
-0. 5C00E- 
-0. 6CC0E- 
-0. 5C00F- 
-0. lOOOE- 
* 

* 

* 

* 

* 



0 , 
0 . 
0 , 
0 . 
0. 
c. 
0 , 



03-0 

03-0 

03-0 

03-0 

03-0 

03-0 

03-0 

03- C 
04 0 
C3 0 
03 0 
03 0 
03 0 
03 0 
03 0 
03 0 
03 0 
03 0 

03 0 

04 0 

04- 0 
04-0 
04- C 
04-0 
04-0 



.1263E- 
.2112E- 
.21h2^~ 
.3225E- 
„3337F“ 
.2870F- 
.2079E- 
.93 84E- 
.5507E- 
.25C8E- 
.45C7F- 
.6225E- 
.75CQE- 
.7 844E- 
.7159E- 
. 55 79F- 
.4089E- 
.2781.E- 
. L673E- 
.6822E- 
.3562E- 
• 3493E- 
.49 04E- 
.4699E' 
.2096E- 



03- 

•03- 

03- 

•03- 

03- 
•03- 
•03- 

04- 
■04 
■03 
03 
•03 
•03 
■03 
•03 
03 
•03 
•03 
03 
•04 
05 
•04- 
04- 
04- 
04- 



0.1160E 
0 . 1897E 
0.2433E 
■0 .2793E 
0.2840F 
0 .2397E 
0 »I689E 
0 c667 LE 
0 .6637E 
0.2364E 
0.4137E 
0 .5643E 
0 .675CE 
0 .6977E 
0 .6378E 
0 .5075E 
0 .3781E 
0.2616E 
0 . 1552E 
0 .7438E 
0 . 8767E 
C .2363E 
■0 .4082E 
0 .4473E 
0 .2918E 



•03-0 
-03-0 
-03-0 
-03-0 
-03-0 
-03-0 
-03-0 
-04-0 
-04 0 
-03 0 
■03 0 
■03 0 
-03 0 
■03 0 
-03 0 
•03 0 
-03 0 
■03 0 
-03 0 
•04 0 
-05 0 
-04-0 
-04-0 
-04-0 
-04-0 



. 1058E-03 
. 168 lF-03 
.2 1C4E-03 
.2362E-03 
.2342E-03 
. 1.925E-03 
. 1299F-03 
. 39 59P-04 
.7767E-04 
.2221E-03 
. 3767E-03 
.5062E-03 
.60C0E-03 
.6 llOE-03 
.5597E-03 
.4572E-03 
.3473E-03 
.2452E-03 
. L632E-G3 
.8055E-04 
.2nCF-04 
.1233E-0A 
.3260E-04 
. 4 24 7 E- 04 
.374CE-04 



-0 



0 .9753E 
0 . 15C8E 
0 . I 8 4 I F 
0 .20 16E 
0.1.945F 
1547E 
0 .9863E 
0. 1789E 
0 .8671E 
0 .21C5E 
0 .3471E 
0 .4596E 
0 .54C0E 
0 .54 16E 
0 .4973E 
0 .4 169E 
0 .3226E 
0.2321E 
0 . 16 1 5F 
0 .8548E 
0 .3C96E 
0.3288E 
0.2603E 
0.4C65F 
0 -4397E 



04-0.90 
03-0 . 13 
03-0 . 16 
03-0 . 17 
03-0 . 15 

03- 0 . 12 

04- C .70 
04 0.20 
04 0.95 
03 0.20 
03 0.32 
03 0.41 
03 0.48 
03 0.47 
03 C.44 
03 0.38 
03 C.30 

0 .22 

03 0.16 

04 C.90 

04 0.40 

05 C.50 
04-C. 20 
04-0.39 
04-C. 50 



-0 



THIS IS N ARRAY 



0 . 

- 0 , 

- 0 , 

- 0 , 



0. 6CC0E 
0. 1 ICOE 
0 . 1390E 
0. 1530E 
0. 1400'" 
7CC0E 
5C00E 
1750E 
3C00E 
-0. 4 ICOE 
-0.4420E 
-0. 4350E 
-0.3740E 
-0. 2800 E 
-0. iSCOE 
-0.45C0E 
0. 4200E 
C. 1 ICOE 
0. 1400E 
0. 1320E 
0. I2C0E 
, 1020E 
75G0E 
4 5 COE 
9C00E 
0 . 



0 , 
0 . 
0 . 
0 



POCC 

P040 

P052 

P025 

P002 

P033 

P004 

P005 

PCC6 

POIC 



C . 
C. 
0 . 
0 . 
0 . 
C . 
0 . 
C. 
0 . 



-04 C 
-03 0 
-03 0 
-03 0 
-03 0 
-04 0 
-04-0 
-03-0 
-03-0 
-03-0 
-03-0 
-C3-0 
-03-0 
-03-0 
-03-0 
-03-0 
-04 0 
-03 0 
-03 0 
-03 0 
-03 0 
-03 0 
-04 0 
-04 0 
-05 0 
1511 
1600 
1600 
2860 
0803 
2860 
0404 
0052 
1660 
2379 



.6000E- 

. 1036E- 

. L278E- 

. 1349E- 

. 1197E- 

.572SF- 

.4618E- 

. 1536E- 

.2598E- 

.3528F- 

.3880E- 

.3785P- 

.33C0E- 

.2520E- 

. 1373E- 

.4373E- 

, 3640E- 

.9728E- 

.1247E- 

. 1193E- 

. I060E- 

.8750E- 

.6355E- 

.3864E- 

. 1053E- 



04 

03 

03 

03 

03 

04 
04- 
03- 
03- 
03- 
03- 
03- 
03- 
03- 

03- 

04- 
04 
04 
03 
03 

03 

04 
04 
04 
04 



0 .6000E- 
0 .9786E- 
0.1176E- 
II86E- 
10 12E- 
4572E- 
42 72E- 
1342E- 
0 .223 8E- 
0 .30C8E- 
0 .3294F- 
0 .3272E- 



^6000E-04 
9280E-04 
. 1087E-03 
. 1042f^-03 
I. 8497E-04 
.3561E-04 
I.3968E-04- 
. 1 172F-03- 



29C0F- 
.2266E- 
. 1257E- 
.4257E- 
0 .3132E- 
0 .8572E- 
.1109E- 
.1077E- 
.9329E- 
.7432E- 
. 5315E- 
.3286E- 
.1191E- 



04 0 
04 0 
03 0 
■03 0 

03 0 

04 0 
04-0 
03-0 
03-0. 19 13E-03- 
03-0.2552E-03- 
03-0.2824E-03- 
03-0.2S23E-03- 
03-0.255CE-03- 
03-0 .2043F-03- 

03 - C. 1 156E-03- 

04- 0.4 156E-04- 
04 0.2667E-04 
■04 0.7561E-04 
03 0.9873F-04 

03 0.9761E-04 

04 0.8217E-04 
04 0.6279E-04 
04 0.44C5E-04 
04 0.278CE-04 
■04 0.13 13E-04 



0.6CC0E-04 0.60 
0.8890'^-04 0.85 
0.1C19E-03 C.95 
0.93C8F-04 0.82 
0.7249E-04 0.60 
0.2780F-04 0.20 
0 .3734E-04-0.35 
0 . 1C41E-03-C .91 
0 . 1667E-03-0 . 14 
C . 22CIF-03-C . 18 
0 .2462E-03-0 .21 
0 .2476E-03-0 .21 
0.228CE-03-C.20 
0 . 1872E-05-0. 17 
0 . 1C78E-03-0 . 10 
0.4C78E-04-0.40 
0.2343E-04 0.20 
0.6780F-04 0.60 
0.8936E-04 0.80 
0.8980E-04 0.82 
0.7358E-04 C.65 
0.5390E-04 0.45 
0.37025-04 0.30 
0.239CE-04 0.20 
0.14C6E-04 0.15 
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POD 


0.2383 


P012 


C.4717 


P012 


0.0250 


P015 


0.5607 


POD 


C. 1340 


P020 


0.C920 


P02 1 


0.C097 


PC22 


0.3576 


P053 


C .C905 


P027 


0.1750 


P055 


C. 1497 


P050 


0.0337 


P051 


0.4380 


P05^ 


0 .C092 


P057 


0.9540 


P02A 


C. 1581 


P05^ 


C .2860 


P026 


0.0803 


POOL 


0 .1750 


P03C 


C . 0404 


P04f3 


0.0920 


P031 


0.0052 


P023 


C.C905 


P034 


0 . 1660 


POOS 


0.2860 


P046 


0.4717 


P032 


0.2379 


P035 


0.2383 


P047 


0. 1497 


P03 6 


0.0097 


P03 7 


C.C250 


P04 L 


0 . 1340 


P04 3 


0.5607 


P044 


C .3576 


P0C7 


0.0337 


POl 4 


0 . C092 


P04 2 


0.1511 


PC 17 


0.1501 


AOOl 




0001 




A002 




A003 




DCC3 




A201 




A006 




A005 




0005 




AOlO 




A014 




AO 16 




A013 




DC13 




AO 17 




0017 




AC22 




A203 




A112 




A06C 




A053 




0053 




A2C6 




A200 




A051 




0051 




A057 




0057 




A202 




A035 




0035 




A026 




A027 




002 7 





A03A 
A033 
DC33 
A03C 
A036 
AO^ 1 
DOA 1 
AOA3 
DOA2 
A0A2 
AlOA 
A052 
A007 
COC7 
A205 
A2CA 
A065 
C065 
A21A 
A06 I 
D061 
A21 1 
A066 
A056 
A02A 
A216 
A213 
A2C7 
A06A 
A210 
A212 



COMPUTER PROGRAM VI 



* PHASE II .INTERACTION FORCE AND MOMENT ARE APPLIED ON 

* BOTH SHI PS 
*^JOB GAL1225 
-iFCRTRAiN LS,GC 

REAL NMAX, NDIM, NSCAL.NBSCAL 
REAL NBAR,N,NBAR1 
DIMENSION Y(30,30) ,N(30,30) 

DIMENSICN POT AO (70) ,P0TVAL(70)r AMPAD(70 ) » AMPVAL(70) 
INTEGER POTAD,AMPAD 
DATA NO/48/ 

C STORAGE OF NOND I HENS I CNAL ARRAY FOR NAND Y 
DC 10 J = l, 25 

10 READ(5,100) (Y( J, I ) , 1^1,6) 

ICC rCRMAT( 6E1 1 .4) 

DC 12 J=L,25 

12 READ(5tl00) ( N ( J , I ) , I - 1 , 6 ) 

X LRITE(6,200) 

200 FCRMAT(* «,'THIS IS Y(Jtl)',//) 

DC 11J=1,25 

X LRITE (6,201) ( Y ( J , I ) , I = 1 , 6 ) 

201 FCRMAT(‘ ‘ , lOX , 6 ( E 1 1 . 4 , 2X ) , / / ) 

11 CONTINUE 

X RRITE(6,300) 

300 FCRMATI' ','THIS IS N(J,I)',//) 

DC 13 J=l,25 

X LRITE (6,301) ( N ( J , I) , I = 1 , 6 ) 

301 FCRMAT(» ' , 1 OX , 6 ( E 1 1 . 4 , 2X ) , // ) 

13 CONTINUE 

-C SECTION CHECKING POT AND AMP VALUES FOR SHIP HI AND H2 
RE AO (5, 10 05) (POT ADI I ) ,POTVAL( I ) , I = 1,N0) 

1005 FORMAT! A4 ,F10. 4) 

DC 10 5 1 = 1, NO 

105 CALL SETPGTI POTAD( I ) , POTVAL ( 1 ) ) 

X WRITE ( 6 , 1015) ( POTADI I ) , POTVALI I ) , 1 = 1 ,N0) 

1015 FCRMAT(A4,F10.4) 

REA0I5, 1025) ( AM P A C ( I ) , I = 1 , 66 > 

1025 FCRMAT(A4) 

CALL RESET (1000) 

DC 125 1=1,66 

125 CALL SCAN! AMPAD( I ) , AMPVAL ( I n 
X WRITE(6,103) (AMPADI I ) ,AMPVAL( I ) , 1=1,66) 

103 FCRMATI A4 , F10.6) 

DC 20 1=1, NO 

20 CALL SC AN { POTADI I ) , POTVAL { I ) ) 

X WRITE (6 ,2005) (POTADI 1 ) , POTVALI I ) ,1 = 1, NO) 

2005 FORMAT! A4,F10. 6) 

C WHEN READY 

30 OUTPUT! 101) 'READY' 

C TRANSFORMATION OF COORDINATES FOR SHIP #1 AND H2 
115 CALL COMPUTE 
I FLAG=0 

111 CALL ADK(0,U1BAR, 1,V1BAR,2,PSIBAR,4,U2BAR,5,V2BAR,6,PS 
6R,3,ALPBAR) 

C RESCALING OF PSIBAR AND PS2BAR 
PSIMAX=0.26 
PS2MAX=C. 26 
PSI=PSIMAX*PSIBAR 
A=SIN( PSI ) 

B=COS( PSI ) 

XC1D0T=U1EAR*B-V1BAR*A 
YC1D0T = U1BAR*A+V1 BAR-^B 
PS2=PS2MAX*PS2BAR 
A2=SIN(PS2) 

B2=C0S( PS2) 

X02D0T=U2CAF-"-B2-V2BAR*A2 
Y02 00T = U2BAR'^A2+V2eAR*B2 

Q T T PD DPI A J T I p_nTM AOq.AY ^A^^ DCTPAD A I O P 

C REMFM EER BETBAR= 1 ‘o'cOP.RL SPONDS BET 166.0 f'EET^*^’ 
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C ALFBAR = ]..0 CORRF.SPONDS ALPHA=--BBO . 0 FRET 
C ALSO BFTBAR AND ALPBAR ARE PAST IN CCRPUTER UMTS 
X URITE(6,ACl) 

AOl FCRMATI* MlOXj'THIS IS AL PB AR- B ET EAR ' , / / ) 

X WPITE(6,A02> ALPBAR, BFTBAR 

A02 FCPMAT( » • , 10X,2F11.6) 

C RESCALING OF ALPHABAR= ALPBAR AND BE T A E AR=BET B A R 
ALPHA=550 . C*ALPBAR 
55S CONTINUE 
X NRITEI6,A03) 

A03 FCRNATf' ‘,lOX,'ThIS IS ALPHA',//) 

X WRITE (6 ,501) ALPHA 

BETA= 100. O^BE TEAR 
X WRITE(6,40A) 

A04 FGPMAT(« ',10X,'THIS IS BETA',//) 

X WRITE! 6,501)BETA 

IF(BETA.LT.50.0) OUTPUTUOL) BE TA ; NSCAL = 0 . 0 ; YSCAL = 0 . Cl 
I F(BETA .GT .100.0 ) OUTPUT ( 10 I ) BET A ; NSCAL=0 . 0 ; VSCAL = 0 .0 ;G 
I F( ALP HA. GT. 5 50.0 I OUT PUT I 10 L ) AL PH A : NSC A L = C . C ; V SC A L = C . 0 
IF( ALPHA. LT .-550 .0 lOUTPUT ( I 0 L > AL FFA ; NSC AL = 0 .0 ;YSCAL = 0. 
C SEARCH FOR IND ET 1 F I C AT I CN OF ARRAY ELEMENT 

C FIND CORRESPONDED Y AND N -LINEAR INTERPOLATION IN 2-DIMEN 
I ei=BETA/lO.O 
IB= IB L-4 
BETAi=I Bl=^10 
WTB=(BFTA-BETAl) /lO.O 
IF( ALPHA. LT. 0.0 ) GO TO I 
IF( ALPHA. GT. 0.0) GO TO 2 
C ALPHA =0 
1 A=13 
WTA=0. 0 
GO TO 3 

C ALPHA LESS THAN ZERO 

1 I Al=-ALPHA/50.0 
ALPHAL=“I Al*50 

WTA=-I ALPHA- ALPHA I )/50.0 
IA=13-I Al 
GO TO 3 

C ALPHA GREATER THAN ZERO 

2 I Al=ALPHA/50.0 
IA=13+! Al 
ALPHAl- I Al^=50 

WTA=( ALPHA-ALPHAl ) /50. 0 

3 IE2=IB+1 
IF(IB.EC.6) 1B2=6 
IA2^IA-M 

IF { lA . EG. 25) I A2 = 25 
X WRITE(6,101) IBl , IB, IB2, lAl 

101 FORMAT!' ' , lOX , ' I B 1= ' , I 3 , 5X , ‘ I B = ' , I 3 , 5X , ' I P 2= ' , I 3 , 5X , ‘ 

X WRITE (6 ,102) lA, IA2 

102 FORMAT!' ' , lOX , • I A= ' , I 3 , 5 X , ' I A2 = ' , I 3 , // ) 

X KRITE(6,202) BETA 1 , WTB , ALPHA 1 

20 2 FORMAT ( ' ' , ' BETAl' , E 1 4 . 8 , 2X , ' NT B ' , E 14 . 8 , 2 X , ' ALPHA 1 ' ,E1 

X WRITE(6,203) WTA 

203 FORMAT! ‘ ' , ' WT A= ' , E 14 . 8 , // ) 

IF! ALPHA.LT.0.0) IA2=IA-1 
IF!IA.EG.l) IA2=1 
IF!WTA.EQ.0.0)IA2=IA 
IF(WTB. EQ.0.0 ) IB2 = IB 
NBARl = N( lA , IB2) 

X WRITE! 6, 60 2) 

602 FORMAT!* ','THIS IS NBARl',//) 

X WRITE!6,502) NBARl, IA,IB2 

NBAR=N! I A, IB) 

X WPITE!6,603) 

603 FORMAT!’ ','THIS IS NBAR',//) 

X KPITE(6,502) NBAR, lA, IB 

BARM = ! NBARl-NBAR) *WT E+NBAR 
X WRITE (6, 604) 

604 FORMAT!' ','THIS IS BARNl',//) 

X WRTT=( 67501 ) OARNl 

I\BAR1: = N ( I A2, IB2) 
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X loPlTE (6 ,605) 

605 FCRf-'AT(« ',‘THIS IS NBARl',//) 

X WRI1F(6,502) NBAR I , IA2, IB2 

N6AR = N( IA2 , IB 1 
X l^'R ITl- ( 6,606) 

606 FCRMATt* *,'THIS IS NBAR',//) 

X V^RITE(6,502) K'BAR,!A2,IB 

BARK2= { N'BAR L-NBAR ) -M-iTB+NBAR 
X WRITE (6,607) 

6C7 FORMAT! • ','THIS IS EARN?.',//) 

X WRITE (6, 50 1) ('ARN2 

NBAR=( BARN2-BARNI ) *WTA+BARN1 
X WR1TE(6,500) 

50C FORMAT! • ‘,'THlS IS FINAL NBAR',//) 

X WRITF(6,50I) NBAR 

5CL FORMAT!' ' , lOX , E 1 A . 8 , / / ) 

YEAPL=Y(IA,IB2) 

X WRITE(6,608) 

608 FORMAT! • ','THIS IS YBARl',//) 

X WRITE(6,502) YBARl, lA, 182 

YBAR=Y( lA, IB) 

X WRITE!6,605) 

609 FORMAT!' ','THIS IS YBAR',//) 

X WR1TE!6,502) YBAR, lA, IB 

BARYl= ( YBARl-YBAR) ''!‘WTfi+YBAR 
X WRITE(6,61C) 

61C FORMAT!' ',‘THIS IS BARYl',//) 

X WRITE!6,601)BARY1 

YEARl=Y!IA2,IB2) 

X WRITE (6,61 1) 

611 FORMAT!' ‘,'THIS IS YBARl',//) 

X WRITE( 6,502) YB AR 1 , I A 2 , I B2 

YBAR=Y! IA2, IB ) 

X WRITE(6,612) 

612 FORMAT! • ','THIS IS YBAR',//) 

X WRJTE(6,502) YBAR,IA2,IB 

B A R Y 2 = ! Y 8 A R 1 - Y B A R ) '■' W T B -F Y B A R 
X WRITE!6,613) 

613 FORMAT !' ','THIS IS BARY2',//) 

X WRITE(6,601) BARY2 

YBAR=! BARY2-BARY 1 ) *WTA+BARY 1 
X WRITE (6, 600) 

6CC FORMAT!' ','THIS IS FINAL YBAR',//) 

X WRITE!6,601) YBAR 

601 FORMAT!' ' , 1 0 X , E 1 4 . 8 , / / ) 

502 FORMAT!' ' , lOX , E 1 A . 8 . 5 X , ' RO W = ' , I 3 , IX , ' 0 OL = ‘ , I 3 ) 
REMEMBER SIMULATION IS DONE IN COMPUTER UNITS SC NO NEED T 
CIVIDE BY 100.0 SINCE SUBR. DAO CONVERTS BY ITSELF IN 
WRITE ! 6 , 55555 ) IFLAG 
^5555 FORMAT!' ' , 5X , ' I FL AG= ' , 1 1 , / / ) 

C MAGNITUDE SCAI. ING OF N 
NMAX=4. 195E + 07 
NCIM.=NBAR* (9 .428E+10) 

IF(IFLAG.EQ.O) NSCAL=NDIM/NMAX 
IF ( IFLAG .EQ. I ) N B S C AL = ND I M/ NMAX 
C MAGNITUDE SCALING OF Y 
YMAX=1. 6C8E+05 
YDIM=YBAR« ! 17 .87E+07 ) 

1 F( IFLAG.EQ.O) Y S C AL= YD I M/ Y MAX 
IF! IFLAG. EQ.l) Y BSC AL = YD I M/ YMAX 
X WRITE (6, 901) NDIM,YDIM 

901 FORMAT! • ‘ ' NDI M = • , E 1 4 . 6 , 5X , ' YD I M= ' , E 1 4 . 6 , / / ) 

X WRITF!6,5001 ) XO ID GT , YO IDOT 

5001 FORMAT!' ' , • XOl DOT = ' , E 14 . 6 , ' YO 1 CCT= ' , E 1 4 . 6 , // ) 

X WRITE!6,5000) X02D0T , Y02D0T 

5000 FORMAT!' ' , '.X02D0T= ' , E 14 . 6 , ' Y02D0T= ' , E 14 . 6 , // ) 

IF! IFLAG. EC. 0) I F L AG= 1 ; AL PH A=-A L PF A ; GO TO 555 
X WRITE(6,900) NSC AL , YSC AL ,NBSCAL , YBSCAL 

900 FORMAT (' ' , 'NSCAL=' ,E14.6, 2X, ' YSCAL=' , E 1 4 . 6 , 2 X , ' N B SC AL 

6' YBSCAL = ' , 014.6, // ) 

8C4 CA'. L O AC! L .XOIDOT , 2 , Yf^ I DOT , 5 , YO 2C0T , 6 , XO 200 T , 3 , N SC A L 
6CAL, 8, YBSCAL) 
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C MANUAL digital SWITCH'D' AN'D 'I' 

15 IF (TEST ( I ) .GT .0 ) GC TO 16 
CALL HOLD 

17 IFITESTd) .LT.O) GC TO 17 
CALL COMPUTE 

16 IFCTEST (3) .GT.O) IFLAG=OlGOTO 111 
CALL RESET(IOOO) 

PAUSE 
GC TO 115 
END 

-.LCAC XP,MAP 
DATA 



THIS IS Y ARRAY 



-0. 1 ACC F- 03-0 . 1263E -03-0. 11 60E-03-0. 10 58 E-03-0 .97 53E-04-0. 
-0. 2AC0E-03-0 .21 12F-03-0 . 1897E-0 3-0. 1681E-03-C. 15COE-03-Q. 
-0. 3 2 00F-03-0 .2 762E-03-0 . 2433 E-0 3--0 . 2 1 0 A F-03- 0 . 1 8 A 1 E-03- 0 . 
-0. 3 800 E- 03-0 .3 22 5E-03-0 .2 793E-03-0.23 62E-03-0 .2C16E-03-0 . 
-0 .AOCOE-0 3-0 . 3337E-03-0 .2840E-03-0.2 3A2E-03-0 . 19A5E-0 3-0 . 
-0. 35C0E-03-0 .2 870E-C3-0.2397E-03-0. 19 25E-03-0 . 1547 E -03-0 . 
-C. 26CCE-03-0 .2079E-03-0 c 168 9E-03-0 , 1 299E-03- 0 . 9 8 6 3 E-04- 0 . 
-0. 1 3 00E-03-0 .9384E-04-0 . 667 1 E-04-0 . 39 59 E-04- 0 . 1789E-04 0 . 
0.4CCCE-C4 0.55C7E-04 0.6637E-04 0.7767E-04 0.8671E-04 0. 
0.2700E-03 0.25C8E-03 0.236AE-03 0.2221E-03 C.21C5E-03 0. 
0.5000E-03 0.45C7E^03 0.4137F~03 0.3767E-03 0.3A71F-03 0. 
0.7000E-03 0.6225F-03 0.5643E-03 0.5062E-03 0.4596^-03 0. 
0.8500E-03 0.7500F-03 0.6750E-03 0.600CE-03 0.5AC0E-03 0. 
C.9CCCE-03 0.7844E-03 0.6977F-03 0.6110E-03 0.5A16E--03 0 . 
C.82C0E-03 0.7159E-03 0.6378E-03 0.5597F-03 0.4973E-03 0. 
0.6250E-03 C.5579F-03 0.5075E-03 0.4572E-03 0.4169E-03 0. 
0.A5CCE-03 0.4089F-03 0.3781F-03 0.3473E-03 0.3226E-03 0. 
0.3000E-03 0.2781F-03 0.2616E-03 0.2452F-C3 0.2321E-03 0. 
0.17CCE-03 0.1673E-03 0.1652F-03 0.1632E-03 0.1615E-03 0. 
0.60C0E-C4 0.6822F-04 0.7438E-04 0.8055E-04 0.8548E-04 0. 
-0.2CC0E-04-0 .3562E-05 0.8767E-05 0.2U0E-04 0.3C96E-04 0. 
-C. 5000 E-04-0 .34<^3 E-04-0 .2363E-04-0. 12 33 E-04-0 .3 288E-05 C . 
-C.6CCGE-04-0 .49 04E-04-0 c4082E-04-0.3 26CE-OA-C . 26C3E-04-0 . 
-0.5000 E-04-0 .4699 E-04-0 .447 3 E-04-0. 4 247 F-04-0 .4C65E-04-0 . 
-0. lCCCE-04^0 .2096E-04-0 .2918E-04-0.37A0E-04-0 .A397E-04-C . 
* 

A 

* 

* THIS IS N ARRAY 

C.6CC0E-04 C.6CC0E-O4 0.60COE-04 0.6CCCE-04 0.6000E-04 0. 
O.llCOE-03 0.1036F-03 0.9786E-04 0.9280E-04 0.8890E-04 C. 
0.1290E-03 0.1278E-03 0.1176E-03 0.L087E-03 0.1C19E-03 0. 
0.I-530E-03 0.1349E-03 0.1L86E-03 0. 10A2E-03 C.93C0E-O4 0. 
O.lAOOE-03 0.1197E-03 0.1012E-03 0.8497E-04 0.72A9E-04 0. 
0.7C00E-04 0.5728E-04 0.4572E-04 0.3561E-04 0.2780E-04 0. 
-0. 5C00E-04-0 .4618E-04-0 . 427 2E-04-0 . 3 9 6 8E-04- 0 . 37 34E-04-0 . 
-0. 1 7 5CE-03-0 .1536E-03-0 . 1342E-03-0 . 1 172E-03-0 . lCAlE-03-C . 
-C. 3C00E-C3-0 .2 598E-03-0 .223 8E-03-C. 19 13E-0 3-0 . 16 67F-03-0 . 
-0 . A lOOE-03-0 .3528E-03-0 .30C8F-03-0 .2552E-03-0 .22C1E-03-C. 
-0 . AA2 0E-03-0.3 880E--03-0 .3294E-03-0.2 8 2A E-03-0 .2 4 62 E -03-0 . 
-0. A350E-C3-0 .3 785P-03-0 .3 272E-0 3-0.2 8 23E-03-0 .24 76F-03-0 . 
•^0. 3 740E-03-0 .33C0E-03-0 . 29 CO E-0 3-0 . 2 5 50 E-03-0 .22 80E-0 3-C . 
-0. 2 800 E-03-0. 2 5 20E -0 3 -0.22 66E- 03-0. 2 04 3 E-03-0 . 1872 E-03-0 . 
-0 . 1500 E-03-0 . 1373E-03-0 . 1257E-03-0. 1 156 E-03-0 . 1C 78 E-03-0. 
-C. A 500 E-03-0 .4373 5-04-0 .42 57 E-04-0 .4 156 E-04-0 .AC78 E-04-0 . 
0.4200E-0A 0.3640E-04 0.3132F-04 0.2687E-04 0.22A3E-04 0. 
O.llOOE-03 0.9728E-04 0 .8572E-04 0.7561E-04 0.678CE-0A C. 
O.lACOE-03 0.12476-03 0.UC9E-03 0.9873E-04 0.8926E-04 0. 
C.122CE-03 C.1192'^-03 0.1077'^-0? 0.'^761‘'-04 n.pcpO'^-P^ 0- 
0.12C0E-O3 0.1060E-03 0.9329E-04 C.8217E-04 0.7258E-04 0. 



90 

13 

16 

17 

15 

12 

70 

20 

95 

20 

22 

41 

48 

47 

44 

28 

30 

22 

16 

90 

40 

50 

20 

39 

50 

60 

85 

95 

82 

60 

20 

25 

91 

14 

18 

21 

21 

20 

17 

10 

AO 

20 

60 

80 

P? 

c 5 
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Oc !C20E“03 0 
0.7500E-04 0 
O.AS'COE-04 0 
C.90C0E-05 0 



P040 


0 .1600 


P052 


C . 160C 


POOO 


0.1511 


P025 


0.2860 


PCC2 


C .0803 


PC23 


C .2860 


P004 


C .0404 


P0C5 


C.C052 


P0C6 


0 . 1660 


POlO 


0.2379 


poll 


C. 2383 


P012 


0.4717 


P013 


0.0250 


P015 


0.5607 


P016 


0 . 1340 


PC2C 


0.C920 


P02 1 


0.C097 


P022 


0.3576 


P053 


C.C905 


P027 


0. 1750 


P055 


0.1497 


P05C 


C .0337 


P05 1 


0 .4380 


P056 


0.0092 


P05 7 


C. 9540 


P024 


0.1581 


P054 


0.2860 


PC2 6 


C .0803 


POOl 


0.1750 


P030 


0.0404 


P045 


C. 0920 


PO? 1 


0 .0052 


P023 


0.0905 


P034 


0. 1660 


P003 


0.2860 


P04 6 


0.4717 


P022 


0.2379 


P035 


0.2383 


PC4 7 


0.1497 


P036 


0.0097 


P037 


0.0250 


P041 


0.1340 


P043 


0.5607 


P0^4 


0.3576 


PC07 


0.0337 


PO 14 


0.0092 


P042 


0.1511 


P017 


0.1581 


AOO 1 




0001 




A002 




A003 




0003 




A201 




A006 




A005 




0005 




AOlO 




A014 




AOU 




AC13 




0012 




AC17 




0017 




A022 




A202 




All 2 




A06 0 





8750E-04 0.7432E- 
6355'’-04 0 .53 15E-- 
3864E-04 0 .3286E- 
1053E-04 0.U91E- 



04 0.6279E-04 
04 0.44C5E-04 
04 0o2780E-04 
0''< 0.1313E-04 



0.3390E-04 0.45 
0.37C2E-04 0.30 
0.2390E-04 0.20 
0.L4C6E“04 0.15 
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AC52 

C053 

A206 

A2C0 

A05 I 

D051 

A057 

C057 

A202 

A035 

D035 

AO 2 6 

AC27 

C027 

A0A4 

A03A 

A033 

D033 

A03C 

A036 

A041 

DCA 1 

AOA3 

DOA3 

AO A2 

AlOA 

A052 

A007 

D007 

A205 

A20A 

A06 3 

C065 

A214 

A061 

C06 1 

A21 1 

A066 

A056 

A02A 

A216 

A213 

A207 

A 0 6 A 

A210 

A212 
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COMPUTER PROGRAM VII 



PEASE II. MANUAL CONTROL IS APPLIED 
JCB GALL225 
FORTRAN LSjGC 

REAL NMAX,NDIM,NSCAL,NBSCAL 
REAL NBAR,N,NBAR1 
DIMENSICN Y( 30,30) ,N( 30,30) 

DIMENSION POTADITO) , POT VAL ( 70 ) , AMPAD ( 70 ) ,AMPVAL(70) 
INTEGER POTADtAMPAC 
DATA N0/A8/ 

C STORAGE OF NOND I MENS I ON AL ARRAY FOR NAND Y 
DC 10 J=l,23 

10 READ! 5, IOC) ( V (J , I ) , I = 1 , 6 ) 

ICO FCPMAT(6E11.4) 

DC 12 J=l,25 

12 READI5, 100) IN(J, I ),I = l,6) 

X NRITE(6,200) 

200 FCPM.AT(» ','THIS IS YU, I)',//) 

DO 11J=1,25 

X WRITE (6,201) (YU, I), 1 = 1, 6) 

201 FCPMAT{‘ S10X,6(E11.4,2X) ,//) 

11 CONTINUE 

X WRITE(6,300) 

300 FGRMAT(‘ U'THIS IS N(J,I)',//) 

DC 13 J=l,25 

X WRITE(6,301) (N( J , I ) , 1 = 1, 6) 

3C1 FCRMAT( ' » , 1 OX , 6 ( E 1 1 . 4 , 2X ) , / / ) 

12 CONTINUE 

C SECTION CHECKING POT AND AMP VALUES FOR SHIP HI AND #2 
READ(5, 1005) (POT AD (I ) ,POTVAL( I ) ,1 = 1, NO) 

1005 FORMAT! A4,F10. 4) 

DC 105 1=1, NO 

105 CALL SETPOT( POTAD ( I ) , PCTVAL { I ) ) 

X WRITE(6, 1015) ( POTAD( I } ,POTVAL( I ) ,1 = 1 ,N0) 

1015 FCtRMATC A4,F10.A } 

READ(5, 1025) ( AMP AD ( I ) , I = 1 , 66 ) 

1025 FORMAT (A4) 

CALL RESETdOOO) 

DC 125 1=1,66 

125 CALL SCAN ( AMPAD( I ) , AMPVAL ( I n 
X WR1TE(6,103) (AMPACn ) ,AMPVAL(I ) ,1 = 1,66) 

102 FORMAT ( A4, FIO. 6) 

DC 20 1=1, NO 

2C CALL SC AN( P0TAD( I ) , POTVAL ( I ) ) 

X WRITE (6, 2005) (POTAD( I ) , POTVAL ( I ) ,1 = 1, NO) 

2005 FCRMATC A4,F10.6) 

C WHEN READY 

30 OUTPUT ( 101 ) 'READY ‘ 

C TRAN'SFORMA-^ ION OF COORDINATES FOR SHIP Hi AND S2 
115 CALL COMPUTE 
I FLAG = 0 

111 CALL ADK( 0,U1BAR, 1 , V1BAR,2 ,PSIBAR,4,U2BAR, 5,V2BAR,6, PS 
6R ,3, ALPBAR ) 

C RESCALING OF PSIBAR AND PS2BAR 
PSIMAX=0.26 
PS2MAX=0.26 
PSI=PSI MAX^PSI BAR 
A=SIN( P SI ) 

B = COS(PSI ) 

XC1D0T = U1 BAR*B-V IBAR’^A 
YO 1DOT = U1BA'5*A4-vibAR*B 
PS2=PS2MAX*PS2BAR 
A2=SIN(PS2) 

B2=C0S ( PS2 ) 

X02D0T=U2BAR*B2-V2BAR^A2 
YC2D0T=U2BAR^:A? + V2BAR*B2 

C INTERPOLATION IN 2-OlM. ARRAY. PASS PETBAR AND A1P3AR 
C P F- V C «rTPAp=l f) r ro pr c ppM nc opfA-IPA P crrTyr'it* 

C ALPBAP. = 1. 6“CCr<RESP0NDS' ALPHA = 550r6 fEET * 
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ALSr BETBAR AND ALPBAR ARE PAST IN CCPPUTER UMTS 
WRITE (6 ,401) 

401 FCRMAT(' ',10X,'THIS IS AL P B AR- E ET B AR ' , / / ) 

X WRITE(6,4021 AL P B AR , BE T 6 AR 

402 FCPMATI ' • ,10X,2FU.6) 

C RESCALING OF ALPHA3AR=ALPEAR AND BE T A E AR=B El B A R 
ALPHA=550. C*ALPBAR 
555 CCNTINUE 
X WRITE ( 6, 403) 

402 FCRMATI' ',10X,'THIS IS ALPHA',//) 

X WRITE! 6 ,501 ) ALPHA 

BETA= 100. 0-RETBAR 
X WRITE (6, 404) 

404 FORMAT!' *,10X,'THIS IS BETA',//) 

X WRITE! 6,501)8ETA 

I F I BETA .LT .50 .0 ) BETA=50.0 
I F ! BETA .GT .100.0 ) BETA = 100.0 
I F! ALPHA. GT .550.0) ALPHA = 550.0- 
I F! ALPHA. LT. -550 .0) ALPHA=-5 50. 0 
C SEARCI- FOR I NDET I F I C AT I CN OF ARRAY ELEMENT 

C FIND CCRRESPQNDFD Y AND N -LINEAR I NTERPOLAT I CN IN 2-DIMEN 
I E1--:BET A/10.0 
I B= IB 1-4 
EETAL=I BT410 
WTB=!BFTA-BETA1)/L0.0 
IF! ALPHA.! T.0.0) GO TO 1 
IF! ALPHA.GT.0.0) GO TO 2 
C ALPHA =0 
I A=13 
WTA=0.0 
GO TO 3 

C ALPHA LESS THAN ZERO 

1 I Al=-ALPHA/50 .0 
ALPHA1 = -I Al^=50 

WT A=-( ALPHA- ALPHA 1 )/50.0 
I A=13-I A1 
GO TO 3 

C ALPHA GREATER THAN ZERO 

2 IAl=ALPHA/50.0 
I A=134-I A1 
ALPHAl= I Al*50 

WTA = ( ALPHA-ALPHA 1 )/50.0 

2 1E2=IB+1 

IFIIB.EQ.6) IB2=6 
IA2=I A+1 

IF(IA.EQ.25) IA2=25 
X WRITE(6,L01) IBl , IB,IB2 ,I AL 

101 FORMAT! ' ' , lOX, ' IB 1= • , 13, 5X , ' IB = ' , I3,5X, ' IB2= ' , 13, 5X,' 
X WP1TE!6,102) IA,IA2 

1C2 FCRMAT!' ' , lOX , ' I A = ' , T 3 , 5 X , ‘ I A2 = ' , 1 3 , // ) 

X WRITE(6,202) B E TA 1 , WT E , ALPH A I 

2 02 FCRMAT ! • * , • B ETAl « , E14 .8 , 2X , ‘ WT B ' , E14 .8 , 2X , 'ALPHA 1 ‘ , El 

X WRITE!6,203) WTA 

202 FCRMA.TC* ' , • WT A= ' , E 14 . 8 , / / ) 

IF(ALPHA.LT.O.O) IA2=IA-1 
IF!IA.EG.l) IA2=1 
IF!WTA.EQ.0.0)IA2=IA 
IFIWTB.EO.0.0) IB2=IB 
NBAR1 = N ! I A , IB2) 

X WPITE!6,602) 

602 FORMAT!' ‘ ' THI S IS NBARl',//) 

X WRITE! 6,502) NEAR 1 , I A , I B2 

KEAR = N! I A, IB) 

X WRITE! 6,603) 

602 FORMAT!' ‘,‘THIS IS NEAR',//) 

X WPITE!6,502) KBAR,IA,IB 

BAPNl=!NBARl-NBARi *WTB+NBAR 
X WPITE!6,604) 

604 FORMAT!' ','THIS IS 8ARN1',//) 

X WRITE! 6, 501) BARNl 

NF ARI = N( I A?, IB2) 

X WRITE! 6, 60 5) 
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605 rORMAK* ',«THIS IS NBARl.*,//) 

X KRITE(6,502) NBARl , IA2 , I B2 

NBAR = N( IA2 , IB ) 

X k RITE (6, 6061 

606 FCPMAK' ','THIS IS NBAR’,//) 

X WRITE(6v502) NBAR,IA2,IB 

BARN2= (NBARL-NBAR) *KTE+NBAR 
X WRITE! 6 ,607) 

607 FORMAT!' ','THIS IS BARN2',//) 

X WRITE (6 ,501 ) BARN2 

NBAR=( BARN2-BARN1 ) *WT A-5- BA RN I 
X WRITE(6,500) 

500 FORMAT!* ','THIS IS FINAL NBAR*,//) 

X WRITE!6,501) NBAR 

501 FORMAT!' ' , 10 X , E 14 . 8 , / / ) 

VBARI=Y! lA, IB2) 

X WRITE!6,608) 

608 FORMAT!' ‘,'TH!S IS YBARl',//) 

X WRITE!6,502) YBARl, IA,IB2 

YEAR-Y! lA, IB) 

X W,PITE(6 ,605) 

605 FORMAT!* ','THIS IS YBAR',//) 

X WRITE!6,502) YBAR, lA, IB 

BAPYI=! YBARl- YBAR)* WTE+YBAR 
X WRITE (6, 610) 

610 FORMAT!' ','THIS IS BARYl',//) 

X WRITE!6,601)BARY1 

YBAR1--=Y I I A2, 182) 

X WPITE!6,611) 

611 FORMAT! * ',‘THIS IS YBARl',//) 

X WRITE! 6,502) YBAR I , I A 2 , I B 2 

YOAR=Y( IA2 , 1 B ) 

X WRITE! 6, 612) 

612 FORMAT!' ','THIS IS YBAR',//) 

X WRITE! 6, 502) YBAR, I A2, IB 

BARY2-! YBARl- YBAR) *WTB+ YBAR 
X WRITE!6,613) 

612 FORMAT!' ','THIS IS BARY2',//) 

X WRITE!6,601) BARY2 

YBAR=! BARY2-BARY1 ) ^'WTA + BARYl 
X WRITE(6,60C) 

60C FORMAT!' ','THIS IS FINAL YBAR',//) 

X WR1TE(6,601) YBAR 

601 FORMAT! • ' , 1 0 X , E 1 4 . 2 , / / ) 

502 FORMAT!' ‘ , lOX , E 14 . 8 , 5X , ' ROW= ' , I 3 , IX , ‘ COL = ' , I 3 ) 

C REMEMBER SIMULATION IS BONE IN COMPUTER UNITS SC NO NEEC T 
C**-** DIVIDE BY 100.0 SINCE SUBR. DAC CONVERTS BY ITSELF IN 
X WRITE(6 ,55555) IFLAG 

55555 FORMAT!' ' , 5X , ' I F L AG = ' , I 1 , / / ) 

C MAGNITUDE SCALING OF N 
NMAX = 4. 195E-^07 
NDIM=NBAR*(9.428E+1C) 

IF! IFLAG. EQ.O) NSCAL = NDIM/NMAX 
IF!IFLAG.EC.l ) NBS C AL=NDIM/NMAX 
C MAGNITUDE SCALING OF Y 
YMAX=1 . 608E+05 
YCI M=YBAK* ! 17 .87E+07 ) 

IF! IFLAG. EQ.O) YSC AL = YD IM/ YMAX 
I F! IFLAG .EQ. 1 ) Y bS C AL = YD I M/ YMAX 
X WRITE!6,501) NDIM,YDIM 

901 FORMAT!' ' , ' NDI M- ‘ , E 14 . 6 , 5X , ' YD I M= • , E 14 . 6 , / / ) 

X WR ITE ! 6 ,5001 ) X 0 1 D CT , YO 1 DOT 

5001 FORMAT!' ' , ‘ XO 1 DOT = ' , F 14 . 6 , ' YOl DGT= ■ , E 1 4 . 6 , // ) 

X NRITE!6,5000) X02 DOT , YO 2D0T 

5000 FORMAT!' ‘ , ' X02 DOT = ' , E 14 . 6 , ' Y02 DOT = ' , E 1 4 . 6 , // ) 

IF! IFLAG. EQ.O) I F L AG= 1 ; AL PH A=-A L Ph A ; GO TO 555 
X WRITE!6,900) NSCAL , YSCAL , NBSCAL , YBSCAL 

90C FORMAT!' • , • NSC AL = ' , E 1 4 . 6 , 2X , ' Y SC AL= ' , E I 4 . 6 , 2X , ‘ N B S C AL 
6'YBSCAL-' ,E14.6,// ) 

8 04 call DAC! 1 -XOIDOT , 2 , Y 0 1 DOT , 5 , YO 2 COT , 6 , XO 2 CO T , 3 , N SC A L ,4 
CCLjetYPSCAL) 

C MANUAL DIGITAL SWITCH'D' AND '1' 
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15 IFUESTd ) .GT.O) 
CALL HOLD 

L7 I RTESr ( 1 ) .LT .0 ) 
CALL COMPUTE 
IC IF (TEST (3) .GT .0) 
CALL RESET(IOOO) 
PAUSE 
GC TO 115 
EKD 

-.LOAD XR,MAP 
-.DATA 
❖ 

A 

A 

A 

* 

A 
A 

-0. lACOF- 
-0 .24C0E- 
-0. 3200E- 
-0.3800E- 
-0. ACCOE- 
-0 . 3500E- 
-0. 26CCE- 
-0. 13C0E- 
O.AOCCE- 
0.27C0E- 
0.5C00F- 
0.70 00E- 
0. 8 5CCE■ 

C.9C0QE- 
0r82C0E- 



GC TO 16 
GO TO 17 

IFLAG = 0;GO TO LU 



THIS IS Y ARRAY 



6250E- 
A5C0E- 
3000E- 
17CCE- 
60CCE- 
-C. 2CCCE- 
-0 .EDOO"-'- 
6CC0E- 
5C00E- 
ICOOE- 



0 

0 

0 

C 

0 



-0 

-0 

-0 

A 

A 

* 



03-0 

03-0 

03-0 

03-0 

03-0 

03-0 

03-0 

03- C 
04 0 
03 0 
03 0 
03 0 
03 0 
03 0 
03 0 
03 0 
03 0 
•03 0 
03 0 
•04 0 

04- 0 
04-0 
04-0 
04-0 
•04-0 
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